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PREFACE 


A sound mathematical training for the modern engineer 
is a sine qua non for attaining new heights in all aspects 
of engineering practice. One of the areas of mathematics 
that plays a big role in the mathematical education of 
the engineer is vector analysis, which is now invariably 
included in the curriculum of higher mathematics in 
engineering colleges. 

The present collection of problems in vector analysis 
contains the required minimum of problems and exercises 
for the course of vector analysis of engineering colleges. 

Each section starts with a brief review of theory and 
detailed solutions of a sufficient number of typical prob¬ 
lems. The text contains 100 worked problems and there 
are 314 problems left to the student. There are also a 
certain number of problems of an applied nature that 
have been chosen so that their analysis does not require 
supplementary information in specialized fields. The 
material of the sixth chapter is devoted to curvilinear 
coordinates and the basic operations of vector analysis 
in curvilinear coordinates. Its purpose is to give the reader 
at least a few problems to develop the necessary skills. 

The exposition in this text follows closely the lines 
currently employed at the chair of higher mathematics 
of the Moscow Power Institute. 

The present text may be regarded as a short course in 
vector analysis in which the basic facts are given without 
proof but with illustrative examples of a practical nature. 
Hence this problem book may be used in a recapitulation 
of the essentials of vector analysis or as a text for readers 
who wish merely to master the techniques of vector 
analysis, while dispensing with the proofs of propositions 
and theorems. 

In compiling this problem book, the authors made 
extensive use of material in published courses of vector 
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calculus and collections of problems. Many problems 
were made up by the authors themselves. 

This collection of problems is designed for students of 
day and evening departments at engineering colleges and 
also for correspondence students with a background of 
vector algebra and calculus as given in the first two years 
of college study. 

We would like to express our sincere gratitude to 
Professor V. P. Gromov (the Moscow Krupskaya Pedagog¬ 
ical Institute), Professor A. V. Efimov and Associate 
Professors I. M. Petrov, B. I. Fridlender, and V. N. Zem¬ 
skov (the Institute of Electronics) for their thorough 
scrutiny of the manuscript of the book and for valuable 
tuggestions and remarks that were made full use of in 
she final editing. 


Moscow-Dubna, 1977. 


M. L. Krasnov 
A. I. Kiselev 
G. I. Makarenko 



CHAPTER I 


THE VECTOR FUNCTION 
OF A SCALAR ARGUMENT 


Sec. 1. The hodograph of a vector 
function 

Definition 1. A vector r is said to be the vector function 
of a scalar argument t if each value of the scalar taken 
from the domain of admissible values is associated with 
a definite value of the vector r. This can be written as 
follows: 

r = r ( t ). 

If the vector r is a function of the scalar argument t, 
r = r (t), 

then the coordinates x, y, z of the vector r are also func¬ 
tions of t: 

x = x (f), y = y (*), z = z (t). 

Conversely, if the coordinates of the vector r are func¬ 
tions of t, then the vector r itself is also a function of t: 
r = x (t) 1 + y (t) j + z (t) k. 

Thus, specifying a vector function r ( t) is the same as 
specifying three scalar functions x (t), y (t), z ( t ). 

Definition 2. The hodograph of the vector function 
r (t) of a scalar argument is the locus described by the 
terminus of the vector r (t), as the scalar t varies, when 
the origin of the vector r (t) is fixed at a point 0 in space 
(Fig. 1). 

The hodograph of a radius vector r = r (0 of a moving 
point is the trajectory L of that point. Some other line L x 
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(Fig. 2) is the hodograph of the velocity v = v (/) of that 
point. Thus, if a material point (particle) is in motion 
around a circle with constant ve¬ 
locity, | v | = constant, then its 
hodograph of velocities is like¬ 
wise a circle with centre at 0 1 
and with radius equal to | v |. 

Example 1. Construct the ho¬ 
dograph of the vector r = ti+ 
t\ + t 2 k. 

Solution. 1°. This construc¬ 
tion may be carried out by using 
points and setting up a table: 



Fig. 1 


/ 

1 0 

i 

2 

3 

> 

3 

0 1 

i + j + k 

2i-H2j + ik 

31 + 3j+9k 

4i + 4j + 16k 


2°. Alternative solution. Denote by x, y, z the coordi¬ 
nates of vector r; we have 

x = t, y = t, z = t 2 . 


Eliminating the parameter t from these equations, we 
get equations of the surfaces y = x, z = x 2 , the line L 



of intersection of which is what defines the hodograph of 
the vector r ( t) (Fig. 3). 
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1. Construct the hodographs of the following vectors: 
(a) r = 2i + *2j — Pk. 


(c) r = cos t-i -f-sin t- j + k. 



Fig. 3 

(d) r = /i + 4< 2 i + 4-* 3k - 

/-x _ 2(i + 2tj + (<*— 2) k 

W r --J* + 2-' 

Sec. 2. The limit and continuity 
of a vector function of a scalar argument 

Suppose a vector function r = r (t) of a scalar argument 
t is defined in some neighbourhood of the value t 0 of the 
argument t, except perhaps for the value t 0 itself. 

Definition 1. A constant vector A is said to be the 
limit of the vector r (f), as t -+■ t 0 , if for any e > 0 there 
is a 6 >0 such that for all t ^ t 0 that satisfy the con¬ 
dition | t — t 0 | <C 6 the following inequality holds true: 

| r (t) — A | < e. 


12 


The Vector Function 0 / a Scalar Argument [Ch. 1 


As in the case of ordinary calculus, we write lim r ( t) = 
= A. 

Geometrically (Fig. 4), this means that the vector 
r (t) tends, as t -*■ t 0 , to the vector A both in length and 



satisfy the condition 1 1 — 
| ct (/) | < e holds true. 

Example 1. Show that the 
infinitesimal when t -*■ 0. 
Solution. We have 


in direction. 

Definition 2. A vector 
a ( t) is said to be infini¬ 
tesimal, as t t 0 , if 
a(t) has a limit, as t-*-t 0 , 
and that limit is zero: 

lima(0 =0, 

0 

or, what is the same 
thing, if for any e >0 
there exists a 6 >0 such 
that for all t =£ t 0 that 
| < 6, the inequality 

vector a ( t) = tl + sin t\ is 


|o(0| = |« + sin *j|<|t| + |sin*|<2|*|. 

From this it is evident that if for every e >0 we take 
6 = e/2, then for | t — 0 | < 6 = e/2 we have | a (t) \ < 
< e. By the definition, this means that a ( t) is an infini¬ 
tesimal vector when t-*~ 0. 

2. Show that the limit of the modulus of a vector is 
equal to the modulus of its limit (if the limit exists). 

3. Demonstrate that for a vector function r (t) to have 
a limit A, as t 1 0 , it is necessary and sufficient that 
r (0 be representable in the form 

r(t) = A +a (t), 

where a ( t) is an infinitesimal vector when t 1 0 . 

4. Show that if the vector functions a (t) and b (t) have 
limits as t 1 0 , 

lima(*) = A, limb(0 = B, 

t-t a «-t 0 
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then their sum a (t) + b ( t ) and their difference a ( t) — 
— b (f) also have limits as t -*■ f 0 , and 
lim [a (t) ± b (f)] = A ± B. 

t-( 0 

5. Let 

lima(t) = A, limb(t) = B. 

(-<0 «-< 0 

Prove that 

lim (a (t) t b(t))-(A f B), 

«-►<# 

where (a ( t ), b (f)) is a scalar product of the vector func¬ 
tions a (t) and b (f). 

6. Let 

r (t) = x (t) I + y (t) j + z ( t) k, A = aji + aj + 

+ a,k. 


Show that if lim r(<) = A, then 
<-<« 

lim x(t) = a lt limi/(f) = a 2 , lim z(t)=a 3 . 

*—<o ‘-‘o <-< 0 

Find the following limits: 

7 . ls (^ 1+ -i=i j+e . k ). 

8 ' ““ (T^ i,+ TW i + k )- 
9 - 1 ,‘“( i T i|+ ' :os ‘- i+ iTir)- 
10. lto(^£ii + i±221ii + i-k). 

“• 1 l S(7=r l + T=T i + 2k )- 

Definition 3. A vector function r = r {t) defined in 
some neighbourhood of the value t = t 0 is said to be 
continuous when t = t 0 if 

limr(0 = r(f 0 ). 
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In other words, r = r (t) is continuous for t = t 0 if for 
every e > 0 there is a 6 > 0 such that for all t that 
satisfy the condition 1 1 — t 0 | < 6 the inequality 
| r (i) — r (t 0 ) | < e holds true. 

The hodograph of a continuous vector function of a 
scalar argument is a continuous curve. 

12. Start with the familiar inequality |a — b | > 
;> | | a | — | b | | and demonstrate that the continuity 
of a vector function implies the continuity of its modulus. 
Is the converse true? 

13. Show that if a ( t ) and b (t) are continuous for t = t 0 , 
then the vector function a (t) ± b ( t ) is also continuous 
for t = < 0 . 

14. A vector function a (t) + b ( t) is continuous for 
t = t 0 . Does it follow from this that the vectors a (t) 
and b ( t) are also continuous when t = t 0 1 

15. Prove that if a (t) and b ( t) are continuous vector 
functions, then their scalar product (a (0, b (*)) and 
vector product [a ( t ), b (t)l are also continuous. 


Sec. 3. The derivative of a vector 
function with respect to 
a scalar argument 

Suppose a vector function r = r (t) is defined for all t 
on the interval (t 0 , t x ). Take some value t £ (£„, tj, then 
give t an increment At such that t + At 6 (*o» h) an d find 
the corresponding increment Ar = r (t + At) — r ( t) in 
the vector function r ( t ). Now consider the ratio At/ At. 

Definition. If, as At 0, the ratio At/A t has a limit, 
then that limit is called the derivative of the vector 
function r = r (t) with respect to the scalar argument t 
for a given value t of the argument and is denoted as 
dr ( t)/dt or r' ( t ) or r (/). Thus, 


lim-g-lim r(t+A A ?~ r(<) . 

* A<-.0 A< At-0 


In this case the vector function i 
differentiable. 


= r ( t) is said to be 



Sec. 3] 


The derivative of a vector function 


15 


16. Show that if the vector function r = r (t) has a 
derivative for some value t of the argument, then it is 
continuous for that value t. 

The derivative of a vector function r ( t) of a scalar 
argument t is a vector directed along the tangent to the 



hodograph of the original vector at the point under con¬ 
sideration (Fig. 5). The vector dr/dt is in the direction 
of the terminus of the vector r (<) as it moves along the 
hodograph when the parameter t increases. 

Suppose r = r (t) is the radius vector of a moving 
point. Then the vector v = dr/dt is the velocity vector 
of that point. 

Suppose 

r (t) = x (t) i + y ( t ) j + z (!) k, 

where the functions z (0, y (0. z (0 are differentiable at 
the point t. Then there exists dr/dt for that value of t and 


dr _ dx . , dy_ 

~dt~~dt X+ ~dt 


k. 


(1) 


Example 1. Find dr/dt if r = la cos t + \b sin t (the 
point is moving in an ellipse). 

Solution. From formula (1), 


■4r = — iasin* + jhcosf. 
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By analogy with the differential of a scalar function, 
the differential of a vector function r — r (t) is a vector dr 
defined by the equality 

where dt = At is the increment in the scalar argument t. 

As in the case of scalar functions, 

Ar = dr + a- At, 
where a = a (t, At) -*■ 0 as At -*■ 0. 

Basic rules for differentiating a vector function 

Assume that all functions being considered (both scalar 
and vector) are continuous and differentiable. 

1°. If c is a constant vector, then dc/dt = 0. 

2°. The derivative of a sum of vector functions is equal 
to the sum of the derivatives of the summands: 

<* (a «)+!>«)) da , dh 

dt dt dt 

3°. Suppose a vector function a (t) is multiplied by 
a scalar function m ( t ) of the same scalar argument. Then 
d{m») da dm 

-di-= m 7r+nr a 



(In this formula, the order of the factors a and b in the 
right-hand member must be the same as that in the left- 
hand member.) 

Let us prove formula 4°. We set <p (t) = (a (t), b (t)). 
Give t an increment At; then, by the distributive prop¬ 
erty, we have for the scalar product 
A<p = <p (f + At) — <p (t) = (a + Aa, b + Ab) — (a, b) 
= (Aa, b) + (a, Ab) + (Aa, Ab), 

whence 
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It is given that the functions a (t) and b (t) have deriva¬ 
tives for the value t of the argument and, hence, are 
continuous for that value of t. Therefore 

lim, lim-^- = -£-, and lim Ab = 0. 

A(-»0 A< dt AI-0 A< dt Af-*0 

Passing to the limit in (2) as At 0, we obtain 



17. Given r = r (t). Find the derivatives: 

( a ) ~7t (r 2 )* < b >-£:( r -ir)• «»£[>£]• 

18. Prove that if the modulus | r | of the vector function 
r = r (t) remains constant for all values of t, then dr/dt J_r. 
What is the geometrical meaning of this fact? 

Prove that if e is a unit vector in the direction of 
the vector E, then 

(e , del = 1*^1. 

20. Suppose 

u = Hi (*. J/. 2. 0 » + u 2 {x, y, z,t) j + u 3 (x, y, z, t) k, 
where u v u s , u 3 are continuously differentiable functions 
of their arguments, and x , y, z are continuously differentia¬ 
ble functions of t. Show t hat 

du du du dx , gu dy . gu dz 

~sr dt ' 6x ir + ~dir~dr + ~dr~dr- 

21. Find the trajectory of motion for which the radius 
vector r ( t) of a moving point satisfies the condition 
dr/dt = [a, r], where a is a constant vector. 

The derivative dr/dt of the vector function r (t) of a 
scalar argument is a vector function of the same argument. 
If there exists a derivative of dr/dt, then it is called the 
second derivative and is denoted dh/dt *. Generally, 



22. Given the radius vector of a point moving in space: 
r{osint, —acost, bt 2 ) 

(t is time, and a and b are constants). Find the hodographs 
of velocity and acceleration. 


!-910 
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23. Given: r = a cos ait + b sin a>t, where to, a, b are 
constants. Prove that 

(,) [ r ’ 

< 2 >Jr+“ ir =°- 

24. Show that if r = ae“‘ + be““ ( , where a and b are 
constant vectors, then d?r/dt — co*r = 0. 

■25. Show that the modulus of the differential of the 
radius vector of a point is equal to the differential of 
the length of the arc described by the point. 

26. Suppose a = a (u) is a vector function of a scalar u, 
where u in turn is a certain scalar function of the basic 
scalar t. Assuming a (u) and u = u (t) to be differentiable 
the necessary number of times, find an expression for the 
derivatives of the composite function daldt, cPa/dt 2 . 

Sec. 4. Integrating a vector function 
of a scalar argument 

Definition 1. We will say that the vector function A ( t) 
is the primitive of the vector function a ( t) when t 0 <t< 
<. if A (t) is differentiable and 

-^T = a(0. *€(* o.M- 

Definition 2. The collection of all primitive functions 
of a (t) is termed the indefinite integral of -the vector 
function of a scalar argument a = a (t). As in integral 
calculus, the indefinite integral of a vector function is 
denoted by the symbol j, and we have 

j a(<)d* = A(t) + C, 

where A ( t ) is one of the primitive functions of a ( t ), and 
C is an arbitrary constant vector. 

The following properties hold true for integrals of 
vector functions: 

1°. jaa(/)<fr = aja(0<fr (a is a numerical constant). 
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Integrating a vector function 


2°. J (a(£) + b(£))d£ = J«(0*± Jb(0*. 

27. Show that if c is a constant vector, and a (I) is a 
variable vector, then 

j (c,a(£))d£ = (c, j a (t)dtj , 

j [c,a(£)]d£ = [c, j a(£)d£]. 


a (t) = fli (t) i + a t (t)j + a 3 (t) k, 

then 


j a(£)d£ = i j a^dt + j j a 2 (t)dt + k j o,(£) dt. (1) 

That is, the integration of a vector function reduces to 
three ordinary integrations. 

Example 1. Find the indefinite integral for the vector 
function a (£) — i cos t + + k. 

Solution. According to formula (1), 

j a (£) dt = i j cos t dt 4- j j <r‘d£ -f k j dt = 

= isin£ — je"‘ + k£ + c, 


where c is an arbitrary constant vector. 

Find the integrals of the following vector functions: 

28. a(£) = £e‘i + sin 2 £-j —-q^g-. 

29- a(t) = T | ir + te t, j + cost-k. 

30. a(£) = cos£e 8,n< -i — £cosi 2 -j + k. 


31. a (£) = -£-£ 2 i —£sin £-j + 2‘k. 

Let a vector function a (£) be defined and continuous 
over a certain interval [£„, T\> which is the range of the 
argument £. 

Definition 3. We define the definite integral of a vector 
function a (£) on the interval [£ 0 , T ) as the limit of the 
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vector integral sums 

o= S a(T h )At fc , t*6[*».<*+! 1 

as the length A t of the largest of the subintervals [f k , f k+1 ] 
(k = 0, 1, .... n — 1) into which the interval U 0 , T\ 
is partitioned tends to zero: 

X n—1 

\ a(<) dt = lim 2 »(**) Af*. 

A, -° ft-0 

The following formula holds true: 

T 

j a(t)dt = A (T) — A (t 0 ), (2) 

to 

where A (t) is some primitive for the function a ( t) on the 
interval U 0 , T 1. 

If 

a (t) = a, (0 1 + a, (t) j + a, ( t) k, 

then 

TXT X 

j a (0 dt = i j a, (t) dt + J j a 2 (t) dt + k j a 3 (0 dt. (3) 

lo <0 to to 

n/2 

Example 2. Compute j a (t)dt, where a(0 —icos/ — 

o 

— jsin 2 f. 

Solution. By virtue of formula (3), 

n/2 n/2 n/2 

j a(f)dt = ij cosfdf — jj sin 2 *df 

= isin *|o -MT-4 )|o =*-Ti- 

Compute the following integrals: 

32. fa {t)dt, where a = sin z * cos <-i + cos 2 <sinfj + k. 
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33. j a (t)dt, where a = — y +-^— + ke*. 
o 

34. j a (t)dt, where a = 3ncosn<-i—+ 2*k. 
o 

n 

35. j fk(t)dt, where a = (2f + n) i + *sinf -j + nk. 

Example 3. An electric current I flows upwards along 
an infinite wire that coincides with the z-axis. Find the 
vector H of the magnetic 
field intensity set up by 
this current at an arbitrary 
point M ( x, y, z) of space 
(Fig. 6). 

Solution. We consider a 
sufficiently small element 
PPx = d£ of the z-axis. By 
the Biot-Savart law, the 
intensity dH of the mag¬ 
netic field set up at point 
M by the current flowing 
through element d£ of the 
wire coincides in direction 
with the vector product 
[d£, rj, where d£ = PP^ 

|d£| = d£, r, = PM (see c 

Fig. 6). By this same law, the Fl 8- 6 

modulus of the vector dH is 

l dH l=-^-sin(d£, N r I )dC, 

where (d£, r x ) is the angle formed by the vectors d£ and r x . 

Since 

I[<*5. M| =r,dCsin (d£,"r,), 

dH = -^-fdg,r 1 ). (4) 



we can write 
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In order to obtain the desired vector H at the point AT 
we have to sum all vectors dH pertaining to distinct 
elements PP x of the wire, that is, we have to integrate 
the expression (4) over the whole z-axis: 


h-Tt'*-' 1 '- 

We have 

r x = OM - OP. 

But 

OM = xi + yi + zk, OP = £k. 


and therefore 


(5) 


so that 


= xl + yj + (z - 0 k 


r .-|ri|-V'*» + y , + (i-0 1 -V'p*+(*—0*. 

where p = ]/ x 1 + y 2 is the distance of point M from the 
axis of the wire. 

For the vector product [d£, rj we have 


K, r,| 


i j k 

0 0 d£ 

x y z—i 


= — ipdC +jxd£, 


and formula (5) takes the form [the point M ( x , y, z) is 
fixed, I = constant] 


H = /(—yi + xi)] [pl+(i 1 0 . (6) 

To compute the integral on the right-hand member of (6), 
make the substitution 
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We then have 

+oo n/2 

f _ K __ f P dt 

J [p* + (z—£)*] 3/t J cos* < [p*+p* tan* t] 3/1 
—n/2 

n/2 

-n/2 

Thus, the intensity vector H of the magnetic field is 
in our case given by 

H = —yi 4-*j) 

or 



where I = 7-k is the current vector, r is the radius vector 
of point M (x, y, z) of the field, and p is the distance of 
M to the axis of the wire. 

Example 4. The motion of an electron in a homogeneous 
magnetic field. 

1°. Suppose a magnetic field H is set up in some region 
of space; let it be constant in magnitude and direction 
(a homogeneous field). Suppose at time t = t 0 , an elec¬ 
tron enters the field with an initial velocity v 0 . Determine 
the path the electron will take. 

Solution. First suppose the vector v 0 is perpendicular 
to H and that the initial position of the electron is at 
point M 0 . Choose the origin 0 at an arbitrary point of the 
plane P passing through M 0 at right angles to the vector H 
(Fig. 7). Let the initial radius vector OM 0 be r 0 , let r be 
the radius vector of the electron at the current instant of 
time t, and let v be the instantaneous velocity at that 
instant. The basic differential equation of motion is 



It will be recalled that the force F acting at time t on the 
electron by the magnetic field is 

F = — e 0 [H, vl, 
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where e 0 is the absolute value of the electron charge. 
Thus, 

m -S- = e o[ v ’ H l- (?) 

At every instant t, the force F is perpendicular to the 
direction of velocity and to the direction of the held H; 



at every instant, it will force the electron to deviate from 
a rectilinear path and to describe a certain curvilinear 
trajectory. 

Let us rewrite (7) as 



and integrate from t 0 to t with respect to t. This yields 

mv — mv 0 = e 0 [r, H) — e a [r 0 , H] 
or 

mv = e 0 [r, H] + (mv 0 - e 0 [r 0 . Hi). (8) 

Now choose a coordinate origin O' such that the term in 
parentheses in the right-hand member of (8) vanishes, that 
is, so that 


e 0 [r 0 , Hi = mv,,. 
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From (9) it follows that the initial vector r 0 must be 
perpendicular to the vector v 0 and must lie on the straight 
line M 0 K, which is perpendicular to the plane of the 
vectors v 0 and H. By virtue of (9), the modulus of the 
vector r 0 must satisfy the relation 

e 0 |r 0 I-1 H | = m | v 0 |, 

whence 

l r °l = e 0 |H| • 

This determines the position of the new origin O'. Relative 
to this origin, equation (8) is rewritten thus: 

mv = e 0 [r, H] (11) 

or 

m£ = e 0 [ r.H], (12) 

From equation (11) it follows that the trajectory of the 
electron is a plane curve lying in the plane P because at 
every instant the vector v is perpendicular to H. Now take 
the scalar product of both sides of (12) by r: 

m (r, -^-)=e 0 ( r < fr, H)). (13) 

The mixed product in the right-hand member of (13) 
is zero, so that 

('•#)-«• 

whence 

-^-(r 2 ) = 0 or -^-(r 2 ) = 0, that is, r 2 = constant. 


This is the equation of a circle lying in the plane P with 
centre at the chosen point O'. The radius of the circle is 
found from formula (10) since the initial point Af 0 must 
also lie on that circle. Thus, we finally have 

r = r «=-5Twf- (10) 


Thus, if an electron enters a homogeneous magnetic 
field H with an initial velocity v„ at right angles to H, 
then it will describe, in that field, a circular trajectory 
lying in the plane P perpendicular to H and passing 
through the initial point. The radius of the circle is 
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given by formula (10) and its centre O' lies on the straight 
line perpendicular to the plane of the vectors v 0 and H; 
note that a rotation from v 0 to H must be seen from point 
O' as a counterclockwise rotation. 

From (10) it is evident that the radius r 0 of the circle 
is inversely proportional to | H |. Thus, the greater the 
intensity of the magnetic field, the greater the curvature 
of the trajectory. 

From formula (11), 

my = e 0 [r, H], 

it is clear that if r is constant in modulus and is all the 
itme perpendicular to H, then also the velocity v of the 



point will be constant in magnitude, 

| v | = v 0 = constant, 

so that the electron is in uniform motion in the orbit. 
The period of revolution T is 


T 


(14) 


This formula does not involve the initial velocity v 0 . 
Thus, irrespective of the initial velocity v 0 which is 
perpendicular to H with which the electron enters the 
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homogeneous magnetic field H, it will perform a single 
orbital revolution and always in the same time T. 

2°. Now suppose an electron enters a homogeneous 
magnetic field H with some initial velocity V that is 
not perpendicular to the vector H. This velocity may 
then be resolved into two components: the vector v 0 
at right angles to the field, and a vector v x parallel to the 
magnetic field. 

From the formula 

F = a (V, HI = e 0 [v 0 , Hi 

it is evident that the “twisting” force F is given only by 
the perpendicular component v 0 and that it imparts to 
the electron a rotational motion about the circle (centred 
at O') discussed above. As for the other component v J , 
the electron will retain it by inertia and, besides having 
a uniform circular motion, it will have a rectilinear and 
uniform motion in the direction of H with a velocity 
v t = | V | cos a. The combination of these motions 
yields a helical curve with axis parallel to the vector H 
and passing through the point O' (Fig. 8). 

Sec. 5. The first and second derivatives 
of a vector with respect to the arc length 
of a curve. The curvature of a curve. 

The principal normal 

Consider a curve L in space. On it, choose a point M 0 
as the origin and also choose a direction along L that will 
be regarded as positive. For a parameter, take the arc 
length s reckoned from M 0 of the curve ,(Fig. 9). Then 
the radius vector of a point M of the curve is 

r = r (s). 

With that choice of parameter, 



where t° is a unit vector directed along the tangent to 
the curve L in the direction of increasing values of the 
parameter s. 
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If the vector r is given by the coordinates 



r = xi + y] + zk, 

then 

’•=£‘+li+£ k 

and 




Since | r° 
vector t°. 

| = 1, the vector dx°/ds is orthogonal to the 


The modulus of the vector dx°lds is 


Here, K is the curvature of the curve L at the point M. 

The straight line having the direction of the vector 
dx°/ds and passing through the point M of the curve is 



termed the principal normal of the curve at the point M. 
Denoting the unit vector of that direction by n°, we have 


dxQ 

ds 


K n°. 


( 1 ) 
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The inverse of the curvature of a curve at a given point 
is called the radius of curvature of the curve at that point 
and is denoted by R: 



Thus, formula (1) may be rewritten as 

dh __ dx° _ n° 

ds* ds ~ R * 

From this, 



Using (2), we can compute the curvature of a curve at any 
point if the curve is specified by parametric equations in 
which the parameter is the arc length s. 

In the particular case of a plane curve, a circle of 
radius a, 

x = a cos i- t 
y = a sin—, 

we have 



and formula (2) yields 

X = -jj- = |/^-jj- cos 2 — -|- —f sin 2 — — —. 

This means that the curvature of a circle of radius a is 
constant and is equal to the inverse of the radius of the 
circle. 
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If the curve L is given by the vector-parametric equa¬ 
tion r = r ( t ), where the parameter t is arbitrary, then 


K = 


i*r 


( 3 ) 


Formula (3) permits computing the curvature of the 
curve at any point prov ided we have an arbitrary para¬ 
metric specification of that curve. 

Example 1. Compute the curvature of the helical curve 
r — a cos t-i + a sin i* j + Atk. 

Solution. Since 


-^-= — asin<-I + acosf.j + Ak, 

jlf 

— acosM—asinf>j, 

the vector product 

I j k 

- a sin t a cos t h 

- a cos t — a sin t 0 

= aA sin <-i — ah cos t • j + a 2 k. 

Consequently, 


r — *** l 

i dt • di> 


By virtue of (3), 


or 

R = gf+fcl = constant. 

a 

Thus, a helical curve has a constant radius of curvature. 
Find the radius of curvature of each of the given curves: 

36. r = In cos M + In sin t -j + V 2t-k. 

37. r = t 2 i + 2f*j. 
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38. r = 3t a l + (3 1 - t*) j + 2k for t = 1. 

39. r = a (cos t + t sin t) i + a (sin t — t cos t) j 
for t = n/2. 

40. r = a cosh M + a sinh t -j + atk at any point t. 


Sec. 6. Osculating plane. Binormal. 

Torsion. The Frenet formulas 

The plane passing through the tangent line and principal 
normal to a given curve IT at a point M is termed the 
osculating plane at the point M. 

For a plane curve, the osculating plane coincides with 
the plane of the curve. 

If the vector r = r (t) has a continuous derivative 
drldt in the neighbourhood of a point t 0 and, besides, a 
Second derivative (Pr (t 0 )/dt 2 such that 

then at the point t =* t 0 there is an osculating plane ‘o 
the curve r = r (t) whose vector equation is 

where p = p (t) is the radius vector of the current point 
of the plane. 

The normal to the curve at the point M, which normal 
is perpendicular to the osculating plane of the curve at 
that point, is called the binormal of the curve at the given 
point M. 

Denote by b° the unit vector of the binormal oriented 
so that the vectors x°, n°, b° form a right-handed trihedral 
(Fig. 10). Then 

b°’ = 1, b° = [T°, n°J. 

For the derivative db°/ds we get 



The vector db°/ds is perpendicular both to the vector t° 
and to the vector b^, that is, it is collinear with the vec- 
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tor n°. Set 

We then have 


The quantity i/T is termed the torsion of the given 
curve, and T is called the radius of torsion of the curve. 

The torsion of a curve 
is given by the formula 

_L_ m(*L *!L\ 

where the symbol (a, b, c) 
denotes a mixed product of 
the vectors a, b, c, that is, 
(a, b, c) = (a, [b, cl). 

For the case where the 
curve is given by the vec¬ 
tor-parametric equation r = 
= r ( t ), we have 



/ dt d 1 r d? i \ 




(1) 


Example 1. Find the torsion of the helical curve 
r = a cos M + a sin t -j + htk. 

Solution. We find the derivatives of the given vector: 
^f= — asinM-f acost-j + Ak, 


— acost-i — asinf-j, 

d*r ... . • 

- 5 - 3 - = a sin / •! — a cost-]. 
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The mixed product of these vectors is 

—asin£ a cost 

/ dt d*r d 3 t \ . . 

(dT'-ST.-sr ) = — cos* a sin t 

a sin* —a cost 

In example 1, Sec. 5, we found that 

it*. 4wr—^+*'»- 

Using (1), we obtain for the torsion 
1 h 
T ~ «»+A»* 


Thus, the torsion of a helical curve is the same at all 
its points. 

Example 2. Write the equation of the osculating plane 
at the point t = 0 of the helical curve 

r = a cos M + a sin t -j + htk. 

Solution. We find the values of the derivatives of the 
given vector and its derivatives dr/dt and cPrldt 2 at the 
point t = 0: 

r(0) = ai, dT Jp- = ai + hk, d *J t T " = ~ ai - 
Consequently (see example 1, Sec. 5), 

KP-. i+a 

The vector equation of the osculating plane is 
(p-r(0). Jjp*. ^)-0 

or 

(p — ai, — ah\ + a 2 k) = 0. 

Since the radius vector of the current point of the osculat¬ 
ing plane p = xi + y\ + zk, it follows that by passing 
to coordinate notation we obtain an equation of the 
desired plane in the form hy — az = 0. 

Formulas expressing the derivatives of the vectors 
t°, b°, n° are called Frenet formulas-. 


dt° 1 0 db° _ 1 0 
ds ~ R " ’ ds T " ’ 


dn° 

ds 


-5- 1° — -5T b°. 


3-910 
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41. Write down the equation of the osculating plane at 
the point t = 2 of the curve 

r = <i-*j + y* 2k - 

42. Write down the equation of the osculating plane at 
the point t = 0 of the curve 

r = e‘i + e-‘j + /2<k. 

43. Find the torsion at the point t = 0 of the curve 

r = e l cos t-i + e* sin t-j + e'k. 

44. Find the torsion at any point t of the curve 

r = a cosh + a sinh t- j + atk. 
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Sec. 7. Examples of scalar fields. 

Level surfaces and level lines 

Definition. If a value of a certain quantity is defined 
at every point of space or a portion of space, then we say 
that the field of the given quantity has been specified.'' 

The field is termed a scalar field if the quantity in 
question is a scalar quantity, that is, if it is fully de¬ 
scribed by its numerical value.- 

Examples of scalar fields are: a temperature field, an 
electrostatic field. 

Specifying a scalar field is accomplished by specifying 
the scalar function of a point M: 

u = / (M). 

If a Cartesian coordinate system xyz is introduced in 
space, we have 

u = f ( x , y, z). 

Geometrically, a scalar field is characterized by a 
level surface; this is a locus of points at which the scalar 
function of the field assumes the same value. The level 
surface of a given field is defined by the equation 

/ ( x , y, z) = C, where C = constant. 

In the case of a temperature field set up in a homo¬ 
geneous and isotropic medium by a point source of heat, 
the level surfaces are spheres centred at the source (this 
is a central-symmetric field). 
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In the case of an infinite uniformly heated wire, the 
level surfaces (isothermic surfaces) are circular cylinders 
whose axes coincide with that of the wire. 

Example 1. Construct the level surfaces of the scalar 
field 

u = x + 2y + 3z. 

Solution. The level surfaces are given by the equation 
x + 2y + 3z = C, where C = constant. 

This is a one-parameter family of parallel planes. 
Example 2. Find the level surfaces of the scalar field 
u = z 2 + y 2 — z 2 . 

Solution. The level surfaces are given by the equation 
I s + j/ 1 - ! s = C, where C = constant. 

For C = 0, we obtain a circular cone. For any C >0, we 
obtain a hyperboloid of revolution of one sheet with the 
axis coincident with the z-axis. For C < 0, we obtain 
a hyperboloid of revolution of two sheets. 

Example 3. Find the level surfaces of the scalar field 

u = arcsin T? . , . . 

K **+tf 2 

Solution. The domain of definition of the given scalar 
field is found from the inequality 

I , 1 Hy that is, 0< . <1, 

I Y*'+v' I **+«* 

whence 0 z 2 ^ z 2 + y 2 - This double inequality shows 
that the field is defined outside a circular cone z 2 = 
= z 2 + y 2 and on it, with the exception of its vertex 
0 (0, 0, 0). 

The level surfaces are found from the equation 

arcsin — * = C, where — 

/ x * +! ,* 2 2 

That is, z/j/z 2 + y 2 = sin C or z 2 = (z 2 + y 2 ) sin 2 C. 
This is a family of circular cones located outside the cone 
z 2 = z 2 + y 2 with a common axis of symmetry, Oz, and 
a common vertex, O (0, 0, 0), at which the given field 
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is not defined; note that the cone itself, z 2 = x 2 + y 2 , 
is also included in the family. 

Example 4. Find the level surfaces of the scalar field 
u = e< a > r >, 


where a is a constant vector and r is the radius vector 
of a point. 

Solution. Here, 

r = {x, y, z} = xi + yj + zk 

and let 

a = {a lt a 3 , a 3 } = + aj + a 3 k 

Then the scalar product 


(a, r) = a x x + a 2 y + o 3 z. 
The equation of the level surfaces is 
*<••'> = £:, C> 0, 


whence 


(a, r) = In C 
a,x + a t y + a 3 z = In C. 


This is a family of parallel planes. 

Find the level surfaces of the following scalar fields: 

45. u = 4 + 4 + ii. 

46. u = x 2 + y 2 — z. 


47. u = 


**+y' 


48. u = 2y 2 + 9z 2 . 

49. u — 3 x * 2V ~ r . 

50. u = I*’ ^ (a, b are constant vectors). 

51. u = In | r !. 

52. u = e< a - h - r > (a, b are constant vectors). 

A scalar field is said to be plane if there is a plane 
such that in all planes parallel to the given plane the 
scalar field is the same. 
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If we take this plane as the xy-plane, then the scalar 
field is given by the scalar function 

“=/(*. y), 

which is to say it is not dependent on z. 

An example of a plane scalar field is the temperature 
field of an infinite, uniformly heated wire. 



Fig. H 

Geometrically, plane scalar fields are characterized 
by level lines; these are loci in which the scalar function 
has one and the same value. 

Example 5. Find the level lines of the scalar field 

u = x t -y\ 

Solution. The level lines are given by the equations 
x 2 — y 2 = C, where C = constant. 

When C = 0, we obtain a pair of straight lines: 
y = x, y = —x. 

For C =£ 0, we obtain a family of hyperbolas (Fig. 11). 
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Find the level lines of the following plane fields: 

53. u = 2x — y. 

54. u = In j/ £ . 

55. u = -^-. 

56. u= e x '~ v '. 

57. Find the level lines of the scalar field u given im¬ 
plicitly by the equation 

u-t-ilnu + y = 0. 


Sec. 8. Directional derivative 

Suppose we have a scalar field defined by a scalar 
function 

u = / (Af). 

In the field, take a point M 0 and choose a direction 
indicated by the vector 1. Then in the field take another 
point Af so that the vector MqM is parallel to 1. Denote 
by A u the difference 

Au = / (Af) — / (M 0 ) 

and by A l the length of the vector M(>M. The ratio Au/Al 
defines the average rate of change of the scalar field per- 
unit of length' in the given direction. Allow the point M 
to move towards the point M 0 so that the vector MqM 
is always collinear with the vector 1. Then A l 0. 

Definition. If, as A1 0, there is a limit to the ratio 
Au/AI, then it is called the derivative of the function 
u = f (Af) at the given point M 0 in the direction of 1 
and is denoted by the symbol duldl so that, by definition, 
we have 

-g- = lim ~ = 1 jm /(M) ~ /(Afo) , M 0 M || 1. 

This definition of a directional derivative is invariant, 
that is, it is not connected with any choice of coordinate 
system. 
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Suppose a Cartesian coordinate system has been intro¬ 
duced in space and suppose the function / ( M) = 
= / ( x , y, z) is differentiable at the point M 0 (x 0 , y 0 , z 0 ). 
Then 

sr I = 4 r-| cos a + 4 ^-1 cosP + 4^1 cosy, ( 1 ) 
dl IM„ dx |m, dy Im, r dz Im, r v ' 

where cos a, cos P, cosy, the direction cosines of the 
vector 

1 = a t i + aj 4- a 3 k, 
are found from the formulas 


cos a = cos P = -jjy , cos y ~ *nV ’ 

|1| = Kaf + a * + a 3- 

The symbols -g- lh « l P"' 
tial derivatives are taken at the point M 0 . 

For a plane field u = f (x, y), the directional derivative 
1 at the point M 0 ( x 0 , y 0 ) is 


du 

dl 


sin a, 


( 2 ) 


where a is the angle formed by the vector 1 and the 
x-axis. 

Remark. The partial derivatives du/dx, duldy, du/dz 
themselves are derivatives of the function u in the direc¬ 
tion of the coordinate axes Ox, Oy, Oz respectively. 

Formula (1)—used to compute the directional derivative 
at a given point—holds true even when the point M 
tends to M 0 along a curve for which the vector 1 is the 
tangent line at the point M 0 . 

Example 1. Find the derivative of the scalar field 
u = xyz 

at the point M 0 (1, —1, 1) in the direction from M 0 to 
M i (2, 3, 1). 

Solution. We find the direction cosines of the vector 
M 0 M, = {1, 4, 0}, the length of which is |M 0 M,| = /l7, 
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cos a = , cos ft = , cos v = 0. 

The values of the partial derivatives of the function 
u = xyz at the point Af 0 (1, —1, 1) are 


dx |m„ ’ dy |M, ’ dz |m„ " 

Using formula (1), we get 


Y 17 ' Y 17 Y 17 * 


Af 0 increases in the given direction. 

Example 2. Compute the derivative of the scalar field 


u = arctan xy 

at the point M 0 (1, 1), which belongs to the parabola 
y = i 2 , in the direction of the curve (in the direction of 
increasing abscissas). 

Solution. The direction of 1 of the parabola y = x 2 at 
the point M 0 (1, 1) is the direction of the tangent to the 
parabola at that point (Fig. 12). 

Suppose the tangent 1 to the curve at M 0 forms with 
the x-axis an angle a. We then have 

y' = 2x, tan a = y' \ x ^ = 2, 


whence the direction cosines of the tangent line are 

cos a = 1 — = —^rzr , cos 6 = sin a 

V 1 + tan* a /5 r 

tan a _2_ 

“ Y 1 + tan* a ~ /5 ' 

The values of the partial derivatives of the given function 
u (x, y) at the point M 0 (1, 1) are 

du | _ 1 | _J_ du_ | 1 I _1_ 

dx |m, — l + *vl*. 2 ’ dy \M,—1 +x*V*\m, 2' 
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Substituting them into (2), we obtain 
du _ 1 1 1 2 3 

« 2 ' /5 + 2 ' /5 ~ 2 VI ‘ 

Example 3. Find the derivative of the scalar field 



u = xz 1 + 2 yz at the point M„ (1, 0, 2) along the circle 
x = i -fcosi, 
y = sin<—1, 
z — 2. 

Solution. The vector equation of the circle is of the 
torm 

r (t) = (1 + cos t) i + (sin t — 1 ) j + 2 k. 

We find the vector t tangent to it at any point M to be 
x = = — sin <-i + cost- j. 

The given point M 0 (1, 0, 2) is found in the xz-plane in 
the first octant and is associated with the value of the 
parameter t = n/2. At this point we have 

*|m.= — sin-j-i + cos y-j= -i. 
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From this we obtain that the direction cosines of the 
tangent to the circle are equal to cos a = — 1 , cos p = 0 , 
cos v = 0. The values of the partial derivatives of the 
given scalar field at the point M 0 ( 1 , 0 , 2 ) are: 

£L=< 2 “ +2 4.= 4 - 

Hence the desired derivative is 

3rL-wL- 4 - ( - i)+4 - 0+4 -°“- 4 - 

In the following problems, it is required to find, for 
given functions, the derivative at the point 
(x 0 , zo) » n the direction of the point (xj, j/„ Zj). 
u- V x 2 + y 2 + z*, Af 0 (1, 1, 1), Mi (3, 2Tl). 
u = x 2 y + xz*-2, M 9 (i, 1 , -1), Mi (2, -1,3). 
u = xe v 4- ye* — z 2 , M 0 ( 3,0,2), Af t (4, 1,3). 

61. u-i-i, M 0 (l,l), M ,(4, 5). 

62. Find the derivative of the scalar field 

u = In (x 2 + y 2 ) 

at the point M„ (1, 2) of the parabola y l = 4x in the 
direction of the curve. 

63. Find the derivative of the scalar field u = arctan ylx 
at the point M 0 (2, —2) of the circle x 2 + y 2 — 4r = 0 
along an arc of the circle. 

64. Find the derivative of the scalar field u = x 2 + y 2 
at the point M 0 (x 0 , y 0 ) of the circle x 2 + y 2 = R 2 in the 
direction of the circle. 

65. Find the derivative of the scalar field u = 2 xy + y 2 
at the point (|/2, 1) of the ellipse x*/4 + y 2 l2 = 1 in the 
direction of the outer normal to the ellipse at that point. 

66 . Find the derivative of the scalar field u = x 2 — y 2 
at the point (5, 4) of the hyperbola x 2 — y 2 = 9 in the 
direction of thb curve. 
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67. Find the derivative of the scalar field u = 
= In ( xy + yz + xz) at the point M 0 (0, i, 1 ) in the 
direction of the circle x = cos t , y = sin t, z = 1 . 

68 . Find the derivative of the scalar field u = x 2 + y 1 + 
+ z 2 at the point M 0 that corresponds to the value of 
the parameter t = jt /2 in the direction of the helical 
curve x = R cos t, y — R sin t, z = at. 


Sec. 9. The gradient of a scalar field 

Suppose we have a scalar field defined by a scalar 
function 

u = / (*, y , z) 

where the function / is assumed to be differentiable. 

Definition. The gradient of a scalar field u at a given 
point M is a vector denoted by the symbol grad u and 
defined by the equation 

* rad “=£*+^+1^ Cl) 

Using formula (1) of Sec. 8 for the directional derivative, 
we have 

^- = (grad u, 1 °), ( 2 ) 

where 1 ° is a unit vector in the direction of 1 , that is, 

1 ° = -ryj- = icosa + jcos p + kcosy. 


Properties of a gradient 

- 1. The gradient is in the direction of the normal to the 
level surface (or to the level line if the field is a plane 
field). 

' 2. The gradient is in the direction of increasing values 

of the function of the field. 

3. The modulus of the gradient is equal to the largest 
directional derivative at a given point of the field: 

maxifo grad “I - V + (%Y + (Ir)’ •* 

* The maximum is taken over all directions of l at the given 
point of the field. 
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These properties yield an invariant characteristic of 
the gradient. They state that the vector grad u indicates 
the direction and magnitude of maximum change of 
a scalar field at a given point. 

Example 1. Find the gradient of the scalar field 
u = x — 2y + 3z. 

Solution. By (1) we have 

grad u = li — 2j + 3k. 

The level surfaces of the given scalar field are the planes 
x — 2y + 3z = C\ the vector grad u = { 1, —2, 3} is the 
normal vector of the planes of this family. 

Example 2. Find the greatest steepn ess (rate) of rise 
of the surface u = x v at the point Af (2727 4). 

Solution. We have 

grad u *= yx«~H -f- x v In xj, grad u \ M = 4i + 4 In 2j, 
(IrLx = le rad u \ = ^y i + (in2) 2 . 

Example 3. Find the unit vector of the normal to the 
level surface of the scalar field 

u = x 2 + y 2 + z 2 . 

Solution. The level surfaces of the given scalar field 
are the spheres 

x* + j,* + z* = C (C > 0). 

The gradient is directed along the normal to the level 
surface so that grad u = 2x-i + 2g-j + 2z-k defines the 
vector of the normal to the level surface at the point 
M (x, y, z). For the unit vector of the normal, we obtain 
the expression 

n o _ gradu _ jl + yj + ^k _ _r_ 

I grad u | Yx* + y* + i* l r l ' 

Example 4. Find the gradient of the field u = (a, b, r), 
where a and b are constant vectors and r is the radius 
vector of the point. 

Solution. Let 

a = {a lt a 8 , n 3 ), b = {6^ b t , fc 3 }, r = {x, y, z}. 
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fc, b 2 b 3 

x y z 

By the rule for differentiating a determinant* we have 


a, a 2 a 3 

du _ 

a 2 a 3 

dx _ 

• dz 

a t a 2 a 3 

b\ b 2 b 3 

b t b 2 b 3 

b t b 2 b 3 

\ 0 0 


© 

o 

0 0 1 


a 2 a 3 a l a 3 a t a 2 

“= b, b, - 6, 6, 1+ 6, 


i j k 

a t a 2 a 3 

b t b 2 b 3 


= [a,b). 


* Given a determinant D (<) whose elements a,j are differenti¬ 
able functions of t: 



<*n M 

<*ll(0 • 

• <*in(0 

D(t)= 

<*«(<) 

<*»» (0 • 

• a tn (0 


<*ni (1) 

<*nl (0 • 

• ■ a nn (t) 


Then the derivative of the determinant, D' ( t ), is found from the 
formula 


<*«'i(0 <* is (0 •••<* i »(0 

<*ll(0 ««(<) ••• <*l*(0 


<*ni (0 <* 

ns U) ••• <*nn W 



<*il (0 <>ll (<) • 

•• <*in (0 

+ 

<*»(<) <*«(<) 

•• a'tnW 


<*ni (0 <*n» (<) • 

• • <*nn (0 



<*11 (0 

<*11 (0 • 

<*!»«) 

+ 

<*11 (<) 

<*11 (0 • 

<*ln(<) 


<*ni «) 

<*ns (0 • 

■ • “kn (0 
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Example 5. Find the gradient of the distance 

r — V (* — *o) 2 T (y - y 0 ) 2 + (2 - z «>) 2 . 

where P (a:, p, z) is the point of the field being studied 
and P 0 (x 0 , y 0 , z 0 ) is some fixed point. 

Solution. We have 


(x-xp) i + (y-y 0 ) j + (z-zo)k _ 0 
/(*— x 0 )* + (y— y,,)* + (z — zo)* 
which is the unit vector of the direction P 0 P. 



Example 6. Let us consider the scalar function 
u = Tj + r 2 , 

where r lt r 2 are the distances of some point P ( x, y) of the 
plane from two fixed points, F, and F 2 , of the plane. 

Solution. The level lines of this function are ellipses. 
We have (see example 5) 

grad (r, + r 2 ) = rj + r“. 

This shows that the gradient is equal to the diagonal of a 
rhombus constructed on the unit vectors of the radius 
vectors drawn to point P from the foci F x and F 2 (Fig. 13). 
Consequently, the normal to the ellipse at some point 
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bisects the angle between the radius vectors drawn to that 
point. 

Physical interpretation: a light ray coming from one 
focus enters the other focus. 

Example 7. Find the angle 6 between the gradients of 
the functions 

u = V x 2 + y 2 and v = x + y + 2 Y xy 


at the point Af 0 (1, 1). 

Solution. We find the gradients of the given functions 
at the point M 0 ( 1 , 1 ): 


grad u |m, = - 


i + i 
yl ' 


' Yz'+y' 

gr«d „ _ [ (1 + /I) 1 + (1 + /|) i] l- 2i+2j. 

The angle 0 between grad u and grad v at the point M 0 
is found from 

2 2 

a (gradu, gradw l Ml> /2 + /2 , 

005 |grad«| M> .|gradv| A1< - i. 2 /S “ 


From this we have 

0 = 0 . 


Example 8. Find the directional derivative of the 
radius vector r for the function u = sin r, where r = | r |. 

Solution. By (2), the directional derivative of the 
given function of the radius vector r is 

= (grad sin r, r°). (3) 

We find the gradient of the function: 

. . dfainr) . . d (sin r) . d (sin r) , 

grad 3in r = _L_±, + j + _i_^L k 

_ d (sin r) dr . , d (sin r) dr . . d (sin r) dr . 
dr dx 1 dr dy * ' dr oi 

=(-S- i +-|' i+ l- k ) cosr=r0cosr - ^ 
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Substituting (4) into (3), we get 

= (r° cos r, r°) = (r°, r°) cos r = cos r. 


Example 9. Find the derivative of the scalar field 
u = f (x, y , z) at the point M 0 (x 0 , y 0 , z 0 ) of the curve l 
specified by the system of equations 
/ (*. y, z) = a, \ 

L a = constant, 
q)(x,p, z)»0 J 


in the direction of the curve. 

Solution. The direction of the curve l is given by the 
direction of its tangent vector t, which, by definition, 
is a vector that is tangent to the surface / (x, y , z) = a. 
The surface f (x, y, z) = a is & level surface of the given 
scalar field u = / (x, y, z). Since 


~ = (grad u, 1 °) = (grad u, t°) 


and the vector grad u is perpendicular to the level surface 
/ (x, y, z) = a, it follows that grad u is perpendicular to 
the unit vector t°, and therefore 

ifL =( g radu ’ t °) | «.= o - 

Example 10. At the point M 0 (1, 1, 1), find the direc¬ 
tion of the greatest change in the scalar field u = xy + 
+ yz + xz and the magnitude of that change. 

Solution. The direction of the greatest change of the 
field is indicated by the vector grad u (M). We find it 
thus: 


grad u (M) = (y + z) i + (x + z) \ + (y + x) k 
and hence grad u (M 0 ) = 2 (i + j + k). This vector 
determines the direction of the greatest increase in the 
given field at the point M 0 (1, 1, 1). The magnitude of 
maximum change of the field at this point is 

max-^- = | grad u{M 0 ) | =21/3. 


69. Find the gradient of the scalar field u = In (x a + 
+y 2 + z 2 ) at the point M 0 (1, 1, —1). 

70. Find the gradient of the scalar field u = ze xi +y t +‘ 2 a t 
the point 0 ( 0 , 0 , 0 ). 


*-910 
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71. Find the angle 9 between the gradients of the func¬ 
tion u = arctan xly at the points M x (1, 1) and 
M t (—1, -1). 

72. Find the angle <p between the gradients of the func¬ 
tion u = (x + y) e* +v at the points (0, 0) and 

M * (1, 1). 

73. Find the angle 9 betw een the gradients of the func¬ 
tions u = Y x? -f y 2 + z a and v = In (i® + y 2 4- z a ) at 
the point M 0 (0, 0, 1). 

74. Find the points at which the gradient of the scalar 
field u = sin (x + y) is equal to i + j. 

75. Find the points at which the modulus of the gradient 
of the scalar field u = In V a? + y 2 + z a is equal to unity. 

76. Let u = u (x, y, z) and v — v {x, y, z) be functions 
differentiable at the point M (x, y, z). Show that 

(a) grad ( hi ) = X grad u, X — constant; 

(b) grad (u ± v) = grad u ± grad v ; 

(c) grad ( uv) = v grad u + u grad v; 

, d) er , d (±.) = » e r.d.-. g r.a„ [ „ #() 

77. Show that 

grad u(9) = -^-grad <p, 

where 9 = <p (x, y , z) is a differentiable function and 
u = u (9) has a derivative with respect to 9. 

Find the gradients of the following scalar fields if 

r = xi + y\ + zk, r- \ r | = V x z -(- y 2 + z 2 , 

and a and b are constant vectors. 

78. u = In r. 

79. u = (a, r). 

80. u = (a, r)-(b, r). 

81. u = | [a, r] | 2 . 

82. Show that 

(grad u (r), r) = u' (r)-r. 

83. Show that 

(grad u ( r), r] = 0. 

84. Let w = f (u, v), where u = u (x, y, z),v = v ( x, y, z). 
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Prove that 

grad w = grad u + grad v 

if /, u, v are differentiable functions. 

85. Suppose G is a convex region in space (that is, a 
region such that if two points M and N belong to G, 
then the whole line segment MN lies in G). Let there be 
given in G a scalar field u ( M) which at all points has a 
gradient that is continuous and bounded in G: 

| grad u ( M) | < A, M £ G, A = constant. 

Prove that for any points M and N of G we have the 
inequality 

| u (N) — u (M) | ^ A | MN |. 

86 . Find the derivative of the function u = xVa 2 + 
+ yVb 2 + z 2 /c 2 at an arbitrary point M (x, y , z) in the 
direction of the radius vector r of that point. 

87. Find the derivative of the function u = i/r, where 
r = | r | in the direction of the vector 1 = cos a-i + 
+ cos p-j + cos y-k. Under what condition is the 
derivative equal to zero? 

88 . Find the derivative of the function u = Mr, where 
r = | r |, in the direction of its gradient. 

89. Find the derivative of the function u = yze* at the 
point M 0 (0, 0, 1) in the direction of its gradient. 

90. Find the derivative of the scalar field 

u = u{x, y, z) 

in the direction of the gradient of the scalar field 
v = v (x, y, z). 

Under what condition is it equal to zero? 

91. For the following scalar fields, find the direction and 
magnitude of greatest change at the given points Af 0 : 

(a) u ( M) = i?y + y 2 z + z 2 x; M 0 (1, 0, 0). 

(b) u (M) = xyz; M 0 (2, 1, —1). 
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Sec. 10. Vector lines. Differential equations 
of vector lines 

Definition 1. We say that a vector field is given if a 
vector quantity a = a (Af) is specified at each point M 
of space or of a portion of space. 

If a Cartesian coordinate system is introduced in the 
space, then specifying the vector field a = a ( M) is equiv¬ 
alent to specifying three scalar functions of the point 
P (Af), Q (Af), R ( M) so that 

a (A/) = P (x, y,z)i + Q (x, y, z) j + R (x, y, z) k. 

Definition 2. A vector line of a vector field a is a 
curve at each point Af of which the vector a is directed 
along the tangent to the curve. 

Let a vector field be defined by the vector 
a = PI + Qi + Rk, 

where 

p = P (*, if, z), Q = Q (*, y. z), R = R ( X , y, z) 

are continuous functions of x, y, z that have bounded par¬ 
tial derivatives of the first order. 

Then the differential equations of the vector lines are 
of the form 

dx _ dy _ dz ... 

p ~ Q — R • 

Integrating this system of two differential equations 
(1) yields a system of two finite equations: 

Vi (*. y, z) = C lt <p, (*, y , z) = C % 
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which, taken together, define a two-parameter family of 
vector lines: 


y, z)=C„ 1 
<p 2 (x, y, z) = C 2 . ) 


( 2 ) 


If the conditions of the theorem of existence and 
uniqueness of a solution are fulfilled in a certain region 
G for system (1), then a unique vector line 


<Pi (*, y. z) = <p t (x 0 , y 0 . 2o). 1 
<p 2 (x, y, z) = <p 2 (x 0 , y 0 . Zo) I 


passes through every point M 0 (x 0 , y 0 , z 0 ) £ G. 

Example 1. Find the vector lines of the vector field 
a = [c, r], 


where c is a constant vector. 

Solution. We have 

c = Cji + Cjj + c 3 k, r = xi + yj + zk 


so that 
a = [c, r] 


i j k 

Cl c 2 c 3 
X y z 


= (c 2 z — Cay) i 4 - (c 3 x—c,z) j + (c,y - c 2 x) k. 


The differential equations of the vector lines are 

dx _ dy _ dz 

*—c»y ~ e a x-c,z ~ c,y — c t x * ' ’ 

Multiply the numerator and denominator of the first 
fraction by x, the second by y, the third by z and add 
termwise. Using a property of proportions, we have 

dx _ dy _ dz _ x dx + y dy + z dz 

c a z—e 3 y ~ c a x—c x z ~ c t y—c t x ~ 0 ’ 

whence 

xdx + ydy + zdz = 0 
and this means that 

x 2 -f y 2 + z* = j4 lt A 2 = constant >0. 

Now, multiplying the numerator and denominator of the 
first fraction of (3) by Cj, the second by e 2 , the third by c. 
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and adding termwise, we obtain 

dx _ dy dz _ Ci dx-\-c t dy-\-c 3 dz 

c t z—c t y ~ CfX—Cii e r y— c t x ~ 0 ’ 

whence 

Cj dx + c t dy + c 3 dz = 0 
and, consequently, 

c x x + c t y + c 3 z = A 2 , A t = constant. 

The desired equations of the vector lines are: 
z 2 4-y 2 + z 2 = A,, 1 
c 1 x + c 2 y + c 3 2=j4 2 . J 

These equations show that the vector lines are obtained 
via the intersection of spheres (having a common centre at 
the origin of coordinates) 
with planes perpen¬ 
dicular to the vector c = 
= c 2 i + cj + c 3 k. From 
this it follows that the 
vector lines are circles 
whose centres lie on a 
straight line passing 
through the coordinate 
origin in the direction of 
the vector c. The planes 
of the circles are per¬ 
pendicular to the indicat¬ 
ed straight line (Fig. 
14). 

Fl 8 - 14 Example 2. Find the 

vector line of the field 

a = -yi + x\ + 6 k, 

which line passes through the point ( 1 , 0 , 0 ). 

Solution. The differential equations of the vector 
lines are 

dx _ dy_ _ dz 
— y ~~x~~' 

whence we find 



x 2 + y 2 = C u C, >0 
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or, introducing the parameter t, 

x^VCiCOst, y=l/C,sini. 


In this case, the equation 

dy _ dz 


takes the form 

/Cl cos t dt 


—=-= —- or dz = bdt, 

V Ci cos t b 

whence we obtain 

z = bt + C t . 

Thus, the parametric equations of the vector lines are 
x=V(T l cost, 

y==VC x sint , (4) 

z = bt + C 2 . 

If we require that the vector line pass through the point 
( 1 , 0 , 0 ), we will have 

1 = 1^1 cost, 0 sin t, 0 = bt + C 2 . 


The first two equations of this system are satisfied for 
t = 2 kn, k = 0, ±1, . . . and for = 1. Taking k = 0, 

we get t — 0 and the last equation of the system yields 

C t = 0. The desired vector line passing through the 
point (1, 0 , 0 ) is 

x = cos t, 1 
y = sin t, > 

z = bt. J 

This is a helical curve. 

Find the vector lines of the following vector fields: 

92. r = si + yj + zk. 

93. a = a x i + aj + a 3 k, where a lt a t , a 3 are constants. 

94. a = (z — y) i + (x — z) j + (y — x) k. 

95. Find the vector line of the field 

a = x 2 i — i/ 3 j + z 2 k. 


which line passes through the point (V 2 , —Vg, 1 ). 



56 


Vector Fields 


[Ch. Ill 


A vector field is said to be plane if all the vectors of a 
are located in parallel planes and the field is the same 
in each of the planes. 

If a Cartesian coordinate system xOy is introduced in 
one of the planes, the vectors of the field will not contain 
any components along the z-axis and the coordinates of 
a vector will be independent of z, that is, 

a = P (x, y) i + Q (x, y) j. 

The differential equations of the vector lines of a 
plane field are of the form 

dx _ dy _ dz 

P(*. V) - <?(*. y)~ 0 
or 

dy Q(z, v) ^ 
dx ~ P (x, y) • 
z = constant. J 


From this it is evident that the vector lines of a plane 
field are plane curves lying in planes parallel to the 
xy-plane. 

Example 3. Find the vector lines of a magnetic field 
of an infinite current conductor. 

Solution. We will assume the conductor is in the direc¬ 
tion of the z-axis and that the current / flows in that 
direction. The intensity vector H of the magnetic field 
set up by the current is 

H=-£-[I,rJ, (5) 

where I = 7-k is the current vector, r is the radius vector 
of the point M (x, y, z), and p is the distance from the 
axis of the wire to the point Af. Expanding the vector 
product (5), we obtain 

The differential equations of the vector lines are 

dx dy dz 

~—y *" — "6”’ 
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x* + y*=R\ 1 

z = C. J 


That is, the vector lines are circles with centres on the 
z-axis (Fig. 15). 

Find the vector lines of the following plane vector 
fields: 

96. a = xl + 2yj. 

97. a = xi + zk. 

98. a = xi — t/j. 

99. a = 2zj + 4pk. 

100 . a = *•! + y 2 j. 

101 . a = zj — yk. 

The differential equations of 

the vector lines 

dx _ dy_ _ dz 

P ~ Q ~ R 


M(x,y,z) 


may be written as 


or, in vector form, as 

dr 


= a (M). 



Fig. 15 


This form of the equations of vector lines turns out to be 
convenient in the solution of a number of problems. 

Example 4. Find the vector lines of the field a = [c, r 1, 
where c is a constant vector. 

Solution. Applying ( 6 ), we get 


£ = (c,rj. 


(7) 


Forming the scalar product of both sides of (7) by c and 
using the properties of a mixed product, we find 


(*’W) = 


-(c, r) = 0 . 


( 8 ) 
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Similarly, forming the scalar product of both sides of (7) 
by r, we obtain 

r >-°- (9 > 
From equation (8) it follows that 
(c, r) = constant 

and from equation (9) it follows that 
(r, r) = constant. 

The vector lines are lines of intersection of the planes 
(c, r) = constant with the spheres r 2 = constant. 

Find the vector lines of the following vector fields: 

102. a = / (r)-r. 

103. a = (a 0 , r) b 0l where a 0 , b„ are constant vectors. 


Sec. 11. The flux of a vector field. 

Methods of calculating flux 

1. The flux of a vector field. Suppose we have a vector 
field 

a (M) = P (x , y, z) i + Q ( x , y, z) j + R (x, y , z) k, 
where the coordinates P ( x , y, z ), Q (x, y, z), R (x, y, z) 
of the vector a ( M ) are continuous (the field a ( M) is 
continuous] in some region G. Let S be a smooth or piece- 
wise smooth two-sided surface in which a definite side 
has been chosen (an oriented surface). 

Definition. The flux IT of a vector field a (M) through 
an oriented surface S is defined as the surface integral 
of the first kind, over the surface S, of the projection 
of the vector a ( M) by the normal n ( M) to that surface: 

n = | j pr„a dS = f j (a, n°) dS, 

s V 

where n° is the unit vector of the normal n to the chosen 
side of the surface S\ dS is the area element of the sur¬ 
face • S. 

In the case of a closed surface, we will always choose 
the outer normal n that is directed outwards from the 
region bounded by the surface S. 
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If a, p, y are the angles that the normal n forms with 
the coordinate axes Ox, Oy, Oz to the surface S, then 
the flux may be expressed in terms of a surface integral 
of the second kind: 

II = j j (a, n°)d£= j J [P{x, y , z)cosa 

+ Q ( x, y, z) cos p + R ( x, y, z) cos y| dS 
or 

n = j J (a, n°)dS= J j P(x, y, z)dy dz 

S B 

+ Q (x, y, z) dxdz+ R (x, y, z) dx dy, 

where 

cos a dS — dy dz, cos P dS = dx dz, cos y dS = dx dy. 

Basic properties of the flux of a vector field 

(a) The flux reverses sign when the orientation of the 
surface is changed (that is, when the orientation of the 
normal n to the surface S is changed): 

f j (a, n°) dS = — f j (a, n°) dS, 
s+ s- 

where S + is the side of,the surface S on which the normal 
n is chosen, and S~ is the side of S on which the normal 
—n is taken (see [7]). 

(b) Linearity: 

[ [ (Xa + nb, n°)d5 = X J j (a, n°)dS + n j j (b, n°)dS, 
's' s s 

where \ and p are constant numbers. 

(c) Additivity: if the surface S consists of several 
smooth parts S lt S 2 , . . S m , then the flux of the vector 
field a ( M) through S is equal to the sum of the fluxes 
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of the vector a (Af) through the surfaces S lt S t , . . ., S m : 

11=2 ( J(a.n«)dS. 

This property permits extending the notion of flux 
to piecewise smooth surfaces. 



Example 1. Find the flux of the vector a = i through 
an area perpendicular to the x-axis and having the shape 
of a rectangle with sides 1 and 2 (Fig. 16) in the positive 
direction of the x-axis. 

Solution. According to the definition of the flux of 
a vector through a surface S, we have 

n = j J (a, n°) dS. 

In our case, a = i, n° = i so that (a, n°) = (i, i) = 1. 
Taking into account that the area of the rectangle is 
equal to 2, we obtain 

n= J j ldS-2. 

s J 

Remark. If we had chosen the unit vector of the normal 
to the area S so that n° = —i, we would have got II = —2. 
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Example 2. Compute the flux of the vector field 
a = r, where r is the radius vector, through a right cir¬ 
cular cylinder of altitude h, base 
radius R and the z-axis. 

Solution. The surface S consists 
of a lateral surface ct ( , an upper 
base a 2 and a lower base o 3 of the 
cylinder. By the additivity proper¬ 
ty, the desired flux II is equal to 
n = Ilj + n 2 -f n 3 , where F^, Il 3 , a 
n 3 are the fluxes of the given field 
through On a 2 , a 3 respectively. 

On the lateral surface o t of the cyl¬ 
inder the outer normal n° is paral¬ 
lel to the zy-plane and therefore 

(a, n°) = (r, n°) = pr n »r = R 



Fig. 17). Hence 


Fig. 17 


n, = j j (a„ n°) <£S = /f j j dS = /?- 2nR/i = 2nRVi. 


On the upper base <j 3 the normal n° is parallel to the 
z-axis and therefore we can put n° = k (see Fig. 17). 
Then 

(a, n 8 ) = (r, k) = pr 0 , r = h 

and so 


11 2 = j j (a, n°)dS = h j j dS = h-nR 2 = nR' i h. 

o, o. 


On the lower base o 3 the vector a = r is perpendicular 
to the normal n° = —k. Therefore (a, n°) = (r, —k) = 0 
and 

n 3 = j j (a, n 8 )dS = 0. 

The desired flux is then 


n = ^(a, n°) dS = 3nR 2 h. 
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Example 3. Find the flux of the vector field 


a 


through a sphere of radius R with centre at the origin 
of coordinates. 

Solution. Since the normal n to the sphere is collinear 
with the radius vector r, we can take n° = r° = r/| r |. 
Therefore 


(*•"") '(TTF'77r)=T7F (r ' r, = T7F = T7i*' 

Rut on the sphere S we have | r | = R, and so (a, n°) = 


= \IR\ 

The desired flux II is 


n = ^) (a, = = 

s s 

since the area of the whole sphere S is equal to IffydS = 
= 4ni? 2 . 

104. Compute the flux of the vector a = 3j through an 

area having the shape of a triangle with vertices at the 
points M 1 (1, 2, 0), (0, 2, 0), M, (0, 2, 2) in the direc¬ 

tion of the coordinate origin. 

105. Find the flux of the vector 

a = ai + 0j + yk, 

where a, (J, y are constants, through an area perpendicular 
to the z-axis and having the shape of a circle of radius R, 
in the positive direction of the z-axis. 

106. Find the flux of the vector a = r through the outer 
side of a circular cone whose vertex lies at the origin of 
coordinates; the base radius is equal to R and the altitude 
is h (the axis of the cone is along the z-axis). 

107. Find the flux of the vector a = / (| r |) r through 
a sphere of radius R with centre at the coordinate origin. 

II. Methods of computing the flux of a vector. 

1°. The method of projection onto one of the coordinate 
planes. Let an open surface S be projected one-to-one 
onto the x^-plane into a region D xu . In this case, the 
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surface S may be given by the equation z = f (x, y) and 
since the area element dS of the surface is 


dS 


dxdy 

| cos YI ’ 


it follows that computing the flux II through the chosen 
side of the surface S reduces to computing a double integral 
via the formula 

n -y<*•*■» (, » 

Here the unit vector n° of the normal to the chosen side 
of the surface S is found from the formula 


gradU—/(*, »)) 

• | grad !*—/(*, 0)1 I ' 


/(*)■+(*)*+ 


( 2 ) 


and cos y is equal to the coefficient of the unit vector k 
in formula (2): 


■/(*)■+(*)*+*' 


( 3 ) 


If (he angle y between the z-axis and the normal n° is 
acute, then in (2) and (3) the plus sign is taken, and if 
the angle y is obtuse, the minus sign is taken. The symbol 


(«. n°) I 

| COSY | |z=/<x >V ) 


signifies that / (z, y) must be substituted for z in the 
integrand. 

If it appears to be convenient to project the surface S 
onto the yz- and xz-planes, use is made of the following 
formulas to compute the flux II: 


n = 


ii 


(a, °°) 

| cos a | 


|x=«p(lf. z) 


dydz 


( 4 ) 
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or 


H — f [ , (a ’ ° 0) . I dxdz. (5) 

£ J I COS P I ly-W*. *) ' ' 

Formula (4) is used when the surface S is projected 
one-to-one into the region D yz of the yz-plane, which means 
that it may be given by the equation x = 9 (y, z); cos a 
is found as the coefficient of the unit vector i in the for¬ 
mula 

n o_i_ grad [a —q>(y, 2)] 

I grad [*-q>(v, *)] I 

That is, 

cos a = ± — — 1 . 

/•+(-£)’+(-&)' 

The plus sign is taken if the angle a between the z-axis 
and the normal n° is acute, and the minus sign if the 
angle a is obtuse. 

Formula (5) is used in the case of a one-to-one projec¬ 
tion of the surface S onto the zz-plane; in this case, 
S may be specified by the equation y = tJj {x, z) and 
then 



n°=± 


grad ly —1 |> (t, a)] 

| gradfi/—i|>(x, z)\ | 




cos f) is equal to the coefficient of the unit vector j in 
this formula, that is, 




If the angle f) between the y-axis and the normal n° is 
acute, we take the plus sign, and if the angle p is obtuse, 
we take the minus sign. 

Remark. When the surface S is specified implicitly by 
the equation $ (x, y, z) = 0, the unit vector of the 
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normal 

n° = i cos a + j cos p + k cos y 
is found from the formula 
-o _ | grad 0> (a, y, z) 

I grad <I> (x, y, 2 ) | 




where the sign on the right is determined by the choice 
of the normal to the surface S. 

To compute the flux II of a vector field a through a 
surface S, it is necessary to project the surface one-to-one 
onto one of the xy-, xz-, yz-planes; this is possible if the 
equation <D (z, y, z) = 0 is uniquely solvable with respect 
to z (z = / (z, y)), y (y = i|) (z, z)) or z (z = q> (y, z)) 
respectively. Then take advantage of one of the formulas 

(1), (4), (5). 

Example 4. Find the flux of the vector field 
a = (z — 2z) I + (z + 3y + z) j + (5z + y) k 


through the upper side of the triangle ABC with vertices 
at the points A (1, 0, 0), B (0, 1, 0), C (0, 0, 1). 

Solution. The equation of the plane in which the triangle 
ABC lies is of the form x + y + z = 1, whence z = 
= 1 — x — y. The triangle ABC is projected one-to-one 
onto the zy-plane into the region D xy , which is the triangle 
OAB (Fig. 18). 

It is given that the normal n° to the plane in which the 
triangle ABC lies forms an acute angle y with the z-axis 
and so we take the plus sign in ( 2 ) and obtain 


grad (z + y+z—1) 

' | grad(*+y + z— 1 ) I 
We find the scalar product 


75 i+ vT j + 7I k ' 


(6) 


(a, n°) = (z- 2z) + (z + 3y + z) + (5z + y) y=- 


/ 3 


5-910 
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From formula ( 6 ) we find that cos y = \l\ 3 >0 and, 
hence, 

dS-^-VSdzdy. 

Using formula (1), we compute the desired flux: 
n= j j (a, n°)dS = | j (7x + 4j/ —z) \ z= i- x - v dxdy 

= j j (8x + 5y-i)dxdy= j dx j (8* + 5y-l)dy = |. 

Dxu 0 0 

Example 5. Find the flux of the vector a = y 2 j + zk 
through the portion of the surface z = x 2 + y 2 cut off by 



the plane z = 2. The outer normal is taken with respect 
to the region bounded by the paraboloid. 

Solution. The given surface (a paraboloid of revolution) 
is projected one-to-one onto the xy-plane into the circle 
D xu (Fig. 19). We find the unit vector of the normal n° 
to the surface S: 

n o _ .j. grad (z— x*—y*) _ —2xi—2yj + k 

^ | grad(z— x*— y*) | /4x»+4y»+l ' 
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It is given that the normal n° forms an obtuse angle y 
with the z-axis and therefore the minus sign is taken in 
front of the fraction. Thus, 

2*l+2yj-k 
/ 4z* + 4y* + l * 


whence 


C ° 3T -- V4..+V + . <0 



dS = TTSyT “ V i±> + 4 » ! +1 dx dy. 
We find the scalar product 


(a, n°) = 


2 y»- 


Y 4* , +4y*+l 
The desired flux is, by (1), equal to 


n = j j (a, n°) dS = j J (2 y*-z) \ z=x . +i) ,dxdy 
s hxu 

= j J (2y*-y*-x*)dxdy. 
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The domain of integration D xy is a circle of radius R = 
= V 2 with centre at the coordinate origin. Introducing 
polar coordinates x = p cos <p, y = p sin 9 , we have 

II = j j (2p 3 sin 3 <p— p 2 ) pdp dqp 

Dxy 

2n 1/2 t v~ 

= j d<p J ( 2 p 4 sin 3 9 —p 3 )dp= — 2n-^-| o 2 = — 2n. 

Example 6 . Find the flux of the vector field a = 
= i — j + xyzk through the circle S obtained by cutting 



Fig. 20 


the sphere x* + y 1 + z* < R* with a plane y = x. Take 
the side of the circle facing the positive portion of the 
x-axis. 

Solution. Since the plane y = x is perpendicular to the 
xy-plane, the circle S lying in that plane is projected onto 
the xp-plane into the line segment A X A % and so the one- 
to-oneness of the projection is disrupted. The circle S 
is projected one-to-one onto the other coordinate planes. 
Projecting the circle onto the xz-plane, we obtain a region 
D xz bounded by an ellipse (Fig. 20). The equation of the 
ellipse can be found by eliminating y from the system 
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of equations 


whence 


x z + y 2 + z 2 = / i 2 t ■> 

y = x, I 


2 x* + z 2 = R* 



2 


= 1. 


It is given that the normal to the circle S forms an 
obtuse angle f) with the y-axis (see Fig. 20) and so we 
take 


n = — grad(y —i) = i —j, 


’ /2 * VI 


From the latter equation we have cos p = —1/^2 < 0. 
The area element dS of the circle is equal to 


dS = 


= V2 dxdz. 


dxdz 
’ I cosp I 

We find the scalar product: (a, n°) = V2. 
The desired flux is, using formula (5), 

Jlit* 


n= 


■.(I 

Dxz 


£j‘ 


2dxdz = 2 dxdz = 2 


"tt 


V2 &n 


since the area Q of the region D xz bounded by an ellipse 
with semi-axes a = R/V% and b = R is equal to 


Q = 


L! 


dx dz = nab 


nR* 

/2 * 


Example 7. Compute the flux of the vector a = 
= xi + yi + zk through the outer side of the lateral 
surface of the circular cylinder x 1 + j/ a = R 2 bounded 
by the planes z = 0 and z = H {H > 0). 

Solution. The given cylinder is projected onto the 
xy-plane into a line, namely, into the circle (Fig. 21) 

-y 2 =fl 2 , 1 

z = 0. } 
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We will therefore project the cylinder onto the other 
coordinate planes, for instance, the yz-plane. Since the 
cylinder does not project one-to-one onto the yz-plane, we 
take advantage of the additivity property of the flux of 
the vector and represent the de¬ 
sired flux II as a sum of 
fluxes: n=nj-|-nj, where IIx is 
the flux of the field through the 
portion S r of the cylinder located 
in the region where y > 0, and 
n s is the flux of that field 
through the portion S 2 of the 
cylinder located in the region 
where y <0. On S x we have 

and so 

n, = j j RdS = R j j dS = RS, 
s, s, 

where S is the area of the por¬ 
tion Sj of the cylinder. Since 
S = nRH, it follows that n x = nR 2 H. 

On S j we again have 

(.,»a 

and so 

n 2 = j j RdS = RS = nR*H. 

The desired flux is II = 2 nR*H. 

Remark. The solution is made simpler if we introduce 
curvilinear coordinates x = R cos <p, y = R sin q>, z = z 
on the cylinder (see item 3° below). 

To find the flux of the vector field a = P (x, y, z) i + 
+ Q (x, y, z) j + R (x, y, z) k through the surface S 
specified by the equation z = / ( x, y) by the method of 
projection onto a coordinate plane, it is not necessary to 


2 * 



X 


Fig. 21 
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find the unit vector of the normal n°, and we can take 
the vector 

n = ±grad[z-/(x, y)] = ± + 

Formula (1) for finding the flux II becomes 

11 = J J (*’ n °) dS = 5 J (a ’ n) U. „) dx dy% (7) 

In similar fashion we obtain formulas for computing 
the fluxes through surfaces given by the equation x = 
= <p (y, z) or y = t|> (x, z). 

Formula (7) is written thus in coordinate form: 

n=±\\{-P[x,y,nx.y) ] £-Q [x, y, f (x, y)l -g- 

Dxu 

+ R\x, y, /(x, y))}dxdy. 


Example 8. Compute the flux of the vector field 
a = xl + pj + Yx 2 -\-y 2 — 1 k 


through the outer side of the hyperboloid of one sheet 
z = /x 2 + y s — 1 bounded by the planes z = 0, z = YZ. 

Solution. The given surface is projected one-to-one onto 
the xy-plane into the region D xy bounded by the circles 


y ‘ = l ' } and ^ + = 

= 0, J z = 0. 


We find the outer normal n: 
n = ± grad (z— Vx z + y*— l) = ± | 


-xt-yj 


1 y x i +!/ »_-i +k )* 
Since n forms an obtuse angle y with the z-axis (Fig. 22), 
we take the minus sign and, hence, 

— *+«i. 

Yx*+y t — 1 

We find the scalar product 
**+]/* 


-k. 


(a, n)= 




V x 2 + y 2 — i = 


/*■+*■-1 
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Using formula (7), we obtain 

Passing to the polar coordinates x = p cos <p and y = 
= p sin <p, we have 

2n 2 

= 2nVp 2 — 1 |j =2Y3n. 

Example 9. Compute the flux of the vector field 
a = yi + zj + xk 


through a closed surface bounded by the cylinder i 2 + y 2 
= R 2 and the planes z = x, z = 0 (z > 0). 



Fig. 22 


Solution. The surface S is piecewise smooth and so we 
take advantage of the additivity property of a flux re¬ 
presenting the desired flux II as a sum of fluxes n x , n 2 , n 3 
through the smooth portions S x (semicircle x 2 + y 2 < R 2 , 
0 ^ x ^ R, z = 0), S t (part of the plane z = x), and 
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S 3 (part of the cylinder x 2 + y 2 = /? a ): II = rij + II 2 + 
+ fl s . Since S is closed, we take the outer normal to it 
(Fig. 23). 

(1) On S lt where z = 0, we have n° = — k and so 
(a, n°) = —x. 

This means that the flux 

II,-j j xdS= - j j xdxdy. 

Passing to the polar coordinates x = p cos <p and y = 



Fig. 23 


= p sin <p, we find 
n, = — j j p cos qpp dp dcp 

n/2 J? 

= — | costpdcp j p 2 dp= — 7? 3 . 

-n/2 0 

(2) On iS 2 , where z = x, we have 

n = ±grad (z — x) = ± (—i + k) 
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and since the normal n to S 3 forms an acute angle with 
the z-axis, we take the plus sign in the right-hand member. 
Thus, n = —1 + k and, hence, (a, n) = x — y. 

Projecting S 2 onto the zy-plane, we get the semicircle 
D xy : 0<ar<]//?2-y 2 . 

Then by (6) we have 

n 2 = J j (a, n)| ^ ^ dxdy , 

Dxy 

and again passing to polar coordinates, we find 

n/2 R 

n 2 = \ (cos <p—sin <p) dtp l p 2 dp = 4- R 3 . 

—n/2 0 

(3) On S 3 , where x* + y 2 = R 3 , that is, on the lateral 
surface of the cylinder, we have 

o _ , grad(i« + y«-J?«) _ al + yj _ , xl+yl 

n ~ == | grad + | ~ = Y^+7' ~ * ’ 

In this case, the normal n forms a right angle with the 
z-axis and therefore cos y = 0 and so the choice of sign 
in the right-hand member is arbitrary. Take the plus 
sign and then 

and so 

n 3 = i-J j (x + z)ydS. 

It is impossible to project the surface S 3 (right cylinder) 
onto the * (/-plane since it projects into a line, a semicircle 
/'(the one-to-one nature of the projection will be upset). 
The same occurs when projecting onto the xz-plane. We 
therefore project the surface S 3 onto the yz-plane, onto 
which it is projected one-to-one into the region D yz 
bounded by the line 

X 2 + y2 = R2, | 
z — x.- I 

Eliminating x from this system, we obtain the equation 
for the projection of this line onto the zy-plane: z 2 + y 2 = 



Sec. 11] 


The flux of a vector field 


75 


= R 2 is a circle. Since 

| cos a I = I cos (nfTi) I = I K - 0 I = | jf \ = -j (x>0), 
we will have 


= \\l2-dydz-2\\ 

Dyz Dyz 


y dy dz. 


Using the polar coordinates y = p cos (p and z = p sin 9, 
we get 

n R 

n ’- 2 SJ p cos cp p dp d<p = 2 j cos <p dtp j p 2 dp = 0. 

0 0 

Thus, 

11= -4n* s + 4*fl s + 0 = 0. 


108. Compute the flux of the vector field a = yi + zj + 
+ xk through the upper side of a triangle bounded by the 
planes 

x + y + z = a, x = 0, y = 0, z = 0. 

109. Compute the flux of the vector field a = xzi through 
the outer side of the paraboloid z = 1 — x 2 — y 2 bounded 
by the plane z = 0 (z ^ 0). 

110. Compute the flux of the vector field a = xi + zk 
thro ugh the lateral surface of the circular cylinder y = 
= v R 2 — X 2 bounded by the planes z = 0, z = h {h > 0)» 

111. Compute the flux of the vector field a = xi + yj + 
+ zk through the upper side of a circl e cut ou t of the 
plane z = h (h > 0) by the cone z = V x 2 + y 2 . 

112. Compute the flux of the vector field a = 3xi — yj — 
— zk through the outer side of the paraboloid x 2 + y 2 = 
= 9 — z located in the first octant. 

113. Compute the flux of the vector field 

a = (x 2 + y 2 ) i + (y 2 + z 2 ) j + (z 2 + x 2 ) k 
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through the portion of the plane z = 0 bounded by the 
circle x 2 + y 2 = 1 in the direction of the unit vector k. 

114. Compute the flux of the vector field a = yzi — xj — 

— yk through the total surface of the cone x 2 + y 2 = z 2 
bounded by the plane z = 1 (0 z 1 ). \ 

115. Compute the flux of the vector field a = 2xi + 
+ (1 — 2y) j + 2 zk through the closed surface bounded 
by the paraboloid x 2 + z 2 = 1 — 2y (y ^ 0 ) and the 
plane z = 0 (z 0 ). 

116. Compute the flux of the vector field a = x®i + 
+ y 2 ] + z*k through the total surface of a pyramid 
bounded by the planes x + y -f z = 1 , x = 0 , y = 0 , 
z = 0 . 

117. Compute the flux of the vector field a = xl + yj + 
+ zk through the sphere x 2 + y 2 + z 2 = fl 2 . 

2°. The method of projection onto three coordinate planes. 
Suppose a surface S is projected one-to-one onto all three 
coordinate planes. Denote by D sv , D xt , D yz the projec¬ 
tions of S onto the xy-, xz-, yz-planes respectively. 

In that case the equation F (x, y, z) = 0 of the surface 
S is uniquely solvable for each of the arguments x, y, z 
so that 


x = x (y, z), y = y (x, z), z = z (x, y). 

Then the flux of the vector 

a = P (x, y, z) 1 + Q (x, y, z) j + R (x, y, z) k 

through the surface S, the unit vector of the normal of 
which is 

n° = cosa-i + cos P-j + cosy-k, 
can be written thus: 
n = j j (a, n°) dS 

= j j [ P (x, y, z) cos a -4- Q (x, y, z) cos p 


+ R (x, y, z) cos y] dS. ( 8 ) 
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We know that 

dS cos a = ±dy dz , 

dS cos p = ±dx dz, (9) 

dS cos y = ±dx dy, 


the sign in each of the formulas of (9) being chosen to 
coincide with the sign of cos a, cos P, and cos y on the 
surface S. Substituting (9) into (8), we get 


n= ± j z), y, z]dydz ±j y (x, z), z\dxdz 

Dyl Dxt 

± j j R [x, y , z (x, y )) dx dy. (10) 

Dxu 

Example 10. Find the flux of the vector 
a = xyi + yzj + xzk 

through the portion of the external side of the sphere 
i 2 + y 2 + z 2 = 1 located in the first octant. 

Solution. We have 

„o _ grad (z»+y» + z«— 1) _ xl+yj+ak _ . , , 

n — | grad (x*+y t +2*—1) | “ ^ x * +y t+ z * -^ + i/J+ z k, 
whence, taking into accqunt that the surface S lies in 
the first octant, we obtain 

cos a = x ;> 0, cos P = y > 0, cos y = z ^ 0. 
We therefore take the plus sign in (10) in front of the 
integrals, and putting 

P = xy, Q = yz, R = xz, 

we obtain 

n= j j xydydz+ j J yzdxdz+ j j xzdxdy. (11) 

Dyl Dxz Dx\) 

From the equation of the sphere x 2 + p* + z 2 = 1 we get 


z = z(x, y) = Vl — x 2 — y 2 , y = y(x,z)= V 1 —x 2 —z 2 , 
x = x(y, Z) = /l-p2_ Z 2. 


Substituting these expressions for x, y, z respectively 
into the third, second, and first integrals on the right of 
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(11), we get 

II = ^ j xY \ — x 2 — y 2 dxdy + ^ J z i —x 2 —z 2 dxdz 

+ j j yV\-y 2 -z 2 dydz. (12) 

Dyz 

Let us compute the first integral on the right and pass to 
the polar coordinates x = p cos cp and y = p sin r|>, 
where 0 < <p < n/2, 0 < p < 1. This yields 

7 ( = j j x V \ — x 2 — y 2 dxdy= j j p 2 V 1 — P 2 cos (p dtp dp 

Djty Dxy 

n/2 1 1 

= j cos<pd<p j p 2 |/l— p 2 dp= j p 2 V^l — P 2 <*p. 

0 0 0 

Setting p = sin t, dp = cos t dt in the last integral, we 

have 

n/2 n/2 

/, = j sin 2 1 cos 2 1 dt = ^ j sin 2 2t dt — . 

The second and third integrals in (12) are computed in 
similar fashion and we obtain 

1 2 = j j Z 1 — X 2 — Z 2 dx dz = -jg- , 

Dxi 

/ 3 = j j yV i-y 2 —z 2 dydz = -^. 

The desired flux is 

n=/ 1 +/ 2 +/ 3 =^-. 

118. Use the method of projecting onto all three coordinate 
planes to compute the flux of the vector field through a 
surface S. 

(a) a = zi — x\ + yk; 

S is the upper side of a bounded portion of the plane 
3a: + 6y — 2z = 6 cut out by the coordinate planes. 

(b) a = (x + y + z) i + (x + y + 2 - 1) j - 2k; 
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S is the upper side of part of the plane x + y + z — 1 = 0 
lying in the first octant. 

(c) a = (x— 1/y — z 2 )i + ] + (Vy~x 2 — z)k; 


S is the outer side of the paraboloid of revolution y = 
= x 2 + z 2 bounded by the plane y = 4 and lying in the 
first octant. 

3°. The method of introducing curvilinear coordinates 
on a surface. In certain cases, when calculating the flux 
of a vector field through a given surface S it is possible 
to choose a simple coordinate system on the surface 
itself to compute the flux instead of projecting onto 
coordinate planes. 

Let us consider some special cases. 

Case ( 1). Suppose a surface S is part of the circular 
cylinder x 2 -f y 2 = R 2 bounded by the surfaces z = 
= fi (*> y) and z = / 2 (x, y), and we have / x (x, y) < 

< 1% (z, y)- 

Setting 

x — R cos cp, y = R sin <p, z = z, 


we have for the given surface 
0 < <p < 2n, fi {R cos <p, R sin q>) < z < 

^ / 8 (R cos <p, R sin <p), 
and for the element of area dS we obtain the following 
expression (Fig. 24): 

dS = R dq> dz. 

Then the flux of the vector field a through the outer 
side of the surface S is computed from the formula 


where 


2n /.(fi cos <p, fi sin q>) 

n = /?jd(p j (a, n°)dz, 

0 /,(H cos q>. R sin q>) 


Example 


o _ grad (j« + y 8 — R*) _ art + yl 

“ I grad (*» + »• —/?*) | “ R ' 

11. Find the flux of the vector 
r = xi + yj + zk 


(13) 
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through the outer side of the lateral surface of the cir¬ 
cular cylinder x 2 + y 2 = R 2 bounded by the planes 
z = 0 and z = H (H > 0). 

Solution. Here we have 
0 < <p < 2 ji; /, (R cos <p, R sin <p) = 0 

f 2 (R cos <p, R sin q>) = H. 



Introducing curvilinear coordinates on the cylinder, 
we get 

x = R cos qp, y = R sin cp, z = z. 

Then the desired flux of the vector r is 

2n U 

n = ^? j dip j (r, n°) dz. 
o o 

But since 

r = xi + yj + zk = R cos <p-i + R sin <p-j + zk 
and the normal n° on the cylinder is 


n xi+y\ cos sin <p-j 

1 ~ R ~ R 


= cos<p*i-t-sin<p-j, 


it follows that the scalar product on the cylinder will be 
(r, n°) = R cos 2 qp + R sin 2 <p = R. 
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Finally, we obtain 

in U 

n = /?2 j d( p j dz=2n/?2//. 

U 0 

Example 12. Compute the flux of the radius vector 
r = xi + yj + zk 

through the lateral surface of the circular cylinder 
x 2 + y 2 = 1 bounded from below by the plane x + y + 
+ z = 1 and from above by 
the plane x + y + z = 2 . 

Solution. Here (Fig. 25) 
we have 

It = 1,/, (x, y) =1 — x- 

— yJt (x,y) = 2 — x — y. 

Passing to coordinates on 
the cylinder 
x = cos <p, y = sin <p, 
z = z, 

we get 

fi (*. y) — 1 ~ cos cp — sin <p, 

/, (x, y) = 2 - cos q> - 

— sin <p. 

According to (13), the flux 
of the vector r is 

2n 2-cos . _ 

II = ( dip j (r, n°) dz. 

But since on the cylinder x 2 + y 2 = 1 we have 
n° = xi + yj = cos q>-i + sin <p-j, 
it follows that 

(r, n°) = x 2 + y 2 = cos 2 <p + sin 2 cp = 1 
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and, hence, 

2n 2-cos (p-sln <f 2it 

n = j dcp | dz—\^d(f — 2 k. 

0 1-cos (p-sln ip 0 

119. Find the flux of the vector 

a = yi + xj — e*« z k 

through the outer side of the lateral surface of the cylinder 
x 2 + y 2 = 4 bounded by the planes z = 0 and x + y + 
+ z = 4. 

120. Find the flux of the vector 

a = x\ — xyj + zk 

through the outer side of the cylindrical surface x 2 + z 2 = 
= R 2 bounded by the planes y = 1 and x + y = 4. 

121. Find the flux of the vector 

a = i 3 ! — y 3 } + xz?k 

through the outer side of the cylindrical surface x 2 4 - y 2 = 
= 9 bounded by the sphere x? + y 2 + z 2 = 25. 

122. Find the flux of the vector field 

a = x\ — y\ — xyz? k 

through the outer side of the lateral surface of the cylin¬ 
der x 2 -f y 2 = 1 bounded by the plane z = 0 and the 
hyperbolic paraboloid z = x 2 — y 2 . 

123. Find the flux of the vector field 

a — (xy — y 2 ) i + (2x — x 2 + xy) j + zk 

through the outer side of the lateral surface of the cylin¬ 
der x 1 + y 2 = 1 bounded by the elliptic cone z 2 = 
= x-/2 + y 2 . 

Case (2). Suppose the surface S is a part of the sphere 
x l + y 2 + z 2 = R 2 bounded by conical surfaces whose 
equations in spherical coordinates have the form 0 = 
= /j (cp), 0 = / 2 (<p) and by the half-planes (p = (p lt 
<P = 

For the points of the given sphere, set 
x = R cos (p sin 0, y •= R sin <p sin 0, z = R cos 0, 
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where tpi ^ <p ^ cp 2 and 0j ^ 0 ^ 0 2 . Then for the 
element of area dS we obtain (Fig. 26) 

dS = R 4 sin 0 dQ dy. 


In this case, the flux of the vector field a through the 
outer part S of the sphere is found from the formula 
<p, e. 

n = R2 j d<p j (a, n°) sin 0d0, (14) 

9 , 0 , 


where 


o _ grad (x i +y*+z t —R t ) _ xt + yj + dt 
" I grad (**+»*+**-«*) | fl 

Example 13. Find the flux of the vector 



a = (* — 2y + 1) i + (2* + y — 3z) j + (3y + z) k 

through a part of the surface of the sphere x 3 + y 1 + z 3 = 
= 1 located in the first octant into a region where 
3 ? + y* + z 2 > 1. 
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Solution. Here we have 

R = 1, cpj m 0, <p 2 = n/2, 0! = 0, 0, = n/2, 
n° = xi + yj + zk, (a, n») = x 2 + y 2 + z 2 + x. 

On the sphere x 2 + y 2 + z 2 = 1 we introduce coordinates 
(p and 0 so that 

x = cos cp sin 0 , y = sin (p sin 0 , z = cos 0 . 

We then have 

(a, n°) = 1 + cos <p sin 0 
and, using (14), we obtain 

n/2 n/2 

H = j d( P I (l + cos<psin 0 )sin 0 d 0 

n/2 n/2 n/2 n/2 

= j dtp j sin 0 d 0 + j cos (p d<p j sin 2 0 d 0 = ^ n. 

124. Find the flux of the vector field 

a = + zk 

through the outer side of that part of the sphere x 2 + y 2 + 
+ z 2 = 1 cut out by the conical surface z 2 = x 2 + j/ 2 , 
z > yv + y 2 - 

125. Find the flux of the vector field 

a = yzi + xzj + xyk 

through the outer side of that part of the sphere x 2 + y 2 + 
+ z 2 = R 2 located in the first octant. 

126. Find the flux of the radius vector 

r = xi + yj + zk 

through the outer side of that part of the sphere x 2 + y 2 + 
+ z 2 = 2 that is bounded by the planes z = 0 and z = y. 

127. Find the flux of the vector 

a = xzi + yzj + z 2 k 

through the outer part of the sphere x 2 + y 2 + z 2 = 9 
cut off by the plane z = 2 (z ^ 2 ). 
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Sec. 12. The flux of a vector 
through a closed surface. 

The Gauss-Ostrogradsky theorem 

Theorem. If in some region G of space the coordinates 
of a vector 

a = P (*, y, z) i + Q (x , y, z) j + R ( x , y, z) k 

are continuous and have continuous partial derivatives 
dPIdx, dQldy , dRIdz, then the flux of the vector a through 
any closed piecewise smooth surface 2 located in G is equal 
to the triple integral of dPIdx + dQldy + dRIdz over 
the region V bounded by the surface 2: 

n=§<.,nv°-nj(if+f-+ir)< to 0) 

z v 

(the Gauss-Ostrogradsky formula). 

The normal n to the surface 2 is taken to be the outer 
normal. 

Example 1. Compute the flux of the vector 
a = xH + y* j + z 2 k 
through the closed surface 

** + y * + z * = z = 0 (z > 0). 

Solution. By formula (1), 

n= J j j (2x + 2y + 2 z)dv. (2) 

The integral (2) is conveniently computed in the spherical 
coordinates r, 0, cp. We have 

x = r sin 0 cos <p, y = r sin 0 sin <p, z = r cos 0 
and the element of volume is 

dv = r 2 sin 0 dr dQ dt p 

so that 

n ' 2 H5 (r sin 0 cos (p + r sin 0 sin <p 


+ r cos 0) r 2 sin 0 dr dQ dtp 
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2n ji/2 H 

= 2 j rfq> j sin 0 (sin 0 cos qH- sin 0 sin <p + cos 0) dQ j r 3 dr 
oo o 

2n ji/2 

= cos0sin0dO = ^. 

0 0 

Example 2. Compute the flux of the vector 
a = 4ri — y\ + zk 

through the surface of a torus. 

Solution. Taking advantage of the Gauss-Ostrogradsky 
theorem we find that the desired flux is 

n = § ( .,■*)*-{JJ (£+£+£)*-«'. 

s V 

where V is the volume of the torus. To compute the 
volume V, let us take advantage of the Guldin theorem 



Fig. 27 


on the volume of a solid of revolution, by virtue of which 
the volume is equal to the product of the area of the 
rotating figure into the path covered by the centre of 
mass of the figure during the rotation. 

Let i?! and R 2 be the inner and outer radii of the torus 
(Fig. 27). The area S of the circle, which during rotation 
forms the torus, is equal to 
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The path length described by the centre of mass (the 
centre of that circle) is the circumference l of a circle 
of radius (R 2 + R t )/ 2, that is, 

l = 2 n ^±h = n(R t + R 2 ). 


Thus, the volume V of the torus is equal to 
V = n 2 n (R 2 + R t ) = ^- ( R 2 - Rtf ( R 2 + R ,). 

The desired flux is 

Example 3. Using the Gauss-Ostrogradsky theorem, 
compute the flux of the vector field 

through the outer side of that part of the surface z = 
= 1 — x 2 — y 2 located above the xy-plane. 

Solution. In order to be able to apply the Gauss- 
Ostrogradsky theorem, close the given surface from 
below with a portion of the xy-plane that is bounded by 
the circle 

* 2 + J/ 2 = l, 1 
2 = 0. J 

Let V be the volume of the resulting solid bounded by 
a closed piecewise smooth surface a consisting of a part 
a x of the paraboloid of revolution z = 1 — x 2 — y 2 and 
a part <j 2 of the plane z = 0 (Fig. 28). 

The flux of the given vector through the surface a is, 
by the Gauss-Ostrogradsky theorem, equal to 


"“§(»■ n °)*=J|J( 


dP_ dQ_ 
dx ' dy 


dv. 


We find the sum 

dP , dQ , dR 2xy , 2x 2*(l + y) __ A 

dx T dy T Oi ~ i+y * 1 + y* 1+y* — 
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Consequently, the flux is 

n = § (a, n°) do — 0. 

By virtue of the additivity of the flux we have 
n = { j (a, n°) do + J J (a, n°) da = 0. 



From this the desired flux is 

n, = j j (a, n°)do-= — j j (a, n°) do. 

The flux IT, of the vector a through the circle x 2 4- y 2 ^ 
^ 1, z = 0 is equal to 

n 2 = j j (a, n°) do. 

Since on the plane z = 0, we have 

a = -j~j- i-f 2xarctan y ■ j— k, n°= — k, 
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and hence (a, n°) = 1 , it follows that the flux II 2 through 
the circle <x 2 is equal to the area of the circle o 2 : 

n 2 = j j da = n. 


The desired flux IT! = —II 2 = — ji. 

By appropriately closing the given unclosed surfaces 
and making use of the Gauss-Ostrogradsky theorem, com¬ 
pute the fluxes of the vector fields through the indicated 
surfaces (the normal to the closed surface is the outer 
normal). 

128. a = (1 — 2x) i + yj + zk; S: x 2 + y 2 = z 2 
(0 < z < 4). 

129. a = z 2 i + xzj + yk; S: x 2 + y 2 = 4 — z (z > 0). 

130. a = (y 2 + z 2 ) i — y 2 j + 2yzk; S: x 2 + z 2 = y 2 
(0 < y < 1). 


Sec. 13. The divergence of a vector field. 

Solenoidal fields 

The notion of the flux of a vector through a closed 
surface leads to the concept of the divergence of a field. 
This concept yields a certain quantitative characteristic 
of a field at each point in the field. 

Let M be a point of the field under study. Surround it 
by a surface 2 of arbitrary shape, for instance, a sphere 
of sufficiently small radius. Let the region bounded by 
the surface 2 be (1^ and its volume V. We consider 
the ratio 


(a, n®) da 

y -• ( 1 ) 


Definition. If the ratio (1) has a finite limit when the 
region (T) contracts to the point M, then this limit is 
termed the divergence of the vector field (the divergence 
of the vector a) at the point M and is designated as 
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div a ( M ). We thus have 

H (a, nO) do 

div a (M) = lim ——„ -. (2) 

(V)-*M V 

Formula (2) yields an invariant definition of diver¬ 
gence. This definition means that the divergence of the 
field a at the point M constitutes the volume density of 
the flux of the vector a at that point. 

The points M of the vector field a (A/) at which div 
a > 0 are termed sources, while the points at which 
div a < 0 are termed sinks of the vector field. 

The divergence of a vector field is a scalar function of 
the points of the field. 

If the coordinates of the vector 

a (M) = P (x, y, z) i + Q (x, y, z) j + R (x, y, z) k 


have continuous partial derivatives dPIdx, dQ/dy, dRldz 
in the neighbourhood of the point M (x, y, z), then, using 
the invariant definition of divergence, we find from the 
Gauss-Ostrogradsky theorem that 


(3) 


All the quantities in (3) are considered at the same point 
M (x, y, z). 

Using (3) for divergence, we can write the Gauss- 
Ostrogradsky theorem (see Sec. 12) in vector form: 


(a, n°) da = j j j div a dv (4) 


Example 1. Using the invariant definition, compute 
the divergence of the vector a = xi at the point O ( 0 , 0 , 0 ) 
taking for the surfaces a surrounding the point 0 the 
spheres a e of radius e centred at that point. 
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Solution. By the definition of divergence, we have at 
the given point 

§) (a. n°) do 

diva(0)= lim —-, 

<« e )-o "e 

where v t is the volume of a ball bounded by a sphere 
o e , or 


div a (0) = lim - 


But since the volume of the ball is v t = 4ne®/3, it fol¬ 
lows that 


^ (a, n°) do 

div a (0) = lim -. 

Let us compute the flux ^>(a, n°) do of the given vector 

through the sphere <x e . The unit vector of the normal n° 
to the sphere a e is directed along a radius of the sphere, 
and so we can put 



where r° is the unit vector of the radius vector r = 
= xi + y\ + zk, or 

n o_ xi + yj + zk 
e 


The desired flux is 

(a, n°) do = ^ dl. do. 
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Passing to coordinates on the sphere o t , 
x = e cos <p sin 0, y = e sin <p sin 0, z — e cos 0, 
we get 


) (a, n°) da=jj c °s»q>sin*8sinOdcpM 
2 n n 

= e 3 j cos 2 <p dtp j sin 3 cp d<p = y i 


Consequently 


diva(0) = lim-2--1. 


Example 2. Compute div r. 

Solution. We have r = xi + yj + zk, so that P = z, 
Q = y, R = z and, hence, by formula (3), 


div r 


dx . dy_ . dz 
dx + dy + dz 


3. 


Example 3. Compute div (u- a), where u ( M) is a scalar 
function and a (M) = P (x, y, z) i + Q (x, y, z) j + 
+ R (x, y, z) k is a vector function. 

Solution. Using formula (3), we get 


div (ua) 

d(uP) d {uQ) d (uR) 
dx ' dy ' dz 


'»%r+ P £ + u %- + <?-d 


I dP . dQ . 

=u (lF + lF + 


u div a + (a, grad u). 


Thus 

div (ua) = u div a + («, grad u). (5) 

Example 4. Find the divergence of the vector 

a ---= q> (r) r° = r, 
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where r = | r | is the distance from the coordinate origin 
to the variable point M (x, y, z). 

Solution. Using formula (5), we obtain 

di va = ^divr+(r, grad -Sil). 

Furthermore, 

div r = 3, grad ^=(^)'grad r = r -il^=^r». 
And so 

diva*^.3 + ( r -2^=^r«,r) 

^3^- + rq>,(r) ~ q>(r) = 2-gfi + y' (r). 

131. For what function ij> (r) will we have div if (r) r = 
= 2t|) (r)? 

132. Find div (r*r). 

133. Find the divergence of the vector field 

a = lc, rl, 

where c is a constant vector. 

134. Find 

div (r (w, r]), 

where w is a constant vector. 

135. For what function i|> (z) will the divergence of the 
field 

a = xzi + yj + \|> (z) k 

be equal to z? 

136. Find the flux of the radius vector r through the 
surface of a sphere. 

137. The electrostatic field of a point charge q is 



Compute div E. 

138. Show that 


$ (r, n°) da = V, 
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where V is the volume bounded by a closed surface 2. 

139. Prove that if 2 is a closed piecewise smooth surface 
and c is a nonzero constant vector, then 

cos (u, c)do = 0 , 
z 

where n is a vector normal to the surface 2 . 

140. Prove the formula 

^ (<pa, n°) da = j j j (cp div a + (a, grad q>)) dv, 
z v 

where cp = (p (x, y, z) and 2 is a surface bounding the 
volume V. Establish the conditions under which this 
formula is applicable. 

141. Prove that if the function u (x, y, z) satisfies the 
Laplace equation 

d*u , d*u , d i u _ n 
dx* + dy * + dz * ~ U ’ 

then 

z 

where du/dn is the derivative with respect to the direction 
of the outer normal to the piecewise smooth closed sur¬ 
face 2 . 

142. Prove that if the function u (x, y, z) is a second- 
degree polynomial and 2 is a piecewise smooth closed 
surface, then the integral 



is proportional to the volume bounded by the surface 2 . 

Find the flux of a vector field through the indicated 
closed surfaces: (1) directly, (2) via the Gauss-Ostrogradsky 
theorem. 


143. a = xi-(- zk; S : 


z = x 2 + y 2 , 
z-= 4. 
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f Z 2 = X 2 + y 2 , 

!44. ._«+«-* ft { z _„ H>Q 

145. a = ii —2j; ft I I” 6 ""** -9 !' 

I z 2 = x 2 - \-y 2 , z^O. 

r x 2 + z 2 = y 2 , 

146. a = j/ 2 i—*|—jfk; 5: { _ t (0<f<1) . 

( z 2 = i 2 -(-u 2 , 

147. a = xi + 2pj —zk; * { , = , 2 + ^ 

148. a = 2xi—(z —1) k; 5: { Xi+ f = \ 

l z = 0, z = l. 

fx 2 +y 2 + z 2 = 4, 

149. —2*l-»j + 2k; ft { 3z = l2+ „2( z> ^+l!). 

150. a = yxi + 2yj — zk; 5: z ! i-y 2 +z 2 =4. 

151. a = (x+z)i + (t/ + x) j + (z + x)k; 


ft { 


x 2 +y 2 = R 2 
z = y, z^O. 


152. a = 3xi 


{ 9 — z = x 2 + y 2 , 
*-0. y = 0, 2 = 
octant). 


0 (first 


153. a=(y —x)i+(z—y)j-l-(x—z)k; 5: ^ x—y+z=l 

154. a = xi— 2yj — zk; 5: 


{ x+y+z=l, 
x-y+z= 1 , 
i = 0 , z= 0 . 


[ l-z = x 2 + P 2 , 

[ z = 0 . 


Solenoidal fields 

Definition. If at all points M of a certain region G 
the divergence of a vector field (specified in G ) is zero, 

div a ( M ) = 0 , 

then we say that the field is solenoidal in that region. 
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Thus, a solenoidal field is, by definition, without 
sources and sinks. 

From the Gauss-Ostrogradsky theorem it follows that 
in a solenoidal field the flux of a vector a = a (M) 



through any closed surface o 
lying in the field is zero: 

$(a, n°) da ~ 0 . 

In a solenoidal field G, the 
vector lines cannot begin or 
end. They can either be closed 
curves or have ends on the 
boundary of the field. 


Fig. 29 The equation 


div a ( M) = 0 


is encountered in hydrodynamics and is termed the con¬ 
tinuity equation of an incompressible fluid. 

In this case, the amount of fluid exiting through some 
closed surface o is always equal to the amount entering, 
and the total flux is zero. 

Which of the following vector fields is a solenoidal 
field? 

155. a =a x (z* — jr*) 1 + y (x* — z s ) j -f z (y 2 — x*) k. 

156. a = yH - (x* + y*) j + z (3 y* + 1) k. 

157. a - (1 + 2xy) i - y*xj + (z*y - 2zy + 1) k. 

158. Show that the field of the vector 


E = £r« (r=V x* + y* + z*) 

is solenoidal throughout the region that does not contain 
the coordinate origin 0 ( 0 , 0 , 0 ). 

159. Under what condition will the vector field a = 
= <p (r) r be solenoidal? 

Suppose we have a field a ( M ), which is not necessarily 
solenoidal. In the field, consider a closed oriented con¬ 
tour L. The surface 2 containing the curve L as its edge 
will be called a surface spanned by the contour L. Let us 
agree to orient the normal n to the surface 2 so that the 
chosen circuit about the contour L will be seen from the 
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end of the normal as being counterclockwise (Fig. 29). 
160. Show that in a solenoidal field the flux of a vector 
a (M) does not depend on the type of surface 2 spanned 
by the given contour L and depends solely on the contour 
itself. 


Sec. 14. A line integral in a vector field. 

The circulation of a vector field 

Suppose we have a continuous vector field a = a ( M) 
and a piecewise smooth curve L on which a positive 
direction has been chosen (in other words, L is an oriented 
curve). 

Definition 1. The line integral from the vector a = 
= a ( M) along the oriented curve L is termed a line 
integral of the first kind (an integral over the arc length 
of a curve) of a scalar product (a, t°): 

j (a, t°) ds, 

where x° = x° ( M) is the unit vector of the vector tangent 
to the curve L whose orientation coincides with that 
of L; ds is the differential of arc length s of L. 

If r = r ( M) is the radius vector of an arbitrary point M 
of the curve L, then the line integral in the field a (M) 
may be written thus: 

j (a,t°)ds= j (a, dr), ( 1 ) 

If a rectangular-'coordinate xyz-system is introduced 
in the vector field, then r = xi + yj + zk, 

a = P (x, y, z) i + Q (x, y,z)} + R (x, y, z) k 

and the line integral (1) can be expressed in terms of a 
line integral of the second kind: 

j (a, dr) = j P ( x , y, z) dx + Q (x, y, z)dy + R (x, y, z) dz. 

When a = a ( M) is a force field, the line integral (1) 
yields the magnitude of the work of that field over the 
curve L. 


7-910 
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Properties of a line integral 

(a) Linearity: 

j (Xa, + pa 2 , dr) = X j (a,, dr) + | (a 2 ,dr), 

where X and p are constants. 

(b) Additivity: 

j (a,dr)* j (a,dr) + j (a,dr), 

Li+L| L\ L2 

(c) The integral reverses sign with a change in the 
orientation of L : 

j (a,dr)«- j (a. dr), 

BA AB 

where A is the initial and B is the terminal point of the 
curve L. 

Calculating a line integral in a vector field 
Let the curve L be specified by the parametric equations 

x = <p (0. y = ^ (0. * = x (0 (*o < * < *i)- 

Here, at the initial point A of L the parameter t assumes 
the value t = t 0 and at the terminal point B of L it 
assumes the value t = t t (the direction on L corresponds 
to increasing values of the parameter t from t 0 to tj); 
the functions <p (t), ( t ), x (0 Lave continuous deriva¬ 
tives on the interval U 0 , fj). Then 

j (a, dr) = j (a,dr)-|{P|q.(0,*(/),X«)lq»'(0 

L AB to 

+ Q I<p (0. ♦ (0. x (01 (<) + #[? (0. ♦ (0. x (01 x' (0} dt. 

If the curve L is given by a system of equations 
y = i|) (i), z = x (*). a < x < b, then 

j (a, dr) = j {P|x,^(«),x(x)] + 9[x,t|)(x).X(«)]V( x ) 


+ B [x, i|> (x), x(*)lx' (x))dx. 
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Similar formulas may also be written for cases where 
the line is specified by one of the following systems 
of equations: 

x = <p ( y ), z = x (y) (y 0 < y < yJ 

or 

x = 9 (z), y = (z) (z 0 < z < Zj). 

Example 1. Find the line integral of the vector a = 
= rl | r |, where r is a radius vector on the line segment 
from point A (r A ) to point B (r fl ). 

Solution. The desired line integral is 

<«•*)-1W- (,) 

AB AB 

From 

d (f> r) = (*, r) + (r, dr) = 2 (r, dr) 

we find 

(r, dr) = 4 d (r, r) = \ d (l r |*) = 2 l r I d |r| = |r| d |r |, 

whence 

^--nn. <2) 

Substituting (2) into the integral (l) t we get 

r H 

J (■•*)“ J d|r| = [ d|r| = |r*| —|r A |. 

AB AB i A 

Note that 

I dr | |r |. 

Find the line integral over the line segment bounded 
by the points A (r,) and B (r t ) for the following vector 
fields: 

161. a = r. 


163. a= , r° is the unit vector. 

164. Compute the line integral over the straight line 
passing through the points 0 (0, 0. 0) and M l (1, 1, 1) 
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from point 0 to point if a = lb, r] , where b is a con¬ 
stant vector. 

165. Prove that 

j (grad u,dr) = u(B) — u(A). 

Example 2. Find the line integral of the vector 
a = zi + + yk 

over an arc L of the heli¬ 
cal curve 

x=Rc.ost, y — Rs\nt, 

t 

Z ~ 2n 

from point A, the point of in¬ 
tersection of the curve with 
the plane z = 0, to point B, 
the point of intersection with 
the plane z = 1 (Fig. 30). 

Solution. Here, the line in¬ 
tegral is of the form 

j (a, dr) 

= j zdx + xdy + ydz. 

L 

The helical curve is located 
Fig. 30 on a circular cylinder z® + 

+ y 2 = R 2 . At point A we 
have t 0 = 0, at point B we have f, = 2 ji. Since 

dx= —Rsint dt, dy = Rcost dt, dz = -^- f 

it follows that the integral is equal to 
2n 

[ (a, dr) = j ( — sin £ + R 2 cos 2 < + sin <) d< 

L 0 

2n 2 n 

= R z j cos 2 tdt — j t sin / dt=nR 2 -\- R 
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because 

2n 

j cos 2 *d* = Ji; 

2n 

j t sin tdt= — 2ji. 
o 

Example 3. Find the line integral of the vector (see 
example 2) 

a = zi + xj + yk 

over the straight line AB (see Fig. 30) in the direction 
from point A to point B. 

Solution. Since the straight line AB (the generatrix 
of the cylinder x 1 + y 2 = R 2 ) is located in the arz-plane 
and passes through the point A ( R , 0, 0), it follows that 
y = 0, x = R, dx — 0, and for the radius vector r of 
the points of AB we will have r = Ri + zk, dr = dz k. 
Therefore the scalar product 

(a, dr) = z dx + x dy + y dz 

on AB will be zero. 

Hence, the desired line integral 

J (a, dr) = j (a, dr) 

L AB 

on AB will be zero. 

From examples 2 and 3 it follows that in the general 
case a line integral in a vector field depends not only on 
the initial and terminal points of the path of integration 
but also on the shape of the path. 

Example 4. Compute the work done by the force field 
F = yi + xi + (x + y + z)k 
along the segment AB of a straight line passing through 
the points Af x (2, 3, 4) and Af 2 (3, 4, 5). 

Solution. The work done by the force field will be 
equal to the line integral along MyM^: 

A= j (F, dr) = j ydx + xdy + (x4-y+ z)dz, 

MlMl M,M, 
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Let us find the canonical equations of the straight line 
M X M Z . We have 

i-2 _ y-3 _ i -4 
1 1 — 1 * 

whence 


y=*x + \A 
2 = 1 + 2, J 
dy = dx, dz = dx. 

Here, x varies from 2 to 3 (since the abscissa of M x is 
equal to 2 and the abscissa of A/ s is 3). The desired work is 

3 

A= | (x-hl+x + x + x + l+x + 2)dx 


= j(5x + 4)dx = -f-. 

166- In the plane vector field 


a = 

Vx'+y' 

compute the line integral along the semicircle 

x — R cos /, y = R sin t (0 ^ t ^ n). 
v*f. In the plane vector field 

a = (I s + y*) i + (** - y) j 

compute the line integral over the curve y = \ x \ from 
the point (—1, 1) to the point (2, 2). 

168. In the plane vector field 

a = (x 2 — 2xy) 1 + (j/ 2 - 2 xy) j 


compute the line integral: 

(a) along the parabola y = x 2 from the point (—1, 1) 
to the point (1, 1); 

(b) along a segment of the straight line joining the points 
(-1, 1) and (1, 1). 

169. Compute the work of the force field F = 2xyi + x*j 
along an arc of the circle x 2 + y 2 = 1 from the point 
(1, 0) to the point (0, 1) counterclockwise. 
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170. Compute the work of the force field 

F = (x 2 + 2xp) i + (x 2 + p 2 ) j 
along the parabola y = x 2 from the point (0, 0) to the 
point (1, 1). 

171. In the vector field 

a ___ ji+yj + ^k _ 

/**+y*+z’-x-y+2z 

compute the line integral along a segment of the straight 
line from the point (1, 1, 1) to the point (4, 4, 4). 

172. In the vector field 

^ a = (p 2 — z 2 ) i + 2pzj — x 2 k 

compute the line integral over the line 

x = t ] 

y=t\ ^ (0<<<i) 

Z = < 2 J 

in the direction of increasing values of the parameter t. 

173. In the vector field 

a = pi + zi + xk 

compute the line integral along a turn of the helical curve 

x — acos t, "j 

y = a sin t, l (0<*<2n) 

z = bt J 

in the direction of increasing values of the parameter t. 

174. In the vector field 

a = x 2 ! + p*J + z*k 

compute the line integral in the direction from the point 
(0, 0, 0) to the point (1, 1, 1) along the line segment 
between the two points. 

Computing the circulation of a vector field 

Definition 2. The circulation C of a vector field a = 
= a ( M) is a line integral taken around a closed oriented 
curve L Thus, by definition, 

C = ^ (a, dr), 
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where the symbol <y denotes the integral around the 

L 

closed curve L. 

If the vector field a = a ( M ) is gfven* in coordinate 
form, 

a = P (x, y, z) i + Q ( x , y, z) j + R ( x , y, z) k, 
then the circulation of the vector field is 
C = ^> P dx + Qdy+ Rdz. 

L 

For the positive direction of traverse around the closed 
curve L we take the direction in which the region bounded 
by the curve is on the left. 

Example 5. Compute the circulation of the vector 
field 

a = —y*l + X 3 ! 

around the ellipse L: a?la* + y a /6 a = 1. 

Solution. By the definition of circulation we have 

C = ^ (a,dr) = ^ — y 3 dx+x 3 dy. (3) 

L L 

The parametric equations of this ellipse are of the form 
x = a cost, ) 

, (0<*<2n), (4) 

p = 6sinM 

whence 

dx = —a sin t dt, dy = b cos t dt. (5) 

Substituting (4) and (5) into (3), we get 

2n 

C - ab j (b 2 sin 4 / + a 2 cos 4 t)dt = -j-nab (a 2 +<*) 


2k 

= (1—C032/) 1 * 

= | J (l-2cos2< + i±fiii-)d( 
0 


since 

2n 
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-t] (t~^ cos ^'*‘T c03 **) ^“T j T^ ( = T 

Similarly, we find that 

2n 

j cos 4 tdt = - 1- ji. 

0 

Example 6. Calculate the circulation of the vector 
field 

a = ye*> i + xe*"] + xyzk 

around the curve L obtained by cutting the cone x 2 + y 2 = 
= (z — i) 2 with the coordinate planes (Fig. 31). 



Solution. The curve L consists of two parts BC and CA 
located in the coordinate yz- and xz-planes respectively, 
x 2 + y 2 = 1 1 

and the arc AB of the circle g __ of • Therefore 
the circulation of the given vector field will be 

C = §( a, dr) = j (a, dr) + j (a, dr) + j (a, dr). 

L BC CA Tb 

(1) On the line segment BC we have 
x = 0, dx = 0; z = 1 — y, dz = —dy\ 1 > y > 0. 
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Consequently, 

j (a, dr) = j ydx^O. 

BC BC 

(2) On the line segment CA we have 

y = 0, dy = 0; z = 1 — x, dz = —dx; 0 < x < 1 
and so 

j (a, dr) = f xdy = 0. 

CA CA 

^ x 2 + y 2 = i,\ 

(3) On the arc AB of the circle g ^ | we have 

2 = 0, dz = 0, which means that 


j (a, dr) = j e xv (y dx + x dy) = j e xu d(xy) 

AB AB AB 


= j d(e x «) 

AB 


e xv 


|B(0. I) 
|A(l, 0) 


= 1 - 1 = 0 . 


The desired circulation of the vector field is zero. 

Example 7, Compute the circulation of the vector 
field 


if 


a = xyi + yz\ + xzk 


L: 

Solution. We have 


^+y*=l, 

x + y + z = 1. 


C = ^ (a, dr) = j xy dx + yz dy + xz dz. 


The curve L is an ellipse obtained by cutting the cylinder 
x 1 + y 2 = 1 by the plane x + y + z = 1. Let us find 
the parametric, equations of the curve. The projection 
of any point of the curve on the xy-plane lies on the 
circle a? + y 2 = 1. From this we obtain x = cos t, 
y = sin t. But the ellipse lies in the plane x + y + z = 1, 
whence z = 1 — x — y or z = 1 — cos t — sin t. Thus 
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the parametric equations of the curve L are: 


x = cos/, 

y = sin /, i (0</<2ji). 

z = 1 — cost — sin/ J 

From this we find 

dx = —sin t dt, dy = cos t dt, dz = (sin / — cos t)dt 
and so the circulation is 

C— j ( — cosfsin 2 f + sin/(l— cosf —sin/)cos/ 

+ cos/(l —cos/ —sin/) (sin/ —cos t)\dt 

2n 

= f (—sill 2 /cos/ + sill 2/ — cos 2 /sin/—cos 2 / + cos 3 /) dt 


Zn 

— — j cos 2 tdt — — ji. 


Compute the circulation C of a vector a around the 
given curve L : 

175. a = (xz + y) i + {yz — x) j — (r 5 + y 2 ) k; 


L: 


fx 2 + P 2 = l, 
( 2 = 3. 


176. a = y 2 i + z 2 j + x 2 k; 


L: 


\ x 2 + y 2 + z 2 = R*, 
[ x 2 + p 2 = Rx 


(*>0). 


177. a = (2x + z) i + (2y — z) j + xyzk. L is the line of 
intersection of the paraboloid of revolution x 2 + y 2 — 
= 1 — z with the coordinate planes. 

178. Show that if in a vector field the circulation of a vec¬ 
tor around any closed circuit is equal to zero, then there 
can be no closed vector lines in such a field. 
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Sec. 15. The curl (rotation) 
of a vector field 

Suppose we have the vector field 
a (M) = P (x, y, z) 1 + Q (x, y, z) j + R (x, y, z) k. 
We will assume that the coordinates P , Q, R of the vector 
a (Af) are continuous and have continuous partial deriv¬ 
atives of the first order with respect to all its arguments. 

Definition 1. The curl (or rotation) of a vector a ( M) 
is a vector (symbolized: curl a ( M) or rot a ( M )) defined 
by the equation 

J_(9Q 


\ dx 


or, in easy-to-remember symbolic form, 

j 


( 1 ) 


( 2 ) 


This determinant is ordinarily expanded in terms of 
elements of the first row, in which case the operations of 
multiplication of elements of the second row by elements 
of the third row are regarded as operations of differen¬ 
tiation; for example, 

. *Q 


dx 


dx 


Definition 2. If in some region G we have curl a = 0, 
then the field of the vector a in G is said to be irrotational. 
Example 1. Find the curl of the vector 

a = (x + z) i + (y + z) j + (x 2 + z) k. 
Solution. Using formula (2), we obtain 
i i k 


curl a = 


dx dy dz 

X + Z y + Z x 2 + z 
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Expanding the determinant in terms of elements of the 
first row and regarding the operation of multiplication, 
say dldy by x 2 + z, as an operation of partial differen¬ 
tiation, we obtain 

curl a = —i— (2x — 1) j. 

Example 2. Find the curl of H, the intensity vector of 
a magnetic field under the conditions of example 3 of 
Sec. 10. 

Solution. The magnetic-field intensity vector H is 

H = f d.rj 

or 


i i k 

H = 4- 0 0/ 

p 2 

x y z 

where p 2 = x 2 -f y 2 , whence, by (2), 

i £ £|-[£b 

2 1>J 21 x n 

**+y* x*+y* 




2 lx \ 

+ !/*> 


(x 2 +y 2 ) 2 
(x* + y 2 ^ 0). 


2 + y*-2 y 2 |._ n 

(* 2 + ^)* J k ° 


Thus, curl H = 0 everywhere except the z-axis, at the 
points of which the last formulas are meaningless (the 
denominator vanishes), that is, the field of the vector H 
is irrotational everywhere outside the points of the z-axis. 
Find the curl of the following vectors: 

179. a = (x 2 + y 2 ) i + (y 2 + z 2 ) j + (z 2 + x 2 ) k. 

180. a = z s i + j/ 3 } + x 3 k. 

181. a = 1 (—y 2 i + x?j). 

182. Show that if the coordinates of the vector a ( M) 
have continuous partial derivatives of the second order, 
then 


div curl a = 0, 
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that is, the field of the vector curl a ( M) is a solenoidal 
field. 

183. Show that 

(a) curl (a ± b) = curl a ± curl b, 

(b) curl (Xa) = X curl a. 
where X is a numerical constant. 

184. Show that if u = u ( M ) is a scalar function and 
a = a ( M) is a vector, then 

curl (ua) = u curl a + [grad u, a). 

185. Show that if a and b are constant vectors, and r 
is the radius vector of the point M ( x , y, z), then 

curl (r, a) b = (a, bl. 

186. Show that 

curl (ra) = —(r, a|, 

wher e a is a constant vector and r = | r | = 
= V x? + y 2 + z 2 . 

187. Show that curl (/ (r) a) = ^ Ir, a), where / (r) 
is an arbitrary differentiable function of its argument and 
a is a constant vector. 

188. Show that the vector field a = / (r) r is irrotational 
that is, curl a = 0. 

189. Show that 

div [a, bl = (b, curl a) — (a, curl b). 

190. Show that the curl of a field of linear velocities v of 
a rotating solid is a constant vector in the direction paral¬ 
lel to the axis of rotation, the modulus of which is equal 
to twice the angular velocity of rotation: curl v = 2o>. 

191. Determine the angular velocity o> of rotation of 
a solid about a fixed axis passing through some point of 
the solid if its linear velocity is 

v = 2ii + + xzk. 

192. Show that the field of the curl of the vector a ( M) 
is free of sources and sinks. 

193. What must the function / (x, z) be so that the curl 
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of the vector field 

a = yzi + / (x, z) j + xyk 

is coincident with the vector k — i? 


Sec. 16. Stokes' theorem 
Suppose the coordinates of a vector 
a (M) = P (x , y, z) i + Q (x, y, z) j + R (x, y, z) k 

are continuous and have continuous partial derivatives. 

Theorem. The circulation of 
the vector a around a closed con¬ 
tour L is equal to the flux of the 
curl of the vector through any 
surface 2 spanning the con¬ 
tour L: 

§ (a, dr) = j j (curl a, n°) da. 


( 1 ) 

It is assumed that the orien¬ 
tation of the normal n° to the 
surface 2 is matched with the 
orientation of the contour L 
so that, when viewed from the 
end of the normal, the contour 
is traversed in the chosen direction counterclockwise. 

Example 1. Compute the circulation of the vector 

a = yi + x*j — zk 



around the contour L : 


x 2 + i/ 2 = 4, 
i z = 3 


(1) directly and 


(2) via the Stokes theorem. 

Solution. (1) The contour L is a circle of radius R = 2 
lying in the plane z = 3 (Fig. 32). We choose the orienta¬ 
tion as shown in the drawing. The parametric equations of 
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the curve L are 

x = 2cos<, \ 

y = 2 sin t, i (0<*<2n), 

2 = 3 J 

so that 

dx = —2 sin t dt , dy = 2 cos t dt, dz = 0. 

For the circulation of the vector a we have 

2 ji 

C= j [2sin<( — 2sinf) + 4cos 2 <2cos< — 3-0)rf/ = — 4n. 
o 

(2) To compute the circulation via the Stokes theorem, 
choose some surface £ spanning the contour L. For £ 
it is natural to take a circle having L as its boundary. 
According to the chosen orientation of the contour the 
normal n u to the circle has to be taken equal to k: n° = k. 
Then 

I i J k I 


curl a = 


dy dz 
X 2 —Z 

Therefore, by the Stokes theorem, 

C= j j (curl a, n°) do = ( j (2x-\)do 


(2x— 1) k. 


■H 


(2pcosq>— 1) pdp= — 2ji-^-| 2 = — An. 


194. Show that the flux of the rotor through an open 
surface spanning a given contour does not depend on the 
shape of the surface. 

Find the circulations of the vectors around the indi¬ 
cated contours (1) directly and (2) via the Stokes theorem. 

f x 2 + y 2 = 4 


= zi + xj + yk; L: 


196. a = yi — xj + zk; L: 


1 2 = 0. 

fx 2 +y 2 -fz 2 = 4, 
j x 2 + y 2 = z 2 (z>0). 
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197. a = 2xzi — yj + zk around a contour formed by 
intersection of the plane x + y + 2z = 2 with the coor¬ 
dinate planes. 

<z = x 2 + y 2 , 

198. a = yi — zj + Oc+jOk; L: | z==1 


199. 

200 . 


fx 2 + J/ 2 i-z 2 = 16 I 
a — z i, L. | x = o i j / = o,z = 0. 

[x = y 2 + z 2 , 

a = zy 2 i + xz 2 j -f x 2 yk; L: j g 


i x 2 "|- y 2 = 9 1 
2M ' W + 4,=5. 

fx 2 -f y 2 + z 2 =\, 
202. a = yi-x) + zk; L: | = z 


203. Given a vector field of velocities v of the points 
of a solid rotating with a constant angular velocity co 
about the z-axis. Calculate the circulation of the field 
around the circle 


{ x = ocosf, 

y = asint, (0<*<2n), 
z =0 


directly and via the Stokes theorem. 

From the Stokes theorem we find that the projection 
of the vector curl a on any direction □ is independent of 
the choice of the system of coordinates and is equal to 
the surface density of the circulation of the vector a 
around the contour of the area perpendicular to that 
direction: 

§ (a, dr) 

pr„ curl a| M = (curla,n°)| = lim —. (2) 

|*f <E)-M s 

Here, (2) is a plane area perpendicular to the vector n; 
S is the area of (2); L is the contour of the area and is 
oriented so that the traverse of the contour is counter¬ 
clockwise as seen from the end of the vector n; the nota- 


8-910 
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tion (2) M means that the area (2) contracts to the 
point M at which we consider the vector curTSTand the 
direction of the normal n to that area is always the same. 

Example 2. Compute the density of the circulation of 
the vector a = yi around the circle 


{ x = acost, 

y = asint, (0^<<2ji), 

2—0 

at the centre of the circle, M (0, 0, 0), in the positive 
direction of the z-axis. 

Solution. Here, (2) is a circle of radius a with centre 
at Af, so that S = n a 2 . 

The desired density of the circulation is 





On the other hand, 



and, hence, 


curl a = 


i 

d 

dx 


y 


j 

j_ 

dy 

0 


k 

d 

dt 

0 


-k 


(curl a. n°) \„ = (-k, k) - -1, 


which, by virtue of (2), confirms the correctness of the 
result. 

204. Compute the density of the circulation of the vector 
a = zl + sj 4- yk around the circle L: {y = a cos t, 
z = a sin t, x = 0, (0 ^ t < 2n)} at the centre of the 
circle, C (0, 0, 0), in the positive direction of the x-axis. 

205. Compute the density of the circulation of the vector 
a = 2yi + 5xzj around the ellipse L : {x = a cos t, y = 
= b sin t, z = 1 (0 ^ t < 2ji)} at the centre C (0, 0, i) 
of the ellipse in the positive direction of the z-axis. 
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Sec. 17. The independence 
of a line integral of the path 
of integration. Green's formula 

Definition. A region G of three-dimensional space is 
said to be simply connected (more precisely* it is a simply 
connected plane region) if any closed contour lying in the 
region can be spanned by, a surface lying entirely in G. 
For example, the whole of .three-dimensional space and 
the interior of a sphere are simply connected regions; the 
interior of a torus and three-dimensional space with 
a straight line deleted are not simply connected regions. 
Theorem. In order that the line integral 

j (a, dr) = j P (x, y, z)dx + Q (x, y, z)dy + R (x, y, z) dz 

should not depend on the form of the path of integration L, 
it is necessary and sufficient that the vector field 
a = P (x, y, z) i + Q (x, y, z) j + R (x, y, z) k 
be irrotational, that is, 

curl a (M) = 0. 

It is assumed here that the coordinates P (x, y, z), 
Q (x, y, z), R (x, y, z) of the vector a have continuous 
partial derivatives of the first order and the domain of 
definition of the vector a (M) is simply connected.- 
In that case, the line integral j (a, dr) will depend 

L 

solely on the position of the initial and terminal points 
of the path of integration L. 

If the theorem is complied with, the circulation of 
the vector a (Af) around any closed contour C located in 
the field of the vector a (M) is equal to zero: 

^(a, dr) = 0. 
c 

Example 1. Show that in the field of the vector 
a = xyh i + x?yzj + -^xV k 



S-5TS- 
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the line integral j (a, dr) is independent of the shape 

i 

of the path of integration L. 

Solution. The coordinates of the vector a are everywhere 
continuous functions so that the domain of definition G of 
the vector a is the entire space (a simply connected region). 
In this region we have 


curl a = 


i i k 

a a _a_ 

dx dy dz 

xy 2 z x 2 yz 


0 . 


Consequently, the line integral 

j (a, dr) = j xy 2 z dx + x*yz dp + y x 2 y 2 dz 

is independent of the shape of the path of integration L. 
In particular, for the plane vector field 


a (M) = P (x, y) i + Q (x, y) j (1) 

we have 


curl a (M) = 


i j k 

a a a 

dx dy dz 

P Q 0 



Therefore, for the plane vector field (1) that is defined 
in a simply connected region G, the condition curl a (M) = 
= 0 is written in coordinate form thus: dPIdy = dQ/dx. 
This means that in order for a line integral 

j^P(x, y)dx + Q(x, y) dy 

a plane field defined in a simply connected region G to 
independent of the shape of the path of integration, it 
necessary and sufficient that the following relation hold : 


dP _ dQ 
dy dx • 

Remark. The requirement that the region G, where the 
vector a = a (AT) is defined, be simply connected is essen- 




Sec. 17] 


Green’s formula 


117 


tial. If the region G is nonsimply connected, then, provid- 
e d t hat curl a (A/) = 0, the line integral may depend 
on the shape of the path of integration. 

Example 2. Let us consider the line integral 
C — ydx . xdy 
J x*+y*^x*+y*‘ 


Solution. The integrand is meaningless at the point 
0 (0, 0); and so we eliminate that point. In the remaining 
part of the plane (which is then a nonsimply connected 
region), the coefficients of dx and dy are continuous and 
have continuous partial derivatives, and the following 
identity holds true: 

17 (i*+7») “IF ( — i*T7*) • 

On the other hand, if we compute the integral around the 
circle L\ x 2 + y 2 = R 2 , then by parametrizing the equa¬ 
tion of the circle we get 

& dy = T* - 9ip, -i + i?,cos -^dt = *\dt = 2n. 

L 0 6 


We find that the circulation is nonzero and, hence, the 
line integral depends on the path of integration. 

Determine in which of the vector fields indicated below 
the integral is independent of the shape of the path of 
integration: 

206. a = z*i + x*j + j/ 2 k. 

xl + yj + zk 


207. a = 


/l + i* + y*+ 2 i 
yi-ij + sk 


Green's formula 

Given in a region D with boundary L a plane vector 
field 

a = P (x, y) 1 + Q (x, y) j, 

where the coordinates P (x, y), Q ( x , y) are continuous 
and have continuous partial derivatives dPIdy, dQldx. 
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Then Green’s formula, 

§Pdx + Qdy = j j(-g— ^dxdy, (2) 

L D 

holds true. Here, the boundary L is traversed so that the 
region D is on the left. ” 



Fig. 33 

The region D may also be nonsimply connected so that 
the boundary may consist of several components (Fig. 33). 
In that case, the integral " '•> 

j P dx + Q dy 

is understood to mean the sum of the integrals over all 
components of the boundary of D. 

Green’s formula (2) is a special case of the Stokes theo¬ 
rem (see Sec. 16). 

In certain cases, Green’s formula permits simplifying 
the computation of the circulation of a vector field. 
Example 3. Compute the circulation of the vector 
a = l/l + * 2 + y 2 1 + y [xy + In (r+/l-|-a: 2 -t-{/ 2 )] j 

around the circle i 2 + y 2 = R 2 . 

Solution. The circulation of the given vector is 

C = £(a, dr) = § \fl + x 2 + y* dx 


+ y [xy + ln(a;+ K1-f-r 2 + y 2 )]dy. 
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Here 

P = 1^1 + x 2 -(- y 2 , Q = xy i + y\n(x + Vl+x Zj cy i )- 

We find the partial derivatives 

9P_ _ v dQ_ _ 1)2J _ y 

dV Vi + x'+y*' ** y ' V V\ + x* + y* ' 

Using Green's formula, we obtain 

C= f f (y 2 H— ■ y — - y - ) dxdy 

V ' /l + ** + »* /l + X* + »« / 

= j j y*dxdy. 

Passing to polar coordinates, 

i = p cos qp, y = p sin q>, 

we have 

C — j j p 2 sin 2 cpp dp dq> = j j p 3 sin 2 <p dp dqp. 

D D 

Since 0<<p<2n,0<p<i?, it follows that 



Compute the circulation of the following vectors around 
the given contours using Green’s formula: 

209. a = (y + x) i + (y — x) j; L: x + y = 1, x = 0, 
y = 0. 

210. a = (x — y 3 ) i + 2 xy j; L: y = x, y = x 3 . 

211. a = x In (1 + y 2 ) i + ^ j; L: x 3 + y 3 = 2x. 

212. a = y 3 l — x*j; L: x + y = -1, x = 0, y = 0. 

, 4a . (3x-y3/i +x » + 4y»)l + (l8y»+x» /l+x»+4y»)l . 

3/l+x*+4»3 

L: x 2 + y 2 =l. 

214. Use Green’s formula to compute the difference be¬ 
tween the integrals 

J i= } (x + y) 2 dx — (x—y) z dy 

AmB 
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and 

/ 2 — j (x + y) 2 dx—(x — y) 2 dy, 

AnB 

where AmB is a line segment joining the points A (0, 0) 
and B (1, 1), and AnB is an arc of the parabola y = i 2 . 

215. Prove that the integral 

§(2x + y)dx+2xdy, 

where L is a closed contour, yields the area of the region 
bounded by that contour. 

216. Using Green's formula, compute the line integral 
j (a, dr) in the vector field 

a = ( e x sin y — y) I 4- (c 1 cos y — 1) j, 

where the curve L is the upper semicircle x 2 + y 1 — 2x 
traversed in the direction from the point A (2, 0) to the 
point 0 (0, 0). 
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Sec. 18. The criterion for the potentiality 
of a vector field 

Definition. A vector field 

a (M) = P (x, y, i) i + Q (x, y, z) j + R (x, y, z) k 
specified in a spatial region V is said to be potential if 
there exists a scalar function q> ( M) such that at all 
points of V the following equality holds: 

a (M) = grad <p ( M). (1) 

The function <p ( M ) = <p (x, y, z) that satisfies (1) in 
V is termed the potential (or the potential function) of 
the vector field a. 

The relation (1) is equivalent to the following three 
scalar equations: 

*) = -§-, Q(x,y,z) = - g-, R(x,y,z) = ZL. 

( 2 ) 

The potential of a field is not defined uniquely but 
only up to an additive constant. 

Remark. For force fields, the function <p (Af) is ordi¬ 
narily called a force function, and the potential is the func¬ 
tion —q> ( M). 

Example 1. (The electrostatic field of a point charge.) 
Show that the field of electric intensity E set up by a point 
charge q located at the coordinate origin, 

, E = -£r, r=]/x 2 + y 2 + z 2 , 

is a potential field. 
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Solution. The problem is posed thus: show that there 
exists a function <p (x, y, z) such that relations ( 2 ) hold. 
In our case, we have 

P(x,y,z) = ^, Q(x,y,z) = %, R(x,y,z) = -g-. 
Since 



and, analogously, 



it follows that the function 


1 dr 
~r*~dx 


<p(s, y, Z) = ~ 7 = ~ , 

is the potential of the given field: 

grad ( — 7 ) - E. 

In this case, the coordinate origin (this is where the 
charge q is located) is a singular point of the field E. 

Theorem. For a vector field a (Af) specified in a simply 
connected region V to be potential, it is necessary and suf¬ 
ficient that the following condition hold at every point of V: 

curl a = 0. (3) 


In other words, for a vector field specified in a simply 
connected region to be potential, it is necessary and suf¬ 
ficient that it be irrotational. 

The potential q> (a:, y, z) of the vector field 
a = P (x, y, z) i + Q (x, y, z) j + R ( x , y, z) k 
is defined by the formula 

(*. v, 2) 

9 (x,y,z)= j Pdx + Qdy + Rdz, (4) 

(*.. Vo. 2,) 

where ( x 0 , y 0 , z 0 ) is some fixed point of the field and 
(x, y, z) is an arbitrary current point. 

Example 2. Show that the field of the vector 
a = i 2 ! + + z 2 k 

is a potential field. 
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Solution. The coordinates P=s?,Q = if,R = z*ot 
the vector a are infinitely differentiable functions through¬ 
out the space so that a is an infinitely differentiable 
vector defined throughout three-dimensional space. We 
have 


i j k 


dx dy dz 
X 2 y 2 Z 2 






By virtue of the theorem given on p. 122 the field of the 
vector a is a potential field. It is readily seen that the 
function 


?(».».«) - “+V+- + c. 

where C is an arbitrary constant, is the potential of the 
given field. 

Check to see whether the following vector fields are 
potential fields: 

217. a = xzi 4- 2yj + xyk. 

218. a = (2z y + z 2 ) i + (2 yz + **) j 4 - (2zz + y 2 ) k. 

219. a = 1 (z»i + y s j + xz*k). 

220 . a = yz cos xy-i + zz cos xy -j -f- sin xy-k. 

221 . a = In (1 + z 2 ) i + In (1 + z 3 ) j + zzk. 

222 ..= (^ + i)l + (i + i), + (i+±)k. 

223. H = -^-(-pi + zj), r 2 = z 2 + y 2 , r^O. 

224. Prove that the field a = /(r)-r, where / (r) is a 
differentiable function, is a potential field. 

225. Show that the vector lines in the potential field 
a = grad <p are perpendicular to the level surface of 
the function cp. 
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Sec. 19. Computing a line integral 
in a potential field 

Theorem. A line integral in a potential field a ( M) is 
equal to the difference between the values of the potential 
<p (M) of the field at the terminal and initial points of 
the path of integration : 

M, 

j(a,dr) = q>(M 2 )-<p(M 1 ). (1) 

if, 


Example 1. Compute the line integral in the field 
of the vector 

r = xl + yj + zk 

along the line segment bounded by the points A/j (—1, 0, 
3) and M 2 (2, -1, 0). 

Solution. We will show that the field of the given vector 
is a potential field. Indeed, 


curl a = 


1 j k 
J_ 9_ J_ 
dx dy dz 


x y z 


0 . 


It is easy to see that the potential of this field is 


(p (x, y, z) = x *+y*+ 2 * -|-C. 
Using formula (1), we obtain 


M, 

\ (a, dr) = <p(2, -1, 0)-<p{- 1, 0, 3)=|-5=-|. 

M, 

Note that it is immaterial what line joins the points M x 
and M t ; for fixed M 1 and the integral 


j (a, dr)= j xdx + ydy + z 


always has the same value. 
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Computing the potential of a field in Cartesian coordinates 
The formula 
9 (x, y, z) 

(*. v. z) 

= j P(x, y, z)dx + Q{x, y, z)dy+ R(x, y, z)dz 

(*., Vo. *o) 

( 2 ) 

may be used to find the potential function <p ( M) = 
= <p (x, y, z) of a specified potential field 

a (x, y, z) — P (x, y, z) I + Q (x, y, z) j 

+ R (x, y , z) k. 

To do this, fix the initial point Af 0 (x 0 , y 0 , z 0 ) and join 
it with the current point M (x, y, z) with a broken line 



MqABM whose segments are parallel to the coordinate 
axes; namely, M^A || Ox, AB || Oy, BM || Oz (Fig. 34). 
Then formula (2) takes the form 
* v 

<p (x, y, z) = j P (x, y 0 , z 0 ) dx+ (x, y, z 0 ) dy 
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+ j R (x, y , z) dz, (3) 

where x, y, z are the coordinates of the current point on 
the segments of the broken line along which the integra¬ 
tion is performed. 

Example 2. Prove that the vector field 

a = (y -f z) i + (* + z) j + (x + y) k 


is a potential field and find its potential. 

Solution. 1st method. A necessary and sufficient condi¬ 
tion for the potentiality of a field a (M) is that curl 
a (M) be zero. In our case, 
i j k 

outl *-li x i 

y+zx+zx+y 


= (l-l)i + (l-l)i + (l-l)k^.O. 


That is, the field is a potential field. The potential of 
the field can be found with the aid of formula (3). For 
the initial fixed point we take the coordinate origin 
0 (0, 0, 0). We thus have 

<p(x, y, z)= j (0 + 0 )dx+ j(x + 0 )dy + j(j + y)dz 
0 0 0 

= xy + xz-f yz. 

To summarize, 

<p (x, y, z) = xy + xz + yz + C, 
where C is an arbitrary constant. 

2nd method. By definition, the potential <p (x, y, z) 
is a scalar function for which grad 9 = a. This vector 
equality is equivalent to three scalar equalities: 


■&=»+*• 

(4) 

W = x+Z ' 

(5) 


(6) 
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Integrating (4) with respect to x, we obtain 

<p(x, y, z)= J (y + z)dx = xy + xz + f (y, z), (7) 

where / ( y , z) is an arbitrary differentiable function of 
y and z. Differentiating both sides of (7) with respect to 
y and taking into account (5), we obtain a relation for 
finding the as yet undetermined function / ( y , z). We 


have 

*L= X +°L 

dy “ dy 


or 

I+2 ” I+ f ■ 


whence 



-*• 

(8) 

Integrating (8) with respect to y, we have 



f(y, z)=^zdy = zy + F(z), 

(9) 


o 


where F (z) is an as yet undetermined function of z. 
Substituting (9) into (7), we get 

<p {x, y, z) = xy + xz + zy + F (z). 

Differentiating this equation with respect to z and taking 
into account (6), we obtain an equation for finding F (z): 

x+y = x+y + ^ r , 

whence dFIdz = 0 so that F (z) = C = constant. 

Thus we have 

<p (x, y, z) = xy + yz + zx + C. 

3rd method. By the definition of the complete differen¬ 
tial of the function q> (x, y, z) we have 

j. i d( P j.. . j. 
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Substituting in place of the partial derivatives dyldx, 
dtp/dy, dtpldz their expressions taken from (4), (5), (6), 
we obtain 

dq> = (y + z) dx + (i -f z) dy + (x + y) dz 


or, after some simple algebra, 

d<p = (y dx + x dy) + (z dx -f x dz) + {y dz + z dy) 

= d (xy) + d (xz) + d (yz ) = d (xy + xz + yz). 

Thus 

dtp = d (xy + yz + zx), 
whence it follows that 


<p ( x , y, z) = xy + yz + zx + C. 


In the following problems, establish the potentiality of 
the given vector fields a (Af) and find their potentials 

226. a = 2xyzi + x 2 zj -f x 2 yk. 

227. a = (yz + 1) 1 -f- xz} + xyk. 

228. a = (2xy + z) i + (x 2 — 2y) j + xk. 

1 + 1+k 

'x + y+z ■ 


229. < 


230. a= yri +“i+y k . 

1+*Vz* 

231. a = y. 

232. a = -V. 


233. a = r-r. 

When the region Q is a star with centre at the coordi¬ 
nate origin 0(0, 0, 0)*, the potential <p(Af)of some 


* The region Q is called a star-shaped region relative to the 
point 0 belonging to Q. if any ray e manatin g from that point cuts 
the boundary of tne region at one point at most. For example, star¬ 
shaped regions in the plane are the plane itself, a parallelogram, 
a circle; in three-dimensional space, the space itself, a parallele¬ 
piped, a sphere. 
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vector field a = a (Af) at the point Af (x, y, z) may 
be found from the formula 

qp(Af) = j(a (Af'), r (M))dt + C, C = constant, (10) 
u 

where r (Af) = xi + y\ + zk is the radius vector of the 
point Af ( x, y, z), and the point Af' ( tx, ty, tz) for 0 ^ t 
< 1 runs over the segment OM of the straight line passing 
through the points 0 and Af. 

Example 3. Find the potential of the vector field 
a = yzi + xz\ + xyk. 

Solution. It is readily seen that curl a = 0, which 
means the given vector field is a potential field. This 
field is defined throughout three-dimensional space and 
is star-shaped with centre at the coordinate origin 
0 ( 0 , 0, 0 ), and therefore to find the potential we take 
advantage of formula (10). Since in this case 

a (Af') = a (tx, ty, tz) = t 2 yzi + t 2 xz\ + t*xyk, 
it follows that the scalar product of the vectors a(Af') 
and r (Af) is equal to 

(a (Af'), r (Af)) = t 2 ( xyz + xyz + xyz) — ZtHyz. 

The desired potential is 

l l 

<p(Af)= j(a(Af'), r (M))dt — xyz^'3t 2 dt + C = xyz + C. 
o o 

Thus, 

(p (Af) = xyz + C. 

Using formula (10), find the potentials of the following 
vector fields: 

234. a = ai + Pj + yk, where a, fi, y are constants. 

235. a = (y + z) 1 + (x + z) j + (y -f x) k. 

236. a = yi + xj + e l k. 

237. a = e* sin p-i + e* cos y- j + k. 
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Sec. 20. The Hamiltonian operator del 

Many operations of vector analysis may be written in 
abbreviated form and in a form convenient for calcula¬ 
tions; this is done through the use of a symbolic operator 
called the Hamiltonian operator del: 

This operator combines both differential and vectorial 
properties. We will regard the formal multiplication of 
d/dx by the function / (x, y, z) as the partial differen¬ 
tiation df Idx. 

Within the framework of vector algebra, we will per¬ 
form the formal operations involving del as if it were 
a vector. Using this formalism, we obtain the following. 

1. If u = u (x, y, z) is a scalar differentiable function, 
then, by the rule of multiplying a vector by a scalar, we 
have 


„ . d 

u== \'i7 




du . . du . . du 

e7+' ir +k ir= 


2. If a = P (x, y, z) i + Q (x, y, z) j + R (x, y, z) k, 
where P, Q, R are differentiable functions, then, by 
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the familiar formula for a scalar product, we have 

<*••)-(*£+i- w+*-k’ «+«+ flk ) 

in particular, (V, c) = 0, where c is a constant vector. 

3. If a = P (x, y, z) i + Q (x, y, z) j + R ( x , y, z)k, 
then 

i j k 

aj= -L-L-L = curl a; (A) 

P Q R 

in particular (V, c) = 0, where c is a constant vector. 

Continuing the formalism of operations with V as 
a vector, we obtain the following from the distributive 
property of scalar and vector products: 

(V, a + b) = (V, a) + (V, b), (5) 

that is, div (a + b) = div a + div b 

(V, a + b] = [V, al + [y, b) (6) 

or curl (a + b) = curl a -f curl b. 

Formulas (5) and (6) may also be interpreted as an 
exhibition of differential properties of the del operator 
(V is a linear differential operator). 

When using the formalism of operations involving the 
del operator regarded as a vector, one must bear in mind 
that del is not a vector, for it has neither magnitude 
nor direction, so that, for example, the vector [V, a] 
will not, generally speaking, be perpendicular to the 
vector a (however, for the plane field a = P (*, y) i + 
+ Q ( x > y) j the vector 

IV.aJ-curl •-(-!—f-)k 

will be perpendicular to the zy-plane and, hence, to the 
vector a). In the same way, the concept of collinearity is 
meaningless with respect to the symbolic vector V- For 
example, the expression [v<p, V^], where (p and ip are 
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scalar functions, formally resembles a vector product of 
two collinear vectors, which product is always equal to 
zero. But this is not true in the general case. Indeed, the 
vector Vq> = grad <p is directed along the normal to the 
level surface <p — constant, while the vector Vi|> = grad i|) 
defines the normal to the level surface t|> = constant. 



and in the general case these normals are not necessarily 
collinear (Fig. 35). On the other hand, in any differentiable 
scalar field cp we have Iv<p, V<pl = 0. These examples 
show that the del operator must be handled with care. 

Besides its vectorial nature, the Hamiltonian operator 
del has a differential aspect to it. Taking into account 
the differential aspect of v, we will agree that the operator 
V acts on all quantities that follow it. In this sense, 
(V, a) (a, v). Indeed, 

(V, a) = div a, 

whereas 

is a scalar differential operator. 
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When applying the del operator to a product of any 
quantities, one must bear in mind the rule for differenti¬ 
ating a product: 


~(uv) = 


du 

~dx 


From this it follows that the del operator must be applied 
in succession to each factor, leaving the other factors 
unchanged, and then the sum of the resulting expressions 
taken. In this procedure we are guided by the following 
rules. 

i°. If the del operator acts on some product, first 
take into account its differential nature and only then 
its vectorial properties. 

2°. In order to note the fact that the del operator 
does not act on some quantity involved in a complicated 
formula, that quantity is labelled with the subscript c 
(constant),' which may be removed in the final result. 

3°. All quantities not acted upon by the del operator 
are placed in front of the operator (that is, to the left) 
in the final result. 

Example 1. Show that 

div (ua) = u div a 4- (a, grad u). 

Here, u is a scalar function and a is a vector function. 

Solution. In symbolic notation, 

div (ua) = (v, ua). 

Taking into account first the differential nature of V, 
we have to write 

(V, ua) = (V, u c a) + (V, ua c ). 

Regarding the expression (V, u c a), we can take the con¬ 
stant factor u c outside the del sign and, as a scalar, outside 
the sign of the scalar product; this yields 

(V, u c a) = (u r v, a) = u c (V, a) = u (v, a). 

(The last step is to drop the subscript c). 

In the expression (V, ua c ), the del operator acts only 
on the scalar function u, and so we can write 


(V, ua c ) = (yu, a c ) = (a c , Vu) = (a, Vu) 
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to get the formula 

(V, ua) = u (V, a) + (a, Vu) 


div (ua) = u div a + (a, grad u). 
Example 2. Show that 

curl (ua) = u curl a — la, grad ul. 
Solution. In symbolic notation, 

curl (ua) = (V, ua). 


Taking into account the differential properties of V, we 
first write 

(V, ua] = [V, u c al + (V, ua c l. (7) 

Then in the first term on the right we take the scalar 
factor u c outside the del sign and outside the sign of the 
vector product, which yields 

(V, u c a] = u e (V, a] = u tv, al. 


In the second term in (7) we refer u to the operator V and 
change the order of the factors so that the vector a c , which 
del does not act on, is in front of V. This yields 
lv, ua c ] = I vu, a c ] = —(a, Vul. 

Thus 

[V, ua] = u [v, a] — («, Vu] 
or 

curl (ua] = u curl a — la, grad u]. 

Example 3. Use the symbolic method to find div [a, b] 
Solution. We have 

div [a, b] = (V, (a, b]) = (V, (a, b c ]) 

+ (v, la c , b]). (8) 


Using the property of cyclic permutation of factors in 
a mixed product, we transform the expression on the 
right of (8) so that all constant quantities are in front 
of the del operator and all variable quantities follow it. 
This yields 

div [a, b] = ([V, a], b c )-(V. [b, a c |) 


-(IV, a], b e )-([V, bj. a c )-(b. [V, a])-(a, [V, b]) 
or div la, b] = (b, curia) —(a, curlb). 
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Remark. The use of the symbolic method enables us to 
avoid cumbersome analytical transformations and obtain 
the final results very quickly. On the other hand, the 
various formal transformations involving the del opera¬ 
tor must be performed with extreme caution in order to 
avoid serious mistakes. For this reason, if there is any 
doubt about the final result, it is wise to verify it by 


the analytical method. 

238. Show that 


(a) 

_ / u \ vVu — uW 

V (J = v. S 

(b) 

vf (u) = /' (u) vu. 

239. Prove that 

the vector [Vu, Vt>l is solenoidal if 


u and v are differentiable scalar functions. 

Use the Hamiltonian operator del (V) to prove the 
following equations: 

240. (a) grad ( uv) = v grad grad v; 

(b) curl (a, b) = (b, V) a — (a, V) b + a div b — 

— b div a. 

241. curl [a, r] = 2a, where a is a constant vector. 

242. Prove that the vector a = u grad v is orthogonal to 
curl a. 


Sec. 21. Second-order differential 
operations. The Laplace operator 

Second-order differential operations are obtained as 
a result of a twofold application of the operator V to 
fields. 

Suppose we have a scalar field u = u (Af). In this 
field, the operator V generates a vector field Vu = grad u. 

In the vector field Vu, the operator V, when applied 
a second time to Vu, yields the scalar field 

(V, Vk) = div grad u (1) 

and the vector field 

(V, Vu) = curl grad u. (2) 

If a vector field a = a ( M) is given, the operator V 
generates in it a scalar field (V, a) = div a. In the scalar 
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ield div a the operator V generates a vector field 

V (V, a) = grad div a. (3) 

In the vector field a = a ( M ), the operator V also gener¬ 
ates a vector field lv, a] = curl a. If the operator V is 
again applied to this field, we obtain the scalar field 
(V, (V, a]) = div curl a (4) 

and the vector field 

lV» lv, all = curl curl a. (5) 

The formulas (1) to (5) define what are called differen¬ 
tial operations of the second order. 

Example 1. Suppose a function u = u (x, y, z) has 
continuous partial derivatives up to second order inclu¬ 
sive. Prove that 

curl grad u = 0. 

Solution. 1st method. Operating formally, we obtain 
curl grad u = (v, V«1 = lv, vl u = 0 

since lv, Vl = 0 being the vector product of two identical 
“vectors”. 

2nd method. Using the expressions for the gradient and 
curl in Cartesian coordinates and taking into account the 
given conditions, we have 


du du du 

— I d* u d * u \ • i / d*u d*u \ . . ( d*u d*u \ , _ ^ 
~\dydz dzdy) l '\dzdx dx dz)^~'\dxdy dydxf ~ ’ 

since the mixed derivatives are equal in this case. 



Sec. 21] 


Second-order differential operation> 


137 


In similar fashion proof is given that for the vector 
field 

a = P (x, y, z) i + Q (x , y, z) j + R (x, y, z) k, 
the coordinates P, Q, R of which have continuous partial 
derivatives of the second order, we obtain div curl a = 0. 

Note particularly the second-order differential opera¬ 
tion div grad u = (V, Vu). Assuming that the function 
u (x, y , z) has second partial derivatives with respect to 
x, y, z, we obtain 


(V, V a) = ( 1F I + ^i + ^k, _i + _j + _k 

d*u . a»u d*u 


8 l 9u\ . 8 / 8u \ . 8 ( 8u\ 8*u . d*u . d*u _. 

Thus, (V, Vu) = Au, where the symbol 


Au 


8 l u.8 t u 8*u 

dz*~*~ dy t+ di* 


is termed the Laplace operator (or Laplacian). It may be 
represented as a scalar product of the Hamiltonian oper¬ 
ator V into itself, that is, 

A«(V, V>-*-£ + $ + £. 

This operator plays an important role in mathematical 
physics. 

Let us examine another second-order operation curl 
curl a. We have curl curl a = [V, (V, a]]. Let us take 
advantage of the formula for a vector triple product writ¬ 
ten as [A, [B, Cll = B (A, C) - (A, B) C. Replacing 
A by V, B by V, and C by a, we obtain 

[V, fV, all = V (V, a) - (V, V) a = V (V, a) — Aa (6) 


or 

curl (curl a) = grad div a — Aa, where Aa = 

= A/M + A<?*5 + Afl-k. 

The following table is a pictorial display of second- 
order differential operations: 




138 


Second-Order Differential Operations 


[Ch. V 



Scalar field u 

J Vector field a 

(trad 

d,v 

curl 

peja 


grad div a 










=3 

div grad u = Au 


div curl a = 0 




lino 

curl grad u = 0 


curl curl a 



= grad div a —Aa 




Example 2. The laws o! the classical theorem of electro¬ 
magnetism are postulated in the form of a system of 
Maxwell equation. 

In the most elementary case of a nonconducting homo¬ 
geneous and isotropic medium and in the absence of 
charges and currents, this system is of the form 


fir-iv. hi. 

(7) 

T-f = IV, E). 

(8) 

(V, E) = 0, 

(9) 

(V, H) = 0. 

(10) 


Here, E and H are vectors of the electric-field and mag¬ 
netic-field intensity; e and p are coefficients of the permit¬ 
tivity and permeability (our assumptions are that e 
and p are constants); c is the velocity of light in empty 
space. 

Since the spatial and time derivatives commute, that is, 

i'V. H|_[v.*]. 
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it follows that by differentiating (7) with respect to t, 
we obtain 



Replacing -4^- taken from (8), we get -j- -4^ = 
= -y(V, [V, EU or 

•S-Jr-lv, IV. E]l. (11) 


By virtue of formula (6), [v, IV, Ell = V(V, E) — AE. 
Since (V, E) = 0, it follows from (11) that AE. 

To summarize, for the vector field E we obtain the 
equation 


This is one of the basic equations of mathematical physics 
and is called the wave equation. 

It is easy to see (check this!) that the vector field H 
satisfies the same kind of wave equation 


Thus, under our conditions, each of the coordinates 
E s , E y , E z and H x , H u , H t of the vectors E and H 
satisfies the equation 


d*u 
dt * 


<z* 


! d*u d*u , SPu \ 

\1)x* ~T~ H* ) • 


Here, a = dYTyi is the velocity of propagation of the 
process. In a vacuum, where e = p = 1, we have a = c, 
that is to say, in a vacuum, electromagnetic processes are 
propagated with the velocity of light. 

243. Show that any solution of the equation (v, (V, All — 
— = 0 that satisfies the solenoidal condition satisfies 

the Helmholtz vector equation 


V S A + *»A = 0. 
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Definition. The scalar field u = u (x, y, z) that satis¬ 
fies the condition Au = 0 is termed a Laplace (or harmon¬ 
ic) field. 

Example 3. An important instance of a h armonic field 
is the scalar field u=k/r, fc=constant, r = V x 2 + y 2 + z 2 . 
This function is the potential of a gravitational 
field generated by a point mass placed at the coordinate 
origin. It is'easy to verify that the function u = k/r is 
harmonic everywhere except at the coordinate origin, 
where it is not defined. Indeed, 

(v.v|) = *(v,vi)=A(v, -TfVr) 

= —k (V, -J-r°).« -k (v r°)-*^-(V, r°) 

= -*(-^Vr, r o) 

= %r( r # , r °) —"r(V. r«) = ^--= 0 
for all r =/= 0 since 

(V, r°)=(v, -f) = ( V T‘ r)+-J-(V, r) 

= (-^.r) + A = __J_(rO, r)+ ± 

= — ~(r°, r®) + -5-=-^ + 4 = 4. 

Example 4. Prove that in the potential field of a vec¬ 
tor a its potential function u (x,^y, z) satisfies the Pois¬ 
son equation 

. 5*u . 3*u . d*u , , ,. o, 

Au =ar*+*r»+a? = p(*’ y ' z )• < 12 > 

where p (x, y, z) is the divergence of the vector a. 
Solution. It is given that 

div a = p. (13) 

Since the field of the vector a is a potential field, 
it follows that a = grad u , where u is the potential of 
the field. Substituting a = grad u into 13), we obtain 
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tliv grad u = p or, since div grad u = A u, we have 
Au = p. 

In the special case of points of the field where the diver¬ 
gence is zero, equation (12) turns into the Laplace equa¬ 
tion Au = 0. The Laplace-Poisson equation permits find¬ 
ing the potential function u by integrating the partial 



differential equation. This turns out to be more convenient 
in some problems. 

In electrostatics, preference is often given to finding 
the function v = —u instead of u. Then a = —grad v . 
Accordingly, in the theory of the electrostatic field, the 
Poisson equation is of the form 


d*t> . . <Pv_ 4np 

dx* ’’ dy* ' dz % e 


(14) 


Let us consider an elementary case where the Poisson 
equation is employed. 

Example 5. Suppose two infinite parallel plates AA X 
and BB 1 with opposite charges have potentials v 1 and v a ; 
for the sake of definiteness, v 1 > v t . Find the field E 
between them (Fig. 36). 

Solution. Send the z-axis at right angles to the plates 
in the direction of decreasing potential, and bring the 
yz-plane to coincidence with the positively charged plate 
AA V We now seek the potential function from the Pois¬ 
son equation. By virtue of the symmetry of the problem 
about the z-axis and due to the infinity of the plates, 
we can conclude that the equipotential surfaces are planes 
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parallel to the plates, and the function v depends solely 
on the variable x. Equation (14) takes the form 

< 15 > 

since space charges are absent throughout the space be¬ 
tween the plates. Integrating (15), we find 

v = C x x + C 2 , (16) 

where C x and C 2 are arbitrary constants. 

We require that for x = 0 the function v take the value 
Uj and for x = d, where d is the distance between the 
plates, that it take the value i> 2 . This yields C t =v lt v 2 = 
= C x d + C 2 , whence C 2 = v lt C x = (i > 2 — v x )ld. Substi¬ 
tuting these values of C x and C 2 into (16), we obtain 

v = v ^x + v x =v x - v -±=px. 


The vector E is found from the formula E = —grad v, 
which yields 

I 

so that the field is homogeneous and in the direction of 
the z-axis. The magnitude of E is equal at every point to 
I E | = (v x — v t )/d, that is, it is equal to the potential 
drop per unit length of the shortest distance between the 
plates. 

244. Suppose a scalar function 9 (M) satisfies the Laplace 
equation. Show that the vector V 9 is solenoidal and irrota- 
tional. 

245. Show that A (uv) = u Av + v Au -f 2 (v«, Vi>). 

246. Prove that if r is the radius vector, then 


A r = 


I 


2 . 

— m space, 

— in the plane. 


247. Check to see whether the following scalar fields 
are harmonic or not: 


(a) u = x 2 + 2xy — y 2 , 

(b) u = x 2 y + y 2 z + z 2 x, 

(c) u = x 2 — y 2 . 
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248. Show that the scalar field 

u = ln-^-, where r= j/x 2 + ^ 2 (r^O), 
is harmonic. 

249. Find all harmonic fields that depend solely on x. 

250. Find the general form of a homogeneous harmonic 
polynomial of second degree in x and y. 

251. Find all solutions of the Poisson equation A u — 
= x"~ 2 that depend solely on x. 

Example 6 . Green's formulas. Let q>, ij) he two scalar 
functions of a point. Set up the vector a = <p grad if. 
Then 

div a = div (<p grad tp) = q> div grad ip + (grad 9 . grad if) 
= qp • Aip (grad <p, grad if). 
Now apply the Gauss-Ostrogradsky formula 

^ (a, n°) da = j j j div a dv. 

Note that in our case 

(a, n°) = (<p grad if, n°) = 9 (grad if, n°) = (p-^-. 

We thus obtain Green's first formula: 

j j j [<p Aif + (grad cp, grad if)l dv = 9 da, (17) 

v z 

which for <p = if turns into 

j j j (tp A(p-f|grad(p| 2 ld[; = § (p|^do. (18) 

If we put <p = 1 in (17) we get 

In (17), interchange <p and if and subtract the resulting 
formula 

j j j [if Aq> + (gradip, grad <p)) dv = ^ if da 
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from (17). This then yields Green's second formula: 

J J J (<P-A*-<1>-Aq>)di> = § (<Pl£ —'I’-fr) da * 

V £ 

It is assumed here that all functions that we have to deal 
with and also all their derivatives that occur in the for¬ 
mulas are continuous in the region under consideration. 
Example 7. Find the surface integral 


£ 

taken over the sphere 2: x 2 + j/ 2 + z 2 = 1 for <p = 
= x 2 + y 2 and tp = p 2 -f z 2 . 

Solution. By Green’s first formula, the desired integral 
is 

/ = j j j IqpAip -f (grad cp, grad \|>)) dv. 


where the region of integration V is a sphere: x 2 + p 2 + 
+ z 2 ^ 1. 

We have Ai|> = 4, grad cp = 2xi + 2pj, grad \p = 

= 2pj + 2zk, (grad cp, grad ip) = 4p 2 and therefore 


I = j j j (4x 2 + 4c/ 2 + 4p 2 )du = 4 J j j (x 2 + 2 y 2 )dv. 


Passing to the spherical coordinates x = r cos cp sin 0, 
p = r sin cp sin 0, z = r cos 0, we obtain 

7 = 4 ^ j j (r 2 cos 2 cpsin 2 0 + 2r 2 sin 2 cpsin 2 0) r 2 sin0drd0dcp 

2n n I 

= 4 j (cos 2 cp + 2sin 2 cp) dcp j sin 3 0d0j r k dr 
o oo 

2 n n 

= j (1 -(- sin 2 cp) dcp j (1 — cos 2 0) sin 0 d0 
0 0 

= n ( — cos 0 +cos 3 0) = n. 
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Example 8 . Find the surface integral 

i 

taken over the surface 2: x 2 + y 2 = -ft 2 , z = 0, z = H, 
H > 0, provided <p = i 2 +i / 2 + a: + z, \J) = i 2 + p 1! + 
+ 2 z + x. 

Solution. By Green’s second formula, the desired inte¬ 
gral is 


/ = 




(<p At|> — 1 |) A<p) dv. 


For the given functions 9 and ip we have Atp = 4, Arp = 4 
and therefore 

/=— 4 j | j z du. 

Passing to the cylindrical coordinates x = p cos tp, y = 
= p sin <p, z = z, we obtain 

2n R H 

1 = — 4 j j j zp dpd(p dz = —4 j dtp j pdp j zdz 


Example 9. Find the surface integral 


' dn 


do 


over a closed surface 2 bounded by the planes x + y + z= 
= 1 , x = 0, y = 0 , z = 0 , provided 9 = e* sin y + 1 . 

Solution. The given function is harmonic since Acp = 
= e x sin y — e* sin y = 0. Therefore, by (18) we have 


■w 


I grad 9 p 


We find the modulus of the gradient of the function c 
grad 9 = e* sin y*i + e* cos p*j, | grad 9 | = e x . 
The desired integral is equal to 
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252. Compute the surface integral / = 


<w> , 

( P-g^do over 


£ 

a closed surface 2: (x 2 + y 2 +z 2 = 1, y = 0, y>0}, 
provided 9 = z 2 , 9 = x 2 4-y 2 — z 2 . 


253. Compute the surface integral/= (<p^— ip ) rfa 

£ 

taken over the entire surface of a closed cylinder 2: 
(x 2 + y 2 = l, z = 0, z = 1}, provided <p = 2x 2 , 9 = x 2 + z 2 . 


254. Compute the surface integral ^ = ^ P ro ~ 

£ 

vided (f=(x + y + z)/y 3 and 2 is a sphere: x 2 +y 2 +z 2 =7? 2 . 


255. Find the surface integral I = do, provided 
£ 

(p^^siny + e^sinx+z and 2 is a triaxial ellipsoid 
(x 2 /a 2 ) + (y 2 /6 2 ) + (z 2 /c 2 ) = l. 


Sec. 22. Vector potential 
Let a vector field 

a = a (M) = P (x, y, z) i + Q (*, y, z) j + R (x, y, z) k 
be solenoidal in the region G, that is, div a ( M ) = 0 in G. 

Definition. The vector potential of a vector field a = 
= a (A/) is a vector b (M)=P 1 (x, y, z) i + Q x (x, y, z) j -f 
+ R x (x, y, z) k that satisfies in G the condition 


curl b (M) = a {M) 
or, in coordinate form, 

Mi._i2i.-p Mi._Mi._rt Ml 

dy dz r ' dz dx v ' dx 


(1) 


= R. (2) 


For a solenoidal vector field a (M) the vector potential 
b (M) is not defined uniquely: the condition (1) is also 
satisfied by the vector B ( M) = b (Af) + grad / (M), 
where / (Af) is an arbitrary differentiable scalar function, 
since curl (grad f (M)) = 0. 

Thus, two vector potentials of the solenoidal field a (M) 
differ by the gradient of the scalar field. 
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Finding the vector potential b ( M) of the solenoidal 
field a ( M) reduces to finding some particular solution of 
system ( 2 ) of three partial differential equations for the 
three unknown functions P 1 (x , y, z), Q t (x, y, z), 
Ri ( x , y, z). 

The vector potential b ( M) may be constructed in the 
following manner. Taking advantage of the arbitrariness 
of choice of the vector b ( M ), we will simplify matters by 
setting P x Or, y, z) = 0, that is, the vector b ( M) will 
be sought in the form b ( M) = Q x (x, y, z) j + 
+ R x (x, y, z) k. Then the system of differential equa¬ 
tions ( 2 ) for finding the unknown functions Q x (x, y, z) 
and R x (x, y, z) takes the form 


BRy_ 

dV 


3Qx P 

dz 



dQi 

dx 


R. 


(3) 


From the second and third equations of this system we 
find 

Ri(x, y, z) = - j Q(x, y, z)dx + C x ( y , z), 

Qi(x, y, *)= j R(x , y, z)dx + C 2 {y, z), 

where C 1 (y, z) and C a {y, z) are any differentiable func¬ 
tions of y and z. For the sake of simplicity, set C 2 ( y , z) = 
= 0 and choose the function C x (y , z) so that the first 
equation of system (3) is satisfied. To do this, we sub¬ 
stitute into the first equation the expressions that were 
found for (?, and R x : 

?(*. y, z ) dx +^f--h\ *(*• y> *)<** 

= P (x, y, z). 

From this we get 


liy~ = 'k I Q(x ' y ' z)dx + ~k\ fl(x ’ y ' z ) dx + p (^V’ z )- 


10* 
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It is easy to verify that the right-hand side of this 
equation does not depend on x; this is because div a ( M) = 
= 0 in G. 

Integrating the last equation with respect to y, we find 

Ci ( y, z) = j [-^ j Q (x, y, z) dx 

+ 17 I l,*)dx+P(x t y, z)] dy 4 C 3 (a). (4) 

Setting C , (z) = 0 in (4) and substituting (4) into the 
expression for R x (x, y, z ), we get a particular solution 
of system (3): 

*1-0, (5) 

Qi = j R (x, y, z) dx, ( 6 ) 

*i = J ["Si I y ' H ( x ’ y ' z ) dx 

+ P ( x, y, z)] dy - j Q (x, y, z) dx. (7) 

The vector b ( M ), whose coordinates P x (x, y, z), 
Qi (*. y, z), R x (x, y, z) are defined by formulas (5), 
( 6 ), (7), is the vector potential since it satisfies the condi¬ 
tion curl b = a. — 

Example 1. Find the vector potential b = b (x, y, z) 
for a solenoidal field given by the vector 
a = 2yi — zj + 2xk. 

Solution. We seek the potential b in the form 
b = b (x, y, z) = <?, (x, y, z) j + R x {x, y, z) k, 

where Q x (x, y, z) and R x (x, y, z) are found from (7) 
and (8). Since in the given case P = 2y, Q — —z, R = 
= 2 x, we obtain 

Q x (x, y, z) = j 2xdx = x 2 , 

R l (x, y, z)= j zdx+ j 2y dy = xz + y 2 . 

Thus 

b (x, y, z) = x 1 j + (xz + y 2 ) k. 
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It is clear, by direct verification, that curl b = a, 
and, hence, this vector is the vector potential of the given 
field. 

Remark. Due to the arbitrariness in choosing the vector 
b, we could require Q x (x , y, z) = 0 or f? x (x, y, z) = 0 
instead of P x (x, y, z) == 0. The system of equations (2) 
and formulas (5), ( 6 ), (7) would naturally change. 

Find the vector potentials of the following solenoidal 
fields: 

256. a = i + j + k. 

257. a = 2yi+ 2zj. 

258. a = (e* - e") k. 

259. a = 6yH + 6zj + 6xk. 

260. a = 3yH - 3x*j - (y 2 + 2x) k. 

261. a = ye** i -f 2yzj — (2xyze** -f z s ) k. 

If the vector field a = a (M) is solenoidal in the re¬ 
gion G y which is star-shaped (see Sec. 19, Chapter IV) with 
centre at the coordinate origin O ( 0 , 0 , 0 ) [the field a ( M) 
may not be defined at the point O ], then one of the vector 
potentials b = b ( M) may be found from the formula 


b(A/)= j [a (A/'). r(M)]tdt, (8) 


where r ( M) = xi + y\ + zk is the radius vector of the 
point M (x, y,' z) and the point M' (tx, ty, tz) runs 
over the line segment OM as the parameter t varies 
from 0 to 1 . 

Example 2. Use formula ( 8 ) to find the vector poten¬ 
tial of the solenoidal field 

a = 2 yl - zj + 2 xk. 

Solution. The given vector field is defined throughout 
three-dimensional space, which is a star-shaped region 
with centre at the coordinate origin, and so we can use ( 8 ) 
to find the vector potential. At the point M' (tx, ty, tz) 
we have 


a (M') = 2ty\ — tz j 4- 2lxk. 
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We find the vector product 

i ) k 

[a (M'), r (M)]= 2ty -tz 2 tx 
x y z 

= [ - (2 xy + z 2 ) I + (2x 2 - 2 yz) j + (2 y z + xz)k) t. 
Using (8), we obtain 

t 

b (AT) = j [ — (2xy -f z 2 ) i + (2x 2 — 2yz) j + (2i/ 2 + xz) k\t 2 dt 
o 

= - -i- (2xy + z 2 ) i + (x 2 - yz) j +• -1 (2|/ 2 + xz) k. 

It is easy to establish that curl b ( M) = a ( M ). 

Remark. In examples 1 and 2 we obtain different vector 
potentials for one and the same solenoidal field a = 
= 2yi - zj + 2xk: 
b, (M) = x z i + (xz + y z )k, 

*h (M) = - ± (2 xy + z 2 ) i + (x 2 - yz) j f ± (2y 2 + xz) k. 

They differ by a term equal to the gradient of some 
scalar field / ( M ). This term plays the role of an arbitrary 
constant (when acted on by the curl). It may be repre¬ 
sented as the gradient of some scalar function, / (M). 
Let us find this function in our example. We have 
grad /(Af) = b, (MJ-b^M) 

= (2 xy + z 2 ) i + i- (x 2 + 2 yz) j + -j- (2xz + j/ 2 ) k. 

To find the scalar field / (M), use formula (3) of Sec. 19, 
in which we take the coordinate origin 0 (0, 0, 0) for 
the point (x 0 , y 0 , z 0 ). This yields - 

xy z 

f(M)= j 0-dx + j Y 3fidx + J -Y@ xz + y z ) dz + C 
oo o 

=- 3 - ( x 2 y + y * z + z2x )+ c • 


where C is an arbitrary constant. 
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Example 3. Find the vector potential b of a magnetic 
field H set up by an electric charge e that is moving with 
a constant velocity v. 

Solution. By the Biot-Savart law, the intensity of the 
magnetic field is 

H W = -fe L (9) 


where r is the distance of point M from the charge e. 

Since H is a solenoidal vector, that is, div H = 0, 
there exists for it a vector potential b such that H = 
= curl b or, taking into account formula (9), 


curlb 


ley, r) _* [v, r] 

4nr* 4 ji r* ’ 


Let us rewrite this formula as 
c<i rl b = 3 j{[«y, -i <] + [«*, 7T l] + [«. -?' 1 ]} 

-£{[«• -t]+[i. -*•]+[* -S 1 ]} 
-M*. *(?)] + [!.*(?)] + [*.*(-?)]}■ 

Employing the readily verifiable equation 

we obtain 


whence 


curl b = -r—curl - 



Using formula (8), find the vector potentials of the 
following solenoidal fields defined in star-shaped regions: 

262. a = I. 

263. a = 6zi — 15yj + 9zk. 

264. a = 5xV — lOxj/zk. 

265. a = 2 cos xz- j. 

266. a = , x 2 +y 2 > 0. 
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Sec. 23. Curvilinear coordinates 


In many problems it is more convenient to define the 
position of a point M in space by three numbers (ft, ft, q 3 ) 
instead of the three Cartesian coordinates ( x , y, z). These 
numbers often prove to be more suitable to the problem at 
hand. 

Let every point M be associated with a definite number 
triple (ft, ft, ft) and, conversely, let every number triple 
correspond to a unique point M. Then the quantities 
ft, q 2 , q 3 are termed the curvilinear coordinates of the 
point M. 

The coordinate surfaces in a system of curvilinear 
coordinates ft, ft, ft are the surfaces 



(1) 

( 2 ) 

( 3 ) 


on which one of the coordinates remains constant. 

The line of intersection of two coordinate surfaces is 
called a coordinate line (axis). 

The coordinates ft and q 3 maintain constant values 
along the line of intersection of the coordinate surfaces 
(2) and (3); it is only the coordinate q 1 that varies. Simi¬ 
larly, on the lines of intersection of the surfaces (i) and (3) 
and (1) and (2), it is q 2 and q 3 that vary respectively. 

Let us introduce the unit vectors e a , e 2 , e 3 directed 
along tangents to the coordinate axes, (ft), (ft), (ft) at 
the point M in the direction of increasing variables g lt ft, 
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q a respectively (Fig. 37). Let us agree to take the unit 
vectors e s , e 3 always in that order so that, taken 
together, they constitute a right-handed trihedral. 

The basic difference between curvilinear coordinates and 
Cartesian coordinates is this. In the Cartesian system, the 
vectors e 2 , e 3 are constant for all points of space 
and are equal, respectively, to i, j, k. In any other 



system, they will, generally speaking, change their direc¬ 
tions when passipg from one point, M, to other points. 

As examples of curvilinear coordinates we consider 
cylindrical and spherical coordinates. 

1°. Cylindrical coordinates. The position of a point Af 
in space is defined in cylindrical coordinates by three 
coordinates: 

< 7 i = p, 0 < p< +oo, 

q-t = <p, 0 sg; <p < 2 ji, (4) 

q 3 = z, — oo<z<+oo. 

The coordinate surfaces are: 
p = constant: circular cylinders with the z-axis; 

<p = constant: half-planes adjoining the z-axis; 
z = constant: planes perpendicular to the z-axis. 

The coordinate lines (or axes) are: 

(p): rays perpendicular to the z-axis and having their 
origin on that axis; 
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(q>): circles with centre on the z-axis and lying in 
planes perpendicular to that axis; 

(z): straight lines parallel to the z-axis (Fig. 38). 

Cartesian coordinates are 
related to cylindrical coor¬ 
dinates by the following 
formulas: 

X = p cos <p, 
y = p sin tp, (5) 
z — i. 

2°. Spherical coordinates. 
In spherical coordinates, 
the position of a point M 
in space is defined by the 
following coordinates: 
ft = r, 0 < r < 4-00, 
q 2 = 9, 0 < 0 < ji, (6) 

q 3 = 9, 0 ^ <p < 2n. 

The coordinate surfaces 
are (Fig. 39): 

38 r = constant: spheres 

centred at 0: 

0 = constant: circular half-angle cones with the 
z-axis; 

cp = constant: half-planes adjoining the z-axis. 

The coordinate lines are: 

(r): rays emanating from the point 0; 

(0): meridians on a sphere; 

(<p): parallels on a sphere. 

Cartesian coordinates are related to spherical coordi¬ 
nates via the following formulas: 

x = r cos <p sin 0, 

y = r sin <p sin 0, (7) 

z = r cos 0. 

A system of curvilinear coordinates is said to be or- 
thogonal if at every point M the unit vectors e lt e 2 , e 3 
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are pairwise orthogonal. In such a system, the coordinate 
lines and the coordinate surfaces are also orthogonal. 
Systems of cylindrical and spherical coordinates are instan¬ 
ces of orthogonal curvilinear systems of coordinates. 



Henceforth we consider only orthogonal systems of coor¬ 
dinates. 

Suppose r = r (q lt q 3 , q 3 ) is the radius vector of 
a point M. Then 

dr = H x dq x ^ + H 2 dq 3 e t + H 3 dq 3 e 3 . (8) 

Here 

'.-/&)*+(*)*+(*r) 2 . 

are the Lam6 coefficients of the given curvilinear system 
of coordinates. 

In cylindrical coordinates, 
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By virtue of (5) we have 

) ! +(ir+(f) ! ='- 

In spherical coordinates, 

9i = r > 9s = ®» 9s = *P- 
By virtue of (7) we have 

«.-"-=/(£) 2 +(£)‘+(£r-!• 

(t) 2 +(!) 2 +M) 2 ='- 

»-*•-/ (^) 2 +(^) 2 +(^r=”-e- 

The quantities 

dZ, = tf, dg„ i = 1, 2, 3 

that appear in formula (8) are differentials of the arc 
lengths of the coordinates lines. In a number of cases, this 
consideration permits of a more simple computation of 
the Lame coefficients. For instance, in the case of cylin¬ 
drical coordinates (4) (see Fig. 38), the differentials of 
the arc lengths of the coordinates lines (p), (q>). (z) are 
d (p) = 1-dp, whence H x = 1; 
d (<p) = p-dcp, whence H 2 = p; 
d (z) = 1-dz, whence H 3 = 1. 

It is just as easy to obtain expressions for the Lame 
coefficients in the case of spherical coordinates (6). 


Sec. 24. Basic operations 
of vector analysis in curvilinear 
coordinates 

1°. Differential equations of vector lines. Suppose we 
have a vector field 

a = <h (? 1 . ? 2 . 7s) e i + a 2 (q x , q 3 , q 3 ) e 2 + a 3 (q u q 3 , q 3 ) e 3 . 
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The vector-line equations in curvilinear coordinates 
<7i, 7 2 , 7 s are of the form 

//, d Qi _ H t dq t _ H 3 dq 3 
“l (9i. 7s. 7 s) <*s (<? 1 . 9s. 7s) “ “s (?i. 7s. 9s) ' 

In particular, in cylindrical coordinates (ft = p, ft = <p, 
7s = z): 

_ P ^ ^ 

MP. <P. *) «s(P. 1 >. 2) fl s (P. <P. 2) ’ ' ' 

in spherical coordinates ( q x = r, q 2 = 0, q 3 = <p): 

dr r dd _ rsin 6d<p 

<*i (r, 0. q>) — a, (r, 0, q>) — a, (r, 6, <p) * 

Example 1. A vector field is given in cylindrical coor¬ 
dinates 

a (M) =e p + <pe„. 

Find the vector lines of the field. 

Solution. It is given that a t = 1, a 2 = <p, a 3 = 0. By 
virtue of formula (1) we have 

dp _ p dtp _ dz 
1 <p — 0’ 

whence 

z = C x , 

P = C*<p, 

which are Archimedean spirals lying in planes parallel 
to the xy-plane. 

2°. The gradient in orthogonal coordinates. Suppose we 
have a scalar field 


u = u (ft, q t , 7 j). 


grad u = -jj- e, + e 2 4- e 3 . 

In particular, in cylindrical coordinates (q 2 = p, q 2 
= cp. 7a = z): 

1 


1 du 


1 du 


grad u = 


- e <» + 


dz 


( 2 ) 


in spherical coordinates (q 1 = r, q 2 = 0, 7 S = cp 

, du , 1 du , 1 du 

et adu = — e r + - 15 e 0 +- 


( 3 ) 
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Example 2 . Compute the gradient of the following sca¬ 
lar field specified in cylindrical coordinates (p, cp, z): 
u = p 4- z cos 9. 

Solution. Using formula ( 2 ), we obtain 

grad u= 1 e p -i-zsin9e (p + cos9e I . 

Example 3 . Find the gradient of the following scalar 
field given in spherical coordinates (r, 0, 9): 

, sinB . a 
u — r- j-sin 6 cos 9. 

Solution. Using formula ( 3 ), we have 

. t. sin0\ , cos8 /I \ sin® 

grad u= (1--— J e r -l--— —cos 9 J ee-|-^-e,. 

3 °. The curl in orthogonal coordinates. Suppose 
a=fli (?i. 9s) + a 2 (q x , q 2 , q 3 ) e 2 + a, ( 9l , q 2 , q 9 ) e 3 . 

Then 

1Tp{ 7 e ‘ TFJU** lTJh e3 
curl a = _J_ 

<ty» ^9s 

a \H 1 a zH 2 

In particular, in cylindrical coordinates (q x = p, q 2 = 
= 9, g 3 = z): 



curl a = (4) 

dp dip dz 

a, pa 2 a 3 

in spherical coordinates (q x = r, q 2 = 0, g 3 = 9): 

1 c 1 c 1 . 

r 3 sin 0 Cr rsin 8 ee r e * 

curl a = (4') 

dr d0 dip 

a, raj rsin0*a 3 
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Example 4. Compute the curl of the following vector 
field specified in cylindrical coordinates: 

a = sin<p-e p + e y —p ze t . 

Solution. Taking advantage of formula (4), we get 


r e P «• 


= 4- e p (0 — 0) — (— z — 0) 


dp d(f dz 

sin<p cos<p — pz| 

. 1 /n ^ cos ffi 

-|—— C z (0 — COS <p) = Zttif -j^-e*. 


4°. The divergence in orthogonal coordinates. Suppose we 
have a vector field 

a=«i (9i. 9»» 9s) Ci + a t (g„ q 2l q 3 ) e s + a, (g„ g 3 , g 3 ) e 3 . 
Then 

diva = _i_r d (a t ff t // 3 ) a(* t fW *(«,1W 1 

//j/f t // 3 L dgj 09* 09* J * 

In particular, in cylindrical coordinates (q 2 = p, 7* = <p, 
9a = z): 

diva _1 a (pa,) 1 1 da t da t . 

a a — p ap + p aq> + a* ’ 

in spherical coordinates (g 3 = r, g 3 = 0, g 8 = <p): 


diva 


1 ^1) 


a (sin 9-a*) 


r sin 0 90 "^rsin 9 0<p * 

Example 5. Show that the vector field 


(5) 


is solenoidal. 

Solution. Using formula (5), we have 

=-M-^)+ 2 sin 6 cos 6 = 0 

wherever r # 0, which means that the field of the vector 
a is solenoidal at all points with the exception of r = 0. 
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267. Find the equations of the following vector fields: 

(a) a = e p + ye, + e,; 

(b) a = pe p + <pe* + ze I ; 

. . 2a cos 0 , asinfl . 

(c) a =-5—e r H-3—e e , a — constant. 

Find the gradients of the following scalar fields. 

(a) In cylindrical coordinates: 

268. u = p 2 + 2p cos <p — e* sin <p. 

269. u = p cos <p + z sin 2 <p — 3 p . 

(b) In spherical coordinates: 

270. u = r 2 cos 0. 

271. u = 3r* sin 0 + e T cos <p — r. 

272. u = p -521®. t p, = constant. 

Compute the divergence of the following vectors. 

(a) In cylindrical coordinates: 

273. a = pe p + z sin ip-e, + e* cos z-e z . 

274. a = <p arctan p-e„ + 2e, — zVe r . 

(b) In spherical coordinates: 

275. 8 = ^ —2cos 2 (p ee + -jjq-j-e v . 

Compute the curl of the following vector fields: 

276. a = (2r + a cos <p) e r — a sin 0-e o + r cos 0-e.p, 
a = constant. 

277. a = r*e r + 2 cos 0-e e — cpe v . 

278. a = cosy e p — Cy+p^e^. 

279. Show that the vector field 


is a potential field. 

280. Show that the vector field 

a = f(r) e r , 

where / is any differentiable function, is a potential field. 

5°. Computing the flux in curvilinear coordinates. Let 
S be a part of the coordinate surface q l = C, where C = 
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= constant, bounded by the coordinate lines 
9i = o„ q t = a* (a x < a,); 

9s = Pn 9s = Pi (Pi < Ps)- 

Then the flux of the vector 

a = «i (9i. 9n 9s) e i + As (9i. 9i. 9s) e 2 + <*3 (?i. 9i. 9s) e 3 
through the surface S in the direction of the vector e x is 
computed from the formula 

“« b 

n = j j a t (C, q 2 , 9s) H z (C, 92* 9s) ^ 3 {C, q 2 , 9s)^9s ^92- 

( 6 ) 

The calculation is similar for the flux of a vector 
through a part of the surface q t = C or through a part of 
the surface q 3 = C, where C = constant. 

Example 6. Compute the flux of the vector field, speci¬ 
fied in cylindrical coordinates, 

a = pe p + ze v 

through the outer part of the lateral surface of the cylin¬ 
der p = 1 , bounded by the planes z = 0 and z = 1 . 

Solution. The cylinder is the coordinate surface p = 
= C — constant and so the desired flux 

2n 1 

n= j j C 2 dzd<f = 2nC 2 , 
o o 

whence for the surface p = 1 we obtain 

n = 2 n. 

Example 7. Find the flux of the vector field, specified 
in spherical coordinates, 

a = r* 0 e r + re 2 0 e e 

through the outer side of an upper hemisphere S of ra¬ 
dius R with centre at the coordinate origin. 

Solution. The hemisphere S is part of the coordinate 
surface r = constant, namely, r = R. On the surface S 
we have 


q t = r = R\ q z = Q, O<0<-2- ; 

9s = <P, 0<<p<2n. 
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Taking into account that in spherical coordinates 
H 1 = H r = 1, H t = H$ = r, H a = Hy = r sin 0, 
we find, via (6), that 

jt/ 2 2n n/2 

11= | dej /i 4 0 sin Qd<p = 2nJP j 0 sin 0 d0 = 2n7? 4 . 

e o o 

Compute the flux of the vector field, specified in 
cylindrical coordinates, through the given surface S. 

281. a = pe p — cos tp-e,, -f ze 2 ; S is a closed surface 
formed by the cylinder p = 2 and by the planes z = 0 
and z = 2. 

282. a = pe p + pepe* — 2ze t ; S is a closed surface 
formed by the cylinder p = l, the half-planes <p = 0and 
q> = n/2, and by the planes z = — 1 and z = 1. 

283. Find the flux of the vector field a = (1 Ir 2 ) e r through 
a sphere of radius R with centre at the coordinate origin. 

284. Find the flux of the vector, specified in spherical 
coordinates, 

a = re r + r sin 0>ee — 3rq> sin 0-e 9 

through an upper, hemisphere of radius R. 

285. Find the flux of the vector, specified in spherical 
coordinates, 

a = r?e r + R 2 cos q>-e,, 
through the sphere r = R. 

286. Find the flux of the vector, specified in spherical 
coordinates, 

a = re, — r sin 0 -e^, 

through a semicircle of radius R located in the half-plane 
cp = n/4 (the flux is taken in the direction of the vector 
e„). 

287. Find the flux of the vector, specified in spherical 
coordinates, 

a = r sin e e + r sin 0 cos <p .e«, 

through the outer side of part of the half-angle cone 
V 3z 2 = x 2 + y 2 , bounded from above by the plane z = 
= ]/3 (0 < z < |/3). 



Sec. 24) 


Basic operations of vector analysis 


6°. Finding the potential in curvilinear coordinates. Giv¬ 
en, in curvilinear coordinates ft, q 3 , q 3 , the vector field 
a (A/) = a 3 (ft, q t , q 3 ) e x + ft (ft, ft, ft) e 8 + 

+ ft (ft, ft, ft) e 3 ; 


this is a potential field in some region Q over which 
the variables ft, ft, ft range, that is, curl a = 0 in Q. 

To find the potential u = u (ft, g 2 , ? 3 ) of this field, 
write the equation a (M) = grad u ( M) as 

. 1 du , 1 da 1 du 

a ‘ e * + = it : ^ e 2 + 777 e3 - 


From this it follows that 


(7) 


This is a system of partial differential equations whose 
integration yields the desired potential u=u (ft, q t , q 3 ) -f- 
+ C, where C is an arbitrary constant. 

The system (7) of differential equations is solved in 
the same way as in finding the potential in Cartesian 
coordinates 

The system (7) of differential equations is of the fol¬ 
lowing form: 

(1) In cylindrical coordinates (ft = p, ft = <p, ft = 2 ), 


du 

dtp 


p a 9 , 


du 

dt 


a z . 


(7') 


(2) In spherical coordinates (ft = r, ft = 


du 


du 


77 = a - -dQ= ra »' -^ = rsin0 -<V 


= <P). 
(V) 


Example 8. Find the potential of the following vector 
field specified in cylindrical coordinates: 

a = jarcUn i + cog ep _ sin <p-e v + e t . 

Solution. By formula (4) we find 

-S U\ 


curia = 


arctan z , In p 

—--(- cos q> — p sin q> 


= 0 (p>0). 
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This is a potential field. The desired potential u = 
= u (p, (p, z) is the solution of the following system of 
differential equations: 


du 

Ip 


du 

— = — psin qp, 

du _ In p 
~dz~ 1 + z* * 


From the first equation, integrating with respect to p 
we obtain 


u = In p<arctan z + p cos <p + C (q>, z). 
Differentiating (8) with respect to <p, we get 
du . dC 

--- psm v + — 

and since du/dy = —p sin q>, it follows that dC/dy = 
that is, C = C x (z). Thus 

u = In p-arctan z + p cos q> + C x (z), 

whence 


( 8 ) 

0 , 


^ = T^r+c;W. 

By virtue of the third equation of the system we have 

T^-dV+c;<’> 


or C[ (z) = 0, whence C x (z) = C = constant. 

To summarize, the potential of the given field is 
u (p, <p, z) = In p • arctan z -f p cos q> + C. 

In the following problems, verify that the vector fields 
given in cylindrical coordinates are potential fields and 
find their potentials. 

288. a = e P H— 


289. a = pe p + -^-e, + ze 2 . 

290. a = (pzep + ze^ + pqpe*. 

291. a = sin q> -e p + y eP cos <p e, + 2ze 2 . 
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292. a = q>cosz-e p + cosz-e v + p(psinz-e I . 

Example 9. Find the potential of the following vector 
field given in spherical coordinates: 

a = -i- e**e r + (pe^ee. 

Solution. Using (4'), we find that 

e r re 0 rsinO-e, 

1 J_ _d_ _d_ 

curla = 7*slne dr 36 *<*> “°- 

-i- e** (plnr-e 9 * Olnr-e 9 * 

This is a potential field in the region where r > 0,0 

(» - 0 , ±1 _). 

The system (7) of differential equations for finding the 
potential u = u (r, 0 , <p) is of the form 

-|£ = cpe 9 * In r, (9) 

~ = Qe&f in r. 

Integrating the first equation of system (9), we obtain 
u = e*f In r + C (<p, 0). (10) 

Differentiating (10) with respect to 0 and taking into 
account the second equation of the system, we have 

(pe 8 * In r = In r + . 

That is, dC/dQ = 0, whence C (<p, 0) = Cj (<p) and 
therefore 

u = e 9 * In r + C x (<p). (11) 

Differentiating (11) with respect to (p and taking into 
account the third equation of system (9), we obtain 
Be 9 * In r = Be 9 * In r + CJ (<p) 

or C\ (<p) = 0, whence C x (q>) = C = constant. The 
desired potential is 

u (r, 0, q>) = e 9 * In r + C. 



166 


Curvilinear Coordinates 


[Ch. VI 


Establish the potentiality of the following vector 
fields, specified in spherical coordinates, and find their 
potentials. 

293. a = 0e r + e 0 . 

294. .-2re, + 7 J r5 e,+J-e,. 

295. —4-1*,+^ %+■?■«»• 

296. a = cos <p si n 0 • c r + cos <p cos 0 • ee — sin <p • e*. 

297. a = e r sin9-e r +-^-e r cosO-ee+ (1 . Hq ,^ rsin0 e»- 

7°. Computing the line integral and the circulation of 
a vector field in curvilinear coordinates. Suppose a vector 
field 

a (M) = ai (ft, q 2 , q 3 ) e x + a 2 (q v q t , q 3 ) e 2 

+ a 3 (?„ q 2 , q 3 )e 3 

is defined and is continuous in a region Q over which the 
orthogonal curvilinear coordinates q lt g 2 , q 3 range. 

As we know [see Sec. 23, (8)1, the differential dr of 
the radius vector r of any point M ( q x , q 2 , q 3 ) 6 Q is 
equal to 

dr = H l dq 1 e 1 + H 2 dq t e 2 + H 3 dq 3 e 3 . 

Therefore the line integral of the vector a ( M) over an 
oriented smooth or piecewise smooth curve L c Q is 

j (a, dr) = j a l tI i dq l -\-a 2 H 2 dq 2 -\-a 3 H 3 dq 3 . (12) 

In particular, for the cylindrical coordinates q x = p, 
g 2 = <p, < 7 3 = z we have 

a = a P (p, <P, z) e p + a,, (p, <p, z) e„ + a t (p, <p, z) e z , 
dr = dp*e p + pd<p-e, + dz-e z , 
and therefore 

j (a, dr) = j a p dp -f a»p d(p + a z dz. 


(13) 
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For the spherical coordinates q x = r, q 2 = 0, q 3 = <p we 
have 

a = a r (r, 0 , <p) e r + a fl (r, 0 , <p) e 6 + a* (r, 0 , (p) e„, 
dr = dr-e, + r dQ-e e + r sin 0 dq> -e,,, 
and consequently 

| (a, dr)= | a r dr + ra e dQ + raisin 6 dtp. (14) 

The circulation C of the vector field a ( M) is computed 
in the curvilinear coordinates g lt q 2 , q 3 via formula ( 12 ) 
in the general case; for cylindrical or spherical coordi¬ 
nates it is computed from (13) and (14) respectively. 

Example 10. Compute the line integral in the vector 
field, given in cylindrical coordinates, 

a = 4p sin (p-e p + zePe, + (P + <?) e z 
along the straight line 



from the point 0 (0, a/4, 0) to the point A (1, a/4, 0). 
Solution. In the given example, 

a„ = 4p sin <p, a* = ze p , a z = p + <p. 

By formula (13), the desired line integral is 
f (a, dr)= j 4p sin <p dp + pzeP dtp -f (p + 9) dz. 

L L 

On the straight line L we have 

<p = -J-, dq, = 0; z = 0, dz = 0; 0<p<1. 
Therefore 

1 

J (a, dr)= j2l/2pdp=l/2 J 2 pdp = l/2. 

L L 0 

Example 11. Compute the line integral in the vector 
field, given in spherical coordinates, 

a = e r sin 0-e r 4- 30 2 sin (p-ee + np0e w 
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along the line 

f r = 1 ’ 

in the direction from the point Af 0 (1, 0, nl 2) to the 
point M x (1, ji/2, ji/ 2) (Fig. 40). 



Fig. 40 


Solution. The curve L is an arc of a circle with centre 
at the coordinate origin and radius r = 1 located in the 
yz-plane. The coordinates of this vector are 

a r = e T sin 9, a e = 30* sin q>, = r<p0. 

By virtue of (14) the line integral is of the form 
j (a, dr)= j e r sin0dr + 30 2 rsin(pd0 + r 2 <p0sin0d(p. 

Taking into account that the following conditions hold 
on L, 

r = l, dr = 0; <p = -£-, d<p = 0; O<0<-£-, 
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we obtain 

n/2 

j (a, dr) = J 302d0 = j 302d0 = -£. 

Example 12. Compute the circulation of the vector 
field, given in cylindrical coordinates, 

a = p sin <p-e p 4- p ze, + p 3 e* 
over the curve 


{ p = sin<p, 
2 = 0 , 




directly and via the Stokes theorem. 

Solution. The coordinates of the vector are 
a p = p sin <p, ay = pz, a t = p 3 . 


The contour L is a closed curve located in the plane z 
(Fig. 41). 

(1) Direct calculation of the circulation. 


0 



Substituting the coordinates of the vector into (13), 
we obtain 

C = p sin <p dp + P 2 z d(f + p 3 dz. 

L 

On the curve L we have 

z = 0, dz = 0; p = sin cp, dp = cos <p d<p, 

0 < <p < ji- 
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Therefore the desired circulation is equal to 

C = ^ p sin <p dp = j sin 2 <p cos cp dip = 0. 

L 0 

(2) Computing the circulation via the Stokes theorem. 
By the Stokes theorem, the desired circulation is equal 


C = (a, dr) = j j (curl a, n°) dS, 

L S 

where S is the surface spanning the contour L. 
We find the curl of the given field: 


e P 

pe* 


d 

a 

a 

dp 

dy 

dz 

psin<p p 2 z 

P 3 


+ ( 2 z—cosq>) e t . 

At points where p = 0 we redefine the value of curl a 
with respect to continuity, setting 


curl a (0, <p, z) = (2z — cos q>) e z . 

Thus, curl a is defined throughout three-dimensional 
space. Since the curve L lies in the plane z = 0, for 
the surface S spanning this curve we take that portion 
of the plane z = 0 that is bounded by the curve L. Then 
we can take the unit vector e z for the unit vector of the 
normal n° to the surface S, that is, n° = e z . We find 
the scalar product: 

(curl a, n°) = (—pe p — 3p 2 e^ + (2z — cos q>) e 2 , e*) 

= 2z — cos cp 


because by virtue of the orthonormality of the basis 
ep> e«i, e* we have 


(e P , e z ) = (e„ e z ) = 0 , (e lf e 2 ) = 1 . 
The desired circulation is 
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Taking into account that z = 0 on S and the element of 
area dS of the coordinate surface z = 0 is equal to 
dS = p dp d<p, 

we finally get 

C = — j j cos <p dS = — j j cos <pp dp d<p 


= — C cos qp d<p j pdp =—y j sin 2 q> cos q) dq> = 0. 
oo o 


Example 13. Compute the circulation of the vector, 
specified in spherical coordinates, 


a = re r + (R + r) sin 0-e„ 
around the circle 



in the direction of increasing values of the angle q>, 
directly and via the Stokes theorem. 

Solution. In this example, 

a T = r, a 0 = 0, a* = (R + r) sin 0. 

(1) Direct computation of the circulation. 

By formula (14) the desired circulation is equal to 


C = r dr (R + r) sin 0 r sin 0 d(p = 


= rdr + r (R + r) sin 2 0dq>. 


On the given circle L, the centre of which lies at the 
coordinate origin, we have 

r=R, dr = 0; 0 = -J; 0<cp<2n, 
and, consequently, 

2 n 

C = 2B 2 §dq> = 2R 2 j dq> = 4 nR 2 . 

L 0 
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(2) Computation of the circulation via the Stokes 
theorem. 

The desired circulation is, by the Stokes theorem, equal 
to 

C= ^ (a, dt)= j j (curia, n °)dS, 

where S is the surface spanning the circle L. 

We find the curl of the given vector: 

e r ree rsin 0 -e, 

l 1 9 d a 

curia - rlsin0 dr m dq) 

r 0 (ffr + r 2 ) sin 2 6 

= -j (if + r) cos 0 e r —(if + 2 r) sin 0 -e e . 

For the surface S spanning the circle L we take, for 
example, the upper hemisphere of radius if: r = if, 
0^0^ ji/ 2, 0<q) < 2 ji. The unit vector of the normal 
n° to the outer side of the hemisphere S is directed along 
the vector e r , and so we take n° =e r . We find the scalar 
product 

(curl a, n°) = cos 0 -e r — R + 2r s i n 0 . Cei Cr J 

= cos 0 

since (e r , e r ) = 1 , (e r , e e ) = 0 . 

Taking into account that r = if on the surface S, we 
obtain the following expression for the desired flux: 

n = ( j 2( - * + r) cos 0 dS = 4 j j cos0dS. 

In spherical coordinates, the element of area dS of the 
coordinate surface r = if, that is, the hemisphere S, is 
equal to 


dS = if 2 sin 0 dQ dq> 
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and, consequently, 
n = 4j j cos0 R 2 sin0 dQ dtp 

n/2 2n 

= 4/? 2 j cos0sin0d0j d(p = 4nR z . 
o o 

By the Stokes theorem we get C = 4 nR 2 . 

For the surface S spanning the circle L we could take 
the lower hemisphere, the unit vector of the normal of 
which, n° = — e M and the result would be the same: 
C = 4ji/f*. 

Note that it is not desirable to take the circle bounded 
by L for the surface S spanning the circle L because the 
circle has a point r = 0 (the centre of the circle) in which 
the curl of the given vector has a discontinuity. 

Compute the line integral along the given curves L in 
the following vector fields specified in cylindrical coordi¬ 
nates. 

298. a = ze p + p<pe„ + cos cp e r ; L is a segment of the 
straight line: {p = a, 9 = 0 , 0 ^ z ^ 1 }. 

299. a = pe p + 2pcpe<p + ze z ; L is the semicircle: {p = 1, 
2 = 0 , 0 < 9 < n). 

300. a = e p cos <p-e p + p sin q)-e v + pe z ; L is a turn 

of the helical curve: {p = R, z = <p, 0 9 ^ 2ji). 

Compute the line integral over the given curve L 
in the following vector fields given in spherical coordi¬ 
nates. 

301. a = e T cos0-e r + 20 cos cp-e e + <pe v ; L is the semi¬ 
circle: {r = 1 , 9 = 0 , 0 ^ 0 ^ ji}. 

3u2. a = 4r 3 tan -y e r -f 09 ee-f-cos 2 9 -e„; Lis a segment 
of the straight line: { ( P = _ jr; 0 = -j-, O^rs^lj. 

303. a = sin 2 0-e r + sin O-ee + ^Oe*; L is a segment of 
the straight line: {q> = -y» r = lf^e • 

Compute the circulation of the following vector fields, 
specified in cylindrical coordinates, about the given con- 
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tours, directly and with the aid of the Stokes theorem. 

304. a = ze p + p ze v -f pe 2 ; L is the circle: {p = 1, 
z = 0 }. 

305. a = p sin <p-e p — p 2 ze„ + p 1 ^; L is the circle: 
{p = B, z = R). 

306. a = z cos <p-e p + pe* + L is the loop: {p = 
= sin <p, z = 1 }. 

Compute the circulation of the following vectors, given 
in spherical coordinates, along the given contours L, 
directly and with the aid of the Stokes theorem. 

307. a = r0e r -f-rsin 0 e„; L is the circle: jr = 1, 0 = 

-*}• 

308. a = r sin 0-e P -f L is the loop: |r = sin <p, 
0 = J, 0 <<p<n} . 

309. a = r<p0e„; L is a contour bounded by the semi¬ 
circle: jr = 7?, (p = -j-, O< 0 <ji| and its vertical 
diameter jcp = -£, 0 o|. 


Sec. 25. The Laplace operator 
in orthogonal coordinates 


If u = u (ft, ft, g 3 ) is a scalar function, then 

. 1 du . 1 du , 1 du 

« rad u =-ff7-^- e *+^-^r e2 +^-^- e 3- 


•=«i (? 1 . 72. 7s) e i + (ft, ft, ft) e 2 + a a (ft, ft, ft)e 8 , 
then 


»-httW[ i 

+ (2) 
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Using formulas (1) and (2), we obtain the following expres¬ 
sion for the Laplace operator A u: 

Au = div grad u = [ ±- ( 

d i du \ d lH_ i£, \ ~| 

8q t \ H t dq t ) + dq a \ H s dq 3 I J * 

In cylindrical coordinates, 

_ 1 d _ 

— p dp \ K dp I p* a<p* 

In spherical coordinates, 

A “ = l4Ee [y( r2si " e f-) +W ( !in e w) 

+ 4" (1ET ) ] = "T 'S' ( r * W) 

1 d*u 
r* sin* 0 dtp* ’ 

Example. Find all the solutions of the Laplace equation 
An = 0 that depend solely on the distance r. 

Solution. Writing the Laplace equation in spherical 
coordinates and taking into account the spherical sym¬ 
metry of the solution (it must not depend on 0 or <p), 
we have 

4 ““ 7 T i ('’£)= 0 <“ = “('». 


so that 

u = -y- + Cz* 

where C x and C 2 are constants. 

310. Given: a scalar field u = u ( M ) in cylindrical coor¬ 
dinates 

u (p, (p, z) = p 2 cp + z*<p 3 — p<pz. 

Find An. 
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311. Given: a scalar field u = u ( M) in spherical coor¬ 
dinates 

u (r, 0 , 9 ) = r *90 + r® 9 2 + 9 -f 0 2 . 

Find Au. 

312. Are the following functions harmonic? 

(1) u = p 2 cos 2<p. 

( 2 ) u = r cos 20 . 

313. Find all possible harmonic functions that 

(a) depend on 0 alone, 

(b) depend on 9 alone 

(in the spherical system of coordinates). 

314. Find all solutions of the Poisson equation 

Au = r "' 1 

in the spherical system of coordinates, provided u = u (r). 



ANSWERS 

f x = 2, 

1. (a) The half-line \ t 

[ y — — z, «/>0, 2^0 is traversed 

twice when — oo < t < -|-oo. (b) When t £(— oo, — 1) 

U (— 1. +oo) the point + (t+V* j 

twice traverses the half-line x + y = 1, 

. x 2 +y 2 =l,l (d) y=-f, *=-f. 

(e) x 2 +y 2 +z 2 =l, x — y = 0. 

7. 1 +k. 8. i + k. 9. -j + 4- k. 

10. -i + k. 11. ei-j + 2k. 12. No. 14. No. 

I 7 . (a) 2(£r); (b)||-| ! +(r, ; <o> [»,#■]. 

21. Circles lying in planes perpendicular to the vector a. 

22. The hodograph of velocity is a helical curve: x = 

= a cos t, y = a sin t, z = 2bt\ the hodograph of accele¬ 
ration is a circle: x = —a sin t, y = a cos t, z = 2b. 
2 fi _ da du . d»a _ d*a / rfu \ 2 da 

d< du dt ’ dt * — du* \ dt ) ' du dt % 

28. (t — l)e‘i + y (f — y sin j —arctan <k + c. 

29. 4- In (l+f 2 ) i + i-e‘*.j + sin tk + c. 

30. e sln 1 • i—|-sin * 2 -j + *k + c. 
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31. £i + (<cos<-sinl)j + -^k + c. 

32. j + nk. 33. (l-e-‘/ 2 )i + (<?i/ 2 -l)j + («-l) k. 
34. — ln2-j + k. 35. 2n 2 i + nj +n 2 k. 36. R = |si ^| ~. 
37. i? = A \t\ (1 + 9Z 2 ) 3 / 2 . 38. /? = 6 - 39. /? = -£-an. 

40. 7? = 2a cosh 2 *. 41. x + y = 0. 

42. x_„_K 5.-0. 43. f =i-. 44. i = 

45. +-^--|--^r = C, a family of triaxial ellipsoids. 

46. x* + y 2 — z = C, a family of paraboloids. 

47. x 2 + y 2 = Cz, a family of paraboloids. 

48. 2y 2 + 9z 2 = C, a family of elliptical cylinders. 

49. x + 2 y — z = C, a family of parallel planes. 

50. A family of planes resulting from the sheaf of planes 
a 3 x + a 2 y + a 3 z = C (b 3 x -f b 2 y + b 3 z) passing through 
the straight line 

a l x + a 2 y + a 3 z = 0, ) 
b l x + b 2 y + b 3 z = 0. j 

via the elimination of the straight line itself. Here, a lt a 2 , 
a 3 are the coordinates of the vector a; b lt b 3 , b 3 are the 
coordinates of the vector b. 

51. x 2 + y 2 + z 2 = C 2 , a family of concentric spheres. 


52. (a, b, r) = C or 


CL\ ^3 

6 , b 2 b 3 


= C, a family of parallel 


planes. 

53. 2 x — y = C, a family of parallel straight lines. 

54. y = Cx, C > 0, x =^0, a family of rays. 

55. y 2 = Cx is a family of parabolas with the vertex 
O ( 0 , 0 ) deleted. 

56. x 2 — i / 2 = C, a family of hyperbolas. 
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57. y = —x In C — C, C > 0, a family of straight lines. 

58. 59. 60. -*yU 3 . 61. -1. 


62. 

67. 


S-)/2. 63. 


2 /§ 


(|/2 + 3). 66. 0. 


- 2 . ( 


no* 


70. k. 71. <p = n. 72. 9 = 0 . 

73. cp = 0. 74. y = —x + 2ren, re 


(i + j-k). 


= 0 , ± 1 , ± 2 , 


75. x z + y z + z* = i. 78.79. a. 80. a(b, r) + b(a, r). 
81. 2|a| 2 r —2(a, r)a. 86. ~ = ^-. 


on to COS (r, 1) to n - , , 

87 - ~er= - W -; ir =0forrj.l. 

88 - l=7T. 89 - T- 1 - 

^=0 if grad u_Lgrad l>. 

91. (a) 1 in the direction of the y-axis; (b) 3 in the 
direction of the vector a= — 1—2j + 2k. 


92. y = e,x; z = e 2 x. 93. y = ^-x+C i \ z = y-x + C 2 . 

94. x*+y* + z2 = C\, 1 95. 1-1 = 1, T+2^" = 4 - 

x + y + z=>C 2 . I 

96* y f 2 = C 2 . 97# z=C\X , y = 0 2 . 98. xy—C|» 2 =^ 2 * 
99. x = C it 2y z —z 2 = C 2 . 100. 1-1 = C„ z = C 2 . 

101. j/ 2 + z 2 = C lt x = C 2 . 102. x = C t y, x = C 2 z. 


103. -£— = ) 

»0l 0 0 t I 

J^ = J^ +C3 ’ I 


12* 
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where b 01 , b 0i , b 03 are cooordinates of the vector b 0 . 
104. n = -3. 105. II = nR*y. 106. II = nR*h. 

107. n = 4ji/?»/(fl). 108. n=-y. 109. II = y. 

no. n=i-n R*h. in. n=nP. 

112 . n = y-it. 113. n = -J. 

114 . n = o. 115 . n = -5-. 116 . n = i-. 117 . n = 4n/?3. 

118. (a) n= -j; (b) 11 = -y; (c) I1 = n. 119. n = 0. 
120. n = 6nfl. 121. n = 0. 

122. n = n. 123. n = 0. 124. 11 = J-n (l. 

125. 

126. n = V2n. 127. II = 45ji. 128. II = -y^n. 

129. n = o. 

130. n = n. 131. i|>(r) = ^-. 132. 7r 4 . 133. 0. 134. 0. 

135. (z) = C — z, C = constant. 136. II = 4 nR 3 . 

137. divE = 0 (r =y&0). 

143. 16n. 144. jiH s. 145. -y-n. 146. 0. 147. y. 

148. 4n. 149. -y-n. 150. -y-n. 151. 2 R*. 

152. -y-n. 153. -1. 154. -n. 

155. Solenoidal field. 

156. Nonsolenoidal field. 

157. Solenoidal field. 

159. «p(r)=-^-, r^=0, C = constant. 

161. 162. ln-^. 163. . 

2 r, r, r, 


164. 0. 
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166. -jB 2 . 167. -y. 168. (a)-^-; (b)-J. 169. 0. 

170. |. 171. 3^3. 172.^-. 

173. — na 2 . 174. 1. 

175. — 2n. 176. 177. y 

179. — 2 (zi + xj -h t/k). 

180. 3 (z 2 —x 2 ) j. 181. (x+y) k. 191. co=lcurl v = --ij. 
193. f(x, z) = xz + x-f z + C, C = constant. 

195. 4n. 1%. —4 ji. 197. y. 198. -2«. 199. -y.. 

200 . 729n. 201 . 0 . 202 . -V 2 n. 

203. 2toJia 2 . Hint: v = [®, r). 

204. n c = l. 205. fi c = 3. 206. Dependent. 

207. Independent. 208. Dependent. 209. —1. 210.-^-. 
211. 0. 212. y. 213. -y. 214. y. 

216. -y. Hint : Supplement the path of integration L 
with the line segment OA of the x-axis. 

217. No. 218. Yes. 219. No. 220. Yes. 221. No. 222. No. 
223. Yes. 226. <p = x 2 yz. 227. 9 = x + xyz. 228. 9 = 
= spy — y 2 + xz. 229. 9 = ln|x + y+ z|. 230. 9 = 
= arctan (xyz). 231. 9 = r. 232. 9 = In r. 

233. 9 = y r 3 . 234. q> = ax + Py -\-yz + C, C — constant. 

235. 9 = xy + yz + zx + C. 236. 9 = xy + e z + C. 
237. 9 = e x sin y + z + C. 247. (a) Yes, (b) No, 
(c) Yes. 

249. u = C x x + C t . 

250. u = Ax 2 + Bxy -f Ay 2 , where A and B are arbi¬ 
trary. 

f 7 ~tt A-CiX + Cz if n^= 1 , 

251. u(x)={ 

l x In |x| -J-CjX + Cz if n = l(x^= 0 ). 
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252. 1 = ^- n. 253. /=--j. 254. / = -|ji7? 2 . 

255. 7 = 0. 

256. b = xj + (y - x) k. 257. b = (y 2 - 2xz ) k. 

258. b = (e x - xe u ) j. 259. b = + (2y* - 6 xz) k. 

260. b-* (x + y 2 ) i + (x 3 + y 3 ) k. 

261. b = — (xz 2 + yz«* 2 ) j — 2 xyzk. 

262. b=|(-zj+yk). 

263. b- 8 yzi + xzj + 7xyk. 

264. b = 2xy 2 zi — 3x 2 i/zj + x*y 2 k. 

265. b = ysinxz-i—^-sinxz-k. 

266. b = fq$*- k - 267 - (*) p-q>+c„ p = 2 + c 2 ; 

(b) p: - inC iq > ■ P = c 2 z ; ( c ) <P = Ci* r = C 2 sin 2 0 . 

268. grad u = 2 (p + cos<j>)e p — ^2 sin <p + -i- e* costp) e, 

— e 2 sin<p*e z . 

269. grad u=-(cos<p — 3 P In 3)e p -f (--j-sin 29 — sin <p) e 9 

+ sin 2 <pe z . 

270. grad u = 2rcos0-e r — rsinO-eg. 

271. grad u = (6rsin0-f-e r coscp —l)e r + 3rcos0-e e 



272. grad u=-fi(^-e r + ^-e e ). 

273. div a = 2 + -^-cos<p — e^sinz. 

274. div a = arctan p + —(z 2 + 2z)e l . 

275. div a = 4r — y cos 2 <p cot 0 + sin e . 

nni > , cos 20 i n a , asin® \ asin 0 

276. curl a = sin9 e r — ^2cos0+ rsin Q Je 9 --—e v . 
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277. curl a =--f cot0.e r + -f-c e + -^^e (p . 

278. curl a = — 2pe p + 31 ° <P e x . 

281. 24n. 282. jn. 283. An. 284. n/? 3 . 

285. 4ji/? 4 . 286. —|fl 3 . 

287. 48. Hint: Write the equations of the surfaces in 
spherical coordinates. 

288. u = p-\-tp-\-z-\-C. 

289. u = y(p 2 + <j>2 + z*) + C. 290. u = p<pz + C. 

291. u = eP sin <p + z* + C. 292. u = p<p cos z + C. 

293. u = rQ + C. 294. u = r* + <p + e + C. 

295. u = y (rep* + 0 2 ) + C. 2%. u = r cos <p sin 0 + C. 

297 . u = e T sin 0 + ln (1 + (p 2 ) + C. 

298. 1. 299. n 2 . 300. 2 nR. 301. n 2 . 302. 1. 

303. T + J T- 1 ’ 

304. 0. 305. — 2nH i . 306. n. 

307. n. 308. 0. 309. 0. 

310. Au = 4<p-^ + -^i + 2<p3. 

311. Au = 6<p0 + 12np 2 + A + 9 cot0 + -^- cot0 + -ii|le • 

312. (1) Yes, (2) No. 

313. (1) u(0) = C 1 ln|tan-|| + C 2 , (2) u (<p) = C,<p + C 2 . 

_Q—+^ + C 2 ,n ¥ ,-l, - 2 , 

314. u(r)= lnr + ^- + C 2 , n--l, (r^=0) 

_ --^- + ^- + C 2 , n=- 2. 



APPENDIX I 


BASIC OPERATIONS OF VECTOR 
ANALYSIS IN ORTHOGONAL 
CURVILINEAR COORDINATES 


1. The scalar field is given in orthogonal curvilinear 
coordinates, u = u (q lt q 2 , q 3 ). Then we have 
1 


grad u = 


H i dqi * H t dq t 


_e 2 + 


H 3 dq 3 


The Laplace operator is 

a.. 1 r 9 ( H t H 3 du \ , d / H x H 3 du \ 

~ H l H 3 H 3 L \ dqi ) + dq t \ H t dq 3 ) 

d ! HjH % du \ ] 

' dq 3 V H 3 dq 3 I J - 
Special cases : (a) The scalar field is given in cylindrical 
coordinates, u = u (p, q>, z). Then we have 
, du . 1 du , du 

en iu= w t l+T -^e,+ lr e,. 

The Laplace operator is 


Au_ p dp ( p dp ) + p» 3q>* + oz- 

(b) The scalar field is given in spherical coordinates, 
u = u (r, 0, 9 ). Then we have 

gr adu = ^e r + li±ee + 7i i 

The Laplace operator is 

= ( r 2 i) ~k {-W sin e ) 

+ 


1 d*u . d*u 


9 <?q> 



Appendices 


185 


2. The vector field is given in orthogonal curvilinear 
coordinates: 

a = ft (ft. 9a. 9s) ft + ft (ft, ft, 9s) ft 

+ ft (9i. 92. 9s) ft- 

Then we have 

J; _. . 1 T d(a,H t H t ) , d(a i H 1 H s ) , ] 

dlV a - H x H t H i L dq x + dq t + d 9i J • 


curl a = 


H,H S e ‘ H x H t ez H x H t e » 

a a a 

aq x dq t dq t 


| a x H\ a^H^ (I3H3 \ 

Special cases: (a) The vector field is given in cylindrical 
coordinates, 


a = «i (P, 9, z) e P + ft (P. 9. z) ft + ft (P. 9. z)e t . 
Then we have 


curl 1 


P 5 <P 
ft - 


dtp dz 

P ft ft 

(b) The vector field is given in spherical coordinates: 
a = ft (r, 0 , <p) e r + a t (r, 0 , 9 ) e e + a 3 (r, 0 , 9 ) e,. 
Then we have 

3(a,r») 1 d (a t sin 8) 1 3a, 

r dr "*■ rsiiiB 30 rsin0 a? ’ 


div a =- 


curl a = 


_3_ 

dr 



_3_ d_ 

30 39 

ra z a 3 r sin 0 
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AREA ELEMENTS 
OF COORDINATE SURFACES 


C nates" 

Coordinate surfaces 

Area elements 

General 

?1» «»■ ?3 

qi = C— constant 
q t = C = constant 
q 3 = C — constant 

dS l =H t (C,q t ,q 3 ) H a (C,q i ,q 3 )dq i dq 3 
dS 3 =H , (9„ C. g s ) H 3 (q, , C, ?,) d ?I dy 3 
dS 3 =Hi (9i, ? t , C) H t (9p 9*» Q <*9i 

Cylind¬ 

rical 

9i = P 

p = C = conslant 

dS = C dtf dz 

<?j = <P 

<p = C = constant 

dS = dpdz 

?, = * 

2 = C = constant 

dS = p dp dq> 

Spherical 



?i = r 

r = C = constant 

dS = C* sin 0 40 dip 

92 = 0 

0 = C = constant 

dS = r sin C dr dq> 

9s = <P 

q> = C = constant 

dS = r dr d0 
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compressible fluid 96 
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oordinate(s) 
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spherical 154 
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curl 108 
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of) 166 

orthogonal 154 
potential in 163 

definite integral 19 
del 130 
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of a vector function 14 
first and second (with respect 
to arc length) 27 
differential equations of vector 
lines 156 

directional derivative 39 
divergence in orthogonal coor¬ 
dinates 159 

divergence of a vector field 89, 90 


equation(s) 
continuity 
Maxwell 133 
Poisson 141 
wave 139 


field(s) 

harmonic 140 
irrotational 108 
Laplace 140 
potential 121 
flux 

in curvilinear coordinates 160 
of a vector 85 
of a vector field 58 
methods of computing 62 
force function 121 
formula(s) 
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Frenet 31, 33 
Gauss-Ostrogradsky 85 
Green’s 115, 117, 143, 144 
Frenet formulas 31, 33 
function 
force 121 
potential 121 
primitive 18 
vector 9 


Gauss-Ostrogradsky formula 85 
Gauss-Ostrogradsky theorem 85 
gradient in orthogonal coordi¬ 
nates 157 

gradient of a scalar field 44 
Green’s formulas 115, 117, 143, 
144 

Guldin theorem 86 

Hamiltonian operator 130 
harmonic field 140 
hodograph of a vector function 9 


indefinite integral 18 
integral 
definite 19 
indefinite 18 

line (see line integral) 97 
irrotational field 108 


Lame coelncients 155 
Laplace field 140 
Laplace operator 130, 135, 137 
Laplace operator in othogonal 
coordinates 174 
Laplacian 137 
level lines 35 
level surfaces 35 


limit of a vector function 11 
line integral 97 
calculating (in a vector field) 98 
computation of (in a potential 
field) 124, 125 
independence of path of inte¬ 
gration of 115 
properties of 98 

Maxwell equations 138 
method of 

introducing curvilinear coor¬ 
dinates 79 

projection on one plane 62 
projection on three planes 76 


normal, principal 28 
operator 

Hamiltonian 130 
Laplace 130, 135, 137 
orthogonal coordinates 
curl in 158 
divergence in 159 
gradient in 157 
Laplace operator in 174 
orthogonal curvilinear coordi¬ 
nates, basic operations in 184 
orthogonal system of curvilinear 
coordinates 154 
osculating plane 31 


Poisson equation 14t 
potential 121 

in curvilinear coordinates 163 
vector 146 
potential fields 121 
potential function 121 
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primitive function 18 isothermic 36 

principal normal 27, 28 level 35 


radius of curvature 29 
radius of torsion 32 
rotation 108 


scalar field 35, 37 
plane 38 

second-order differential opera¬ 
tions 130, 135, 136 
sink 90 

solenoidal fields 89, 95 
source 90 

star-shaped region 128 
Stokes' theorem 111 
surface(s) 


theorem 

Gauss-Ostrogradsky 85 
Guldin 86 
Stokes' 111 
torsion 31, 32 
radius of 32 

vector field 52 
flux of 58 
plane 56 
vector function 9 
integrating a 18 
vector line 52 
vector potential 146 

wave equation 139 
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